INTRODUCTION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

Autopolyploid genomes, containing many different homologous chromosomes that can pair and recombine in most or all possible combinations, have been underexplored at the molecular level due to the special problems they pose in genetic and molecular analyses (Sreenivasan et al. 1987; Burner 1997). The importance of autopolyploidy is highlighted by its prominence among cultivated crops, including sugarcane (8-18×), sugar beet (3×), ryegrass (4×), bermuda grass (3-4×), cassava (4×), potato (4×), alfalfa (4×), red clover (4×), Grande Naine banana (3×), apple cultivars (3×), and many ornamentals (Zeven 1979).

Among the world's leading crops with an annual production at a projected record 97 million metric tons in 1999/2000 (FAS Online, 1999), sugarcane is a classical example of a complex autopolyploid. Cultivated sugarcane varieties have ~80-140 chromosomes, comprising 8-18 copies of a basic × = 8 or × = 10 (D'Hont et al. 1995, 1998; Ha et al. 1999; Irvine 1999). Most chromosomes of cultivated sugarcane appear to be largely derived from Saccharum officinarum (Irvine 1999); however, in situ hybridization data suggest that about 10% may be derived from S. spontaneum (D'Hont et al. 1996). S. officinarum commonly has high sucrose content, low fiber content, thick stalks, little pubescence, rare flowering, and limited tillering. S. spontaneum does not accumulate sucrose, and is fibrous, thin-stalked, pubescent, profusely flowering, and abundantly tillering.

Like other vegetatively propagated plant species, cultivated sugarcane (Saccharum spp. hybrids) and its wild relatives are highly heterozygous. Pure inbred lines do not exist due to the difficulty of self-pollination and the random pairing of multiple homologous chromosomes. The segregating populations used in genetic studies are the progenies (first generation) derived from crosses between two cultivated varieties (Kang et al. 1983; Milligan et al. 1990) or cultivated varieties and wild species (Guimarães et al. 1997; Ming et al. 1998). Chromosome transmission is normal for most crosses, yielding n×+×n progeny (Burner 1997), but 2n×+×n transmission predominates in S. officinarum (2n = 80) × S. spontaneum F₁ and BC₁ crosses, a phenomenon known as "female restitution" (Bremer 1923; Price 1957).

The most abundant restriction fragment length polymorphisms in sugarcane are "single-dose restriction fragments" (SDRFs) showing 1:1 segregation in doubled haploid and interspecific F₁ populations (Wu et al. 1992; Da Silva et al. 1995; Guimarães et al. 1997; Ming et al. 1998). SDRFs represent 70% of the detectable polymorphic loci resulting from the segregation of alleles of different dosages (Da Silva 1993). The random chromosome pairing that is characteristic of autopolyploids makes it necessary to construct linkage maps for each parent of a cross, unlike diploid species, in which allelism permits a unified map of both parents to be generated. Studies using molecular markers have begun to resolve the genetic complexity of sugarcane by analysis of SDRFs (Da Silva et al. 1995; D'Hont et al. 1995; Grivet et al. 1996; Dufour et al. 1997; Guimarães et al. 1997; Ming et al. 1998). Comparative mapping has shown striking colinearity among the genomes of grasses (Hulbert et al. 1990; Ahn and Tanksley 1993; Ahn et al. 1993; Lin et al. 1995; Moore et al. 1995; Paterson et al. 1995), and even distantly related species (Paterson et al. 1996). A comparative approach has greatly expedited sugarcane genome analysis, in particular using the small diploid genome of closely related sorghum as a guide (Dufour et al. 1997; Guimarães et al. 1997; Ming et al. 1998).

Genetic tools for sugarcane have only recently become adequate to quantify the effect of many genomic regions on a trait. Two prior studies reported the association of DNA markers with disease resistance and flowering time in sugarcane. Daugrois et al. (1996) identified a putative major gene for rust resistance linked at 10 cM with a RFLP marker CDSR0029 in sugarcane cultivar `R570.' Guimarães et al. (1997) found an RFLP marker associated with short-day flowering. However, the mapping populations used in these two studies were too small (83 and 100 individuals, respectively) for comprehensive quantitative trait loci (QTL) analysis.

We report here the results of our study into the genetic basis of variation in sugar content among sugarcane genotypes using single-dose DNA markers. The fundamental complexity of autopolyploid genetics resulting from heterozygosity and lack of preferential pairing is further complicated by the fact that sugar content is a complex industrial trait influenced by variation in carbon fixation, photosynthate partitioning into sucrose, transportation and accumulation of sucrose (Berding et al. 1997; Moore et al. 1997) in harvestable biomass, and extractability of sucrose from biomass (Legendre and Henderson 1972). Our primary objectives were to determine the number and location of QTLs for sugar content in sugarcane, which is arguably the most important trait for the sugarcane industry; to investigate the molecular basis of phenotypic buffering that may contribute to the success of autopolyploid crops; and to investigate the possibility that candidate genes for QTLs affecting carbohydrate metabolism in biomass crops might be identified based on discrete mutations affecting seed development in other crops.

	RESULTS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

Sugar Content QTLs

GG × IND progeny values ranged from 39.4 to 249.0 lb sugar per ton of harvested biomass (lb/ton) (see Methods), a range that was ~40.1% wider than the (albeit large) difference between the parents (IND = 53.2, GG = 202.8). A full model comprised of 14 QTLs, eight from GG and six from IND, explained 65.5% of phenotypic variation (PV). The eight GG QTLs alone explained 38.6% of PV, while the six IND QTLs alone explained 36%. Among the 14 QTLs found, only one was inconsistent with the expected parental phenotypes, a GG QTL near pSB0279d that reduced sugar content by about 20.3 lb/ton. Also, a putative (P < 0.006) IND QTL near pSB0044d increased sugar content by 14.4 lb/ton. These two loci together account for part of the transgressive variation. A total of 18 putative QTLs (0.003 &lt P < 0.01) were found, and five (27.8%) of these were associated with significant QTLs on homologous chromosomes, or with candidate genes (Table 1, Fig. 1).

View larger version (61K): [in this window] [in a new window]   Figure 1   Comparative mapping of sugar content QTLs. Bars and whiskers indicate 1 and 2 LOD likelihood intervals. Solid lines connect homologous loci on different sugarcane and sorghum linkage groups. Arrows indicate the inferred locations of markers used to align the homologous linkage groups among maize, sorghum, and sugarcane based on published (Ming et al. 1998) and unpublished data (M. McMullen and E. Coe, pers. comm.). Individual sorghum linkage groups are represented by LGs A to J. Sugarcane linkage groups from four parental varieties are indicated by the last letter of the marker name: G (Green German); M (Muntok Java); I (IND 81-146); P (PIN 84-1). Approximate map positions of double-dose (#) markers are inferred by the method of Da Silva (1995). The letters in parentheses following the marker name represent the sorghum linkage groups where the marker mapped, if different from the corresponding location shown. Only regions that contain, or are homologous to, QTLs are shown. Twenty-eight of the 36 sugar content QTLs are shown on this figure, along with nine putative QTLs. The markers associated with the remaining eight sugar content QTLs could not be mapped in sorghum.

Sugar content of PIN × MJ progeny ranged from 49.0 to 65.2 lb/ton, a range about 30.4% wider than the difference between the parents (PIN = 37.2, MJ = 50.4). Negative sugar content values reflect the lower ratio of pol to brix. If this ratio is below 35% (varies slightly at different factories), the calculated sugar content value will be negative, indicating sucrose can not be separated from other soluble solids in cane juice (Pol value is explained in Methods). A full model comprised of 22 QTLs, 18 from MJ and four from PIN, explained 68.3% of PV (Table 1). The 18 MJ QTLs alone explained 45.7% of PV, while the four PIN QTLs alone explained 33.4% of PV. Allele effects of all QTLs were consistent with expected parental phenotypes, except one putative PIN QTL near pSB1368d increased sugar content by 9.4 lb/ton, accounting for part of the progeny transgression of parental phenotypes. Three DNA probes (CDSC0005, CDSC0042, CSU0428) were each diagnostic of two MJ QTLs at unlinked loci, and one DNA probe (CDSB0032) was diagnostic of two PIN QTLs at unlinked loci. A total of 23 putative QTLs (0.003 < P < 0.01) were found, and four (17.4%) were associated with significant QTLs on homologous chromosomes, or with candidate genes (Table 1, Fig. 1).

Some sugar content QTLs showed clear patterns of association with other traits. Sugar content was positively correlated with Pol (r = 0.64 in GG × IND, r = 0.92 in PIN × MJ) and stalk weight (r = 0.50 in GG × IND, r = 0.53 in PIN × MJ), and negatively correlated with ash content (r = 0.66 in GG × IND, r = 0.55 in PIN × MJ, Table 2). Most IND sugar content QTLs were associated with increased fiber and ash, and reduced stalk weight. Most MJ sugar content QTLs were associated with increased Pol and stalk weight, and reduced ash. Most PIN QTLs were associated with reduced Pol and increased ash.

Dosage Effects of Individual QTLs

In the four cases where single DNA probes detected sugar content QTLs at each of two or more unlinked loci, it was possible to investigate whether the dosage (zero, one, or two `copies') of the chromosomal region(s) containing the favorable allele(s) had nonadditive (i.e., nonlinear) effects on phenotype (Fig. 2), as described in the Methods section. All four showed nonlinear tendencies suggesting less-than-additive effects, but in only one case (CSU0428b, dM) did the regression line have a significant nonlinear (in this case, quadratic) component. Other traits for which significant effects were linked to larger numbers of loci detected by common probes provided a test of higher dosages. For example, two DNA probes each detected three loci associated with plant height in MJ, and another two DNA probes each detected four loci associated with plant height in MJ (R. Ming et al., in prep.). In all four cases, the regression lines showed less-than-additive gene action, with significant (P < 0.05) quadratic trends in three cases, and a significant quartic trend in one case.

View larger version (58K): [in this window] [in a new window]   Figure 2   Dosage effects of QTLs for sugar content and plant height (HT). Lowercase letters after the probe name represent the loci detected by this probe; the capital letter M or P represents parental variety MJ or PIN. Phenotypic effects of allele substitution at each locus are shown next to the locus (letter); the asterisk indicates the significant level at P < 0.01. Pl, Pq, Pc, and Pqt are the probabilities that linear, quadratic, cubic, and quartic components of the dosage curve are equal to zero.

Comparative Analysis of QTLs

Alignment with the high-density sorghum linkage map has enabled us to fill gaps in the sugarcane map (Ming et al. 1998), and also to evaluate correspondence of sugarcane QTLs to structural genes, phenotypic mutants, and QTLs affecting carbohydrate metabolism. Major genes or QTLs affecting maize seed carbohydrate status or levels of carbohydrate-metabolizing enzyme activities were placed on sorghum linkage groups (LGs) based on DNA markers mapped on both the 1998 maize linkage map (Davis et al. 1999), and our sorghum map (Ming et al. 1998). A total of 11 candidate genes or QTLs from maize in nine genomic regions were evaluated for correspondence to sugar content QTLs (Table 3). Four of these corresponded to at least six (and possibly as many as eight) sugar content QTLs in MJ, GG, and IND (none corresponded in PIN), and one corresponded to a putative MJ QTL. Among the 18 MJ QTLs, two corresponded to candidate genes, a level of correspondence that could be explained by chance in 94.6% of cases. Among the eight GG QTLs, at least three (possibly four) corresponded to candidate genes, a level of correspondence that could be explained by chance in only 4% of cases. Among the six IND QTLs, at least one (possibly two) corresponded to candidate genes, a level of correspondence that could be explained by chance in 33.5% of cases. The uncertainty in the number of GG and IND QTLs that correspond to candidate genes is related to the uncertain location of the probe CDSB0010 in the sorghum genome, where it could not be directly mapped due to lack of polymorphism. Its proximity to anchor markers that could be mapped in both sorghum and sugarcane suggest two possible locations, one near the bottom of LG D (with no candidate gene) and the other on LG G near the location corresponding to the su2 mutation of maize. CDSB0010 is being further investigated by physical mapping-if it should prove to map near su2, then the rates of correspondence of GG and IND QTLs to candidate genes would be explicable by chance in only 0.5% and 8.3% of cases, respectively.

Comparative QTL analyses were summarized in reference to sorghum linkage groups as follows:

Sorghum Linkage Group A

Sorghum Linkage Group B

Sorghum Linkage Group C

Sorghum Linkage Group D

Sorghum Linkage Group F

Sorghum Linkage Group G

Sorghum Linkage Group H

Sorghum Linkage Group I

Sorghum Linkage Group J

	DISCUSSION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

The autopolyploidy of sugarcane was reflected in a high level of apparent duplication of QTLs, reflected by both correspondence of QTLs from different genotypes and by segregation for QTLs at multiple, apparently homologous locations in individual genotypes. Across the 36 genomic regions that showed significant association with variation in sugar content in the four genotypes, single homologous genomic regions accounted for three QTLs in three cases (CDSR0067a-pSB0044 on LG A, CDSC0042-CDSR0029 on LG F, and pSB0106-UMC0114a on LG I), and two QTLs in five cases (CDSR0035 on LG C, UMC0044-CDSR0125 and CDSR0046-CSU0423 on LG D, CDO0202-CDSB0032 and SG0322-CSU0063b on LG G). In one of these cases (CDSR0067a-pSB0044), plus at least three additional cases (CSU0460-pSB0101 and SHO0038-CDSB0007 on LG B, CDSR0133 on LG J) putative QTLs that fell slightly below our stringent significance threshold also corresponded to significant QTLs or to each other.

The 36 sugar content QTLs (Table 1) correspond to only eight nonoverlapping regions of the sorghum genome (on LGs A, B, C, D, F, G, and I). This suggests that the observed QTLs may be accounted for by a much smaller number of ancestral genes that have been multiplied by the rapid duplication of chromosomes that has characterized sugarcane genome evolution since its divergence from a common ancestor shared with sorghum (Ming et al. 1998).

In six (75%) of these eight nonoverlapping regions, we find both QTLs from high sugar content parents (GG or MJ) that increase sugar content, and also QTLs from low sugar content parents (IND or PIN) that decrease sugar content. While such a result would simply reflect allelism in a diploid, this explanation is not adequate for an autopolyploid F1 population. The mapping of an interspecific F1 population means that the S. officinarum and S. spontaneum alleles have little chance to pair, as the map is based on heterozygosity and recombination that occurs in the S. officinarum and S. spontaneum parents, respectively (and therefore we made maps of each parent; see Ming et al. 1998). Therefore, the effect of each allele is estimated based on the average phenotype(s) of individuals that differ by the presence/absence of the allele from one parent, which is completely independent of the presence/absence of the allele from the other parent. The discovery of QTLs in the same genomic region(s) from S. officinarum (GG or MJ) that increase sugar content, and from S. spontaneum (IND or PIN) that decrease sugar content, provides independent confirmation of the importance of these genomic regions in the control of this trait. Confirmation from different varieties and/or species increases the level of confidence that a QTL exists in the region, and also suggests that DNA markers linked to the QTL may be useful in other germplasm.

One QTL (pSB0279dG) and two putative QTLs (pSB0044dI, pSB1368dP) showed phenotypic effects that were the opposite of what would be predicted based on the phenotypes of the parents contributing these alleles. The discovery of exceptional QTLs from low-sugar wild genotypes that increase sugar content, and QTLs from high-sugar cultivars that reduce sugar content, confer added incentive for incorporating marker-assisted selection into sugarcane breeding programs. Deleterious QTLs from the high-sugar parent could be purged, and favorable QTLs from exotic sources could be introgressed. The phenotypic effects of these unexpected alleles explain part of the observed transgressive segregation for sugar content.

Dosage Effects of Individual QTLs

Multiplex segregation at QT loci may be partly responsible for the phenotypic buffering that is one factor in the success of many autopolyploid crops, as reported in alfalfa for physiological measurements such as net CO₂ exchange, acetylene reduction, activities of ribulose-1,5-bisphosphate carboxylase, and leaf tissue concentrations of buffer-soluble protein, chlorophyll, and DNA (Leps et al. 1980; Pfeiffer et al. 1980; Meyers et al. 1982a,b; Molin et al. 1982).

In four cases, two or more loci detected by the same DNA probe were each associated with variation in sugar content, enabling us to investigate the possibility of such phenotypic buffering in sugarcane. Three DNA probes (CDSC0005, CDSC0042, CSU0428) were each diagnostic of two MJ QTLs at unlinked loci, and one DNA probe (CDSB0032) was diagnostic of two PIN QTLs at unlinked loci. Such associations may reflect multiplex segregation at orthologous genetic loci, or perhaps just coincidence of different QTL alleles at nearby loci, but in either case permit us to evaluate the net consequences of stacking multiple copies of a genomic region each associated with common phenotypic effects. "Stacking" of multiple doses of chromosomal segments containing favorable QTLs generally yielded diminishing effects on phenotype, especially in cases where high-order duplication could be tested (Fig. 2). This is similar to the results reported from stacking unlinked QTLs in a diploid, tomato, which were attributed to epistasis (Eshed and Zamir 1996). Epistasis in sugarcane is complicated by the possibility of nonlinear interactions between loci at homologous sites (such as we report), in addition to nonlinear interactions between unrelated loci (Eshed and Zamir 1996).

Detecting this type of phenotypic buffering provides strategic information for marker-assisted selection in autopolyploid crops. Although diagnostic DNA markers enable us to pyramid multiple QTLs in a polyploid, incorporating any one copy of the multiple alleles may obtain most of the desired effect in the breeding population.

Nonadditive gene action in multiple-dose QTLs may also confer evolutionary opportunities. If a single copy of a gene/QTL is physiologically sufficient, the extra copies are free to collect mutations, often becoming nonfunctional, but perhaps occasionally resulting in a distinctive new function which improves fitness.

An important future investigation regards the contribution of multilocus QTL genotypes to stability of performance across different environments. Sugar content is a trait of relatively high heritability (Kang et al. 1983); however, a role of multiple-dose QTLs in enhancing environmental stability would be of potentially great importance for less heritable traits.

It was curious that the highest-sugar genotype, Green German, showed only eight sugar content QTLsfar fewer than the 18 found in Muntok Java, which had much lower sugar content. The high ploidy of Green German (2n = 97-117) would make it possible, at least in principle, that additional favorable QTLs may be present in Green German but in so many doses that most progeny have several copiesand consequently phenotypic variation cannot be associated with marker segregation. This may suggest that our experiment has only detected a subset of the QTLs that are responsible for sugar contentspecifically, overlooking those QTLs that have large additive effects and have been driven to high frequencies by selection. This notion could be tested by crossing GG × IND progeny back to IND, and doing further QTL mapping. If this notion is true, then the tendency of multiple-dose QTLs for showing nonlinear dosage effects may be representative only of the subset of alleles that have not yet reached high frequency in improved sugarcane populations.

Candidate Genes for Sugar Content

The complexity of measuring sugar content, and large number of genes influencing the trait, suggest that it will be very difficult to identify the underlying genes using positional approaches alone. An overlapping genetic basis for variation in seed carbohydrate metabolism of grain crops and variation in stem carbohydrate accumulation of biomass crops would be a useful aid in the identification of candidate genes for QTLs affecting sugar content, a complex industrial trait. Perturbations in seed carbohydrate metabolism often result in discrete visible phenotypes, and several underlying genes have been cloned. A growing collection of maize mutants promises to provide additional candidate genes.

Several of the sugarcane QTLs we mapped correspond approximately to the genomic locations of maize mutants: two to sus1, and at least one (perhaps as many as three) to su2. Four additional sugarcane QTLs corresponded approximately to previously mapped maize QTLs that modify the activities of key sugar-metabolizing enzymes: two to ADP glucose pyrophosphorylase (Qagp3), one to sucrose phosphate synthase (Qsps1), and one putative QTL to Qsps3. The rice sps1 gene was also mapped on rice chromosome 1 and corresponded to maize chromosome bin 3.09 (Fig. 1; Sakamoto et al. 1995; Davis et al. 1999). A rice cDNA clone, R1966, which is homologous to the barley sus gene, was mapped on rice chromosome 6 and corresponded to maize chromosome bin 9.03 and sugarcane sus QTLs (Kurata et al. 1994). Because three copies of rice Sucrose Synthase (Wang et al. 1995) and two copies each of maize Sucrose Synthase (Causse et al. 1995b; Miller and Chourey 1995) and Sucrose Phosphate Synthase (McCarty et al. 1986; Binh et al. 1995) have been reported, it is possible that some additional sugar content QTLs in sugarcane might encode or modify the same enzyme but locate on different genomic regions. Even if the candidate genes identified to date prove helpful, other approaches will clearly also be necessary to reveal most (or even much) of the molecular basis for variation in the sugar content of sugarcane. Candidate genes were found for only 17.5% of the sugar content QTLs, and half of the candidates were themselves QTLs. Over all loci examined, the extent of association between sugar content QTLs and candidate genes was strong for Green German, suggestive but equivocal for IND, weak for MJ, and nonexistent for PIN. While sus1 and su2 are strongly implicated as candidates for a direct role in the genetic determination of sugar content, no candidates were found in several genomic regions that showed much stronger evidence of a role in sugar content based on numerous QTLs corresponding to these regions. For example, a sugarcane genomic region corresponding to the central region of sorghum LG A was associated with three QTLs, and one additional putative QTL, but no candidate genes have been found. Similarly, a region near CDSC0042 on sorghum LG F corresponds to four sugarcane QTLs but not to any candidate genes. The recent completion of a database of nearly 300,000 ESTs for sugarcane (http://sucest.lbi.dcc.unicamp.br/en/) and identification of genes involved in most steps of carbohydrate metabolism provide a valuable starting point.

	METHODS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

Mapping Populations

Two interspecific segregating populations were studied, each made by P. Tai, USDA-ARS, Canal Point, Florida. Due to the heterozygous nature of the sugarcane varieties and species, the progenies of these interspecific crosses were segregating like F₂ intercross populations in diploid species based on two populations investigated in this experiment. (1) 264 plants from S. officinarum `Green German' (GG, 2n = 97-117) × S. spontaneum `IND 81-146' (IND, 2n = 52-56) (GG × IND). (2) 239 F₁ plants from S. spontaneum `PIN 84-1' (PIN, 2n = 96) × S. officinarum `Muntok Java' (MJ, 2n = 140) (PIN × MJ). The chromosome numbers of a sampling of the progenies from these two crosses were 2n = 73-85 for GG × IND and 2n = 99-121 for PIN × MJ, indicating n + n transmission (Burner 1997).

Each individual plant of the progenies was vegetatively propagated through cuttings for three replications. Both populations were grown from November 1994 to February 1996, as randomized complete block designs with rows 1.5 m apart and plants 0.6 m apart in the row, at the Texas A&M Agricultural Research and Extension Center. As significant replication effects were found at frequencies lower than the Type I error rate, phenotypes of the three replications were averaged for analysis.

The parents of these crosses were chosen for their differences in sugar content, and also based on the inference that two were S. officinarum, and the other two were S. spontaneum (in order to maximize DNA polymorphism). The levels and patterns of DNA polymorphisms among the four genotypes generally supported this inference. However, determining the taxonomic affinity of sugarcane genotypes is often complex. For example, Green German and Muntok Java are listed as S. officinarum in the world catalog of sugarcane genetic stocks (Rao and Vijayalakshmi 1962). The chromosome numbers of original Green German and Muntok Java were 2n = 80 (Bremer 1923; Rao and Vijayalakshmi 1962). A recent investigation revealed the chromosome numbers of modern accessions of GG as 2n = 97-117 and MJ as 2n = 140 (Burner 1997). This raises the possibility that the clones of these two varieties used in the mapping project may be hybrids of unknown ancestry, although morphologically they resemble S. officinarum.

Phenotyping

Sugar content was expressed in the industry standard units of pounds of sugar per ton of cane, equivalent to the content of sucrose at 96% purity, calculated based on Brix and Pol values as described (Legendre and Henderson 1972). Brix is the percentage of all soluble solids, mostly sugars, minerals, and organic acids, in the sugarcane stalk. Pol is the level of sucrose in stalk juice determined by polarimetry; a "clarified" juice sample from which optically active nonsugar compounds have been removed (Birkett and Seip 1975) is placed in a standard optical cylinder and polarized light is passed through the cylinder. The degree of rotation of the plane of light exiting the tube is recorded. Sucrose and glucose are dextro-rotatory, while fructose is levo-rotatory. In sugarcane juice, glucose and fructose levels are usually similar and small, so cancel each other out. Percentages of phenotypic value were calculated from the range of phenotypic value in the segregating population divided by the difference of parental phenotypic value.

Genotyping

RFLP analysis used laboratory methods as previously described (Chittenden et al. 1994). DNA probes used for QTL mapping were selected based on preliminary analysis of 1255 single-dose RFLP markers (i.e., alleles segregating in simplex segregation ratios) for association with phenotypes in 85 plants (Ming et al. 1998); additional probes were picked at 20 cM or smaller intervals for a more comprehensive search of particular regions of the genome showing even tenuous associations with sugar content in the subpopulations. A total of 186 probes were mapped in both populations using methods described (Ming et al. 1998), and generated 243, 232, 122, and 138 single-dose markers for GG, IND, MJ, and PIN, respectively.

Data Analyses

Single-factor ANOVA was conducted (SAS/GLM, SAS Institute 1989) to determine the associations between RFLP markers and sugar content in sugarcane. Correlations among traits were calculated using SAS/CORR. When flanking markers were available, MAPMAKER/QTL version 1.1 was used to calculate LOD scores by interval mapping. Because single-dose markers of sugarcane were segregating in a 1:1 ratio, the same as a backcross population in diploid species, MAPMAKER/QTL for backcross populations was chosen for map construction and QTL analysis (Da Silva et al. 1995; Grivet et al. 1996; Guimarães et al. 1997; Ming et al. 1998). Significance thresholds of LOD gt; 2.5 (interval mapping) or P < 0.003 (analysis of variance) were used to declare QTLs, based on the genome size and marker density in our sugarcane maps (Lander and Botstein 1989). Markers associated with sugar content at P < 0.01 (LOD > 2.0) were deemed "putative QTLs," but shown on maps only when they corresponded to a genomic region that contained a significant QT locus or loci for sugar content, or a candidate gene in maize. The coefficient of determination R² was calculated using SAS/GLM for each marker or QTL, as the percentage of phenotypic variation explained by each marker or QTL. The total phenotypic variance explained was estimated by including all significant single- and multiple-dose QTLs in a full model for multiple regression analysis. The allele effect of each single-dose QTL was the average difference in phenotype of individuals differing by one copy of the indicated allele (single-dose versus zero-dose).

When two or more loci detected by the same probe were each associated with a trait at P < 0.01 by single-factor ANOVA, these loci were combined to investigate the effects of double-dose (0-2 copies), triple-dose (0-3) or quadruple-dose (0-4) genotypes for the genomic region(s) containing the favorable alleles using trend analysis (Gomez and Gomez 1984). The average phenotypic value of each class (0 copy, 1 copy, 2 copies, etc.) was used for testing the dosage effect. Linear, quadratic, cubic, and quartic trends for the effect of marker dosage on phenotype were tested using the CONTRAST statement of the SAS/GLM procedure. Strictly additive dosage effects would result in a significant (nonzero) linear trend, but nonsignificant higher-order trends. Significant higher-order trends reflect nonadditive QTL dosage effects. To declare correspondence between candidate genes and sugar content QTLs, the probability of these corresponding QTLs occurring by chance was calculated using the hypergeometric probability distribution as previously described (Lin et al. 1995; Paterson et al. 1995). Specifically, this used the following equation:
where n = the number of intervals which can be compared (defined as 30 cM, approximating a QTL likelihood interval); m = the number of `matches' declared; l = the total number of genes/QTLs found in the larger sample, and s = the number of genes/QTLs found in the smaller sample. The average genomic interval that correspondence between QTLs could be assessed was ~30 cM, and our sugarcane mapping data afforded comparisons across 1700 cM between maize, GG, and IND, and 1300 cM between maize, MJ, and PIN linkage maps. The number of intervals that could be compared were therefore estimated at 56 (= 1700/30) for GG and IND, and 43 (= 1300/30) for MJ and PIN. A match was declared between maize and sugarcane when corresponding intervals harbor candidate genes/QTLs in each. In the case of unlinked markers associated with sugar content QTLs, a match was declared when a marker associated with sugar content corresponded to the maximum-likelihood location of a maize sugar mutant/QTL within ~5 cM (since most QTLs spanned ~30 cM genomic region on average).