Effect of genotype and environment on branching in ... - Doust Lab

3 downloads 0 Views 603KB Size Report
processes determine the eventual shape of the plant and allow for dynamic reshaping ... plants, and strongly influence resource acquisition, com- petitive ability, and ..... or axillary branching depending on growing conditions. Factor. d.f.. MS. F.
Molecular Ecology (2006) 15, 1335–1349

doi: 10.1111/j.1365-294X.2005.02791.x

Effect of genotype and environment on branching in weedy green millet (Setaria viridis) and domesticated foxtail millet (Setaria italica) (Poaceae)

Blackwell Publishing Ltd

A N D R E W N . D O U S T and E L I Z A B E T H A . K E L L O G G University of Missouri-St Louis, Department of Biology, One University Boulevard, St Louis, MO 63121, USA

Abstract Many domesticated crops are derived from species whose life history includes weedy characteristics, such as the ability to vary branching patterns in response to environmental conditions. However, domesticated crop plants are characterized by less variable plant architecture, as well as by a general reduction in vegetative branching compared to their progenitor species. Here we examine weedy green millet and its domesticate foxtail millet that differ in the number of tillers (basal branches) and axillary branches along each tiller. Branch number in F2:3 progeny of a cross between the two species varies with genotype, planting density, and other environmental variables, with significant genotype–environment interactions (GEI). This is shown by a complex pattern of reaction norms and by variation in the pattern of significant quantitative trait loci (QTL) amongst trials. Individual and joint analyses of high and low density trials indicate that most QTL have significant GEI. Dominance and epistasis also explain some variation in branching. Likely candidate genes underlying the QTL (based on map position and phenotypic effect) include teosinte branched1 and barren stalk1. Phytochrome B, which has been found to affect response to shading in other plants, explains little or no variation. Much variation in branching is explained by QTL that do not have obvious candidate genes from maize or rice. Keywords: barren stalk1, branching, foxtail millet, genotype–environment interaction, QTL, teosinte branched1 Received 2 June 2005; revision accepted 4 October 2005

Introduction Plant architecture is a result of the initiation and differential elongation of growth axes (branches) (Bell 1991). These processes determine the eventual shape of the plant and allow for dynamic reshaping of its architecture under the influence of environmental stimuli such as shading, mechanical damage, and resource limitation (Thomas 2000). Production and differentiation of branches occur throughout the vegetative and reproductive phases of growth in most plants, and strongly influence resource acquisition, competitive ability, and reproductive success. In grasses, vegetative branches form at various levels above the ground (Gould & Shaw 1983; Clark & Fisher 1987). Most grasses form branches from the short basal nodes on the primary stem (culm). These are commonly Correspondence: Andrew N. Doust, Fax: 341-516-6233; E-mail: [email protected] © 2006 Blackwell Publishing Ltd

known as tillers and may produce adventitious roots. Some grasses also produce branches in the axils of leaves along the stem (cauline axillary branches, called here simply axillary branches). Often these axillary branches only elongate after the meristem of the main culm has transformed from a vegetative meristem to an inflorescence meristem (Perreta & Vegetti 2004). Further axillary branches can be initiated on previously initiated axillary branches, leading to a much-branched ‘bushy’ plant. The number and placement of the branches affect the distribution of the leaves and thus influence both light acquisition and shading of other plants. In many grasses, particularly annuals, all branches terminate in an inflorescence, so the number of branches that are initiated and grow may directly influence the number of seeds the plant produces. Branch production is controlled by a complex interplay of environmental inputs, hormonal responses, and genetic activity (McSteen & Leyser 2005). Light, particularly the red– far-red ratio, affects plant architecture via signals mediated

1336 A . N . D O U S T and E . A . K E L L O G G by the phytochromes, particularly phyB (Schmitt et al. 2003), and plant hormones such as auxin (Dharmasiri et al. 2005; McSteen & Leyser 2005). In maize, genes such as barren stalk1, knotted1, and teosinte branched1 affect axillary meristem formation, maintenance, and outgrowth, respectively (Sinha et al. 1993; Sinha & Hake 1994; Hubbard et al. 2002; Gallavotti et al. 2004). The signal transduction pathways that link all these components are beginning to be elucidated in a few model systems, including rice and maize (Takano et al. 2001; Sawers et al. 2002, 2004, 2005), but much remains to be learned about how environment and genetics interact in agricultural crops and weeds. The developmental plasticity shown by vegetative branching responses in weeds under different growing conditions is an important part of the weedy strategy (Perrins et al. 1992; Dekker 2003). Such genotype–environment interactions (GEI) have been found to be important in both vegetative and reproductive architecture in a number of species (Pigliucci 1998; Juenger et al. 2000; Schlichting 2002, 2003; Schmitt et al. 2003; Ungerer et al. 2003; Koornneef et al. 2004; Weinig & Schmitt 2004; Weinig 2005). In many weedy grasses, axillary branches continue to be initiated as long as favourable conditions for plant growth continue, allowing for continuous flowering and seed production (Dekker 2003). Domestication effectively reverses some of the branching characteristics of weeds, and in domesticated grasses such as maize or millet has led to a sharp reduction in vegetative branching to a few tillers and axillary branches (Doebley & Stec 1991; Doust et al. 2004). In our recent work, we have investigated differences in branching patterns between a serious weed of temperate croplands, Setaria viridis (green millet) (Dekker 2003), and its domesticated derivative, Setaria italica (foxtail millet) [treated as subspecies by several authors, including Wang et al. (1995) and Mabberley (1987)]. The sister relationship between these two species has been confirmed by chromosomal fluorescent in situ hybridization (FISH) and ribosomal, nuclear and chloroplast sequence data (Benabdelmouna et al. 2001; Doust & Kellogg 2002; Doust, Penly & Kellogg, unpublished). It is likely that foxtail millet was domesticated from a native ruderal variant of green millet that would have displayed variable branching, but that would have been under minimal selection for weedy responses to human farming activities. Thus, the differences seen between the two species today can be viewed as the result of divergence from a common stock under different selection pressures. The morphological changes seen between the two species include a shift from many tillers and axillary branches, each terminating in a short inflorescence with relatively few orders of branching (typical weed architecture), to few tillers and no axillary branches, with long inflorescences that have many orders of branching (typical domesticate architecture) (Harlan 1992; Doust & Kellogg 2002; Doust et al. 2004, 2005). The loss of axillary

branches in domesticated foxtail millet correlates with a decrease in the time over which inflorescences develop and mature, a characteristic that is important for efficient management and harvesting (Harlan 1992). This effectively reverses the broad dispersal period characteristic of the weedy progenitor. At the genotypic level, plasticity might be achieved by having multiple genes affecting a particular trait, with each gene responding to environmental inputs in a different way, or by having alternative suites of genes activated in different environments. Such genes can be identified using a quantitative trait loci (QTL) approach to locate regions of the genome that are differentially affected by the environment. Different combinations of alleles at these ‘environmentally responsive regions’ should lead to different reaction norms in the plants that bear them. Our previous studies (Doust et al. 2004, 2005) identified loci that were minimally affected by environment, and we interpreted these as the likely sites of selection during domestication. Here we return to some of the same data, but consider loci with a significant interaction between genotype and environment; these are important for the study of plasticity. In addition, we consider the possibility that ‘domestication’ loci might also be considered ‘weediness’ loci. In this interpretation, the S. italica allele at a locus may have been artificially selected for domestication, but the S. viridis allele at the same locus may have been selected for weediness.

Materials and methods Mapping, plant growth, and phenotype evaluation We used an F2 mapping population derived from a cross between foxtail millet (Setaria italica) and green millet (Setaria viridis) that had previously been used to construct a genetic map, using 257 restriction fragment length polymorphism (RFLP) probes from rice, foxtail millet, pearl millet and wheat (Devos et al. 1998; Wang et al. 1998). As detailed previously, we used 119 of these markers for our QTL analysis, chosen to cover the genome at approximately 10-cM intervals (Doust et al. 2004,2005). An additional 23 RFLP markers from maize plus 6 known genes were also added to the original map. We grew F3 offspring selfed from 120 of the original 127 F2 plants in four separate trials, trials 1 and 2 at high density with 5 representatives per family, and trials 3 and 4 at low density with 15 representatives per family (trials 3 and 4 were labelled ‘1’ and ‘2’ in Doust et al. 2004). Soil, fertilizer, water, and day length were standardized, but natural light intensity and average temperatures differed slightly among the trials, as trials 1 and 3 were grown in early summer and trials 2 and 4 in late summer. In the low density trials, replicates of families were planted one to a pot, and © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

E N V I R O N M E N T A L E F F E C T S O N B R A N C H I N G I N M I L L E T S 1337 randomized with respect to position in the greenhouse. In the high density trials, the five representatives of each family were planted in a single pot, and therefore were not randomized in position in the greenhouse, and are not true replicates. Thus for analyses that compared high and low density trials, we used least square mean values only. Means were also used for QTL analyses, because these are the best approximation to the phenotypes of the F2 plants used to construct the genetic map (Wang et al. 1998). Plants were harvested after the seeds had ripened. We counted the number of tillers (branches coming from the base of the plant) and the number of axillary branches (in the axils of culm leaves). In general these two classes of branching were easy to distinguish, although a few plants produced multiple axillary branches from new tillers late in the growing period. This seems to be a phenomenon separate from the production of tillers and axillary branches during ‘normal’ growth, but could not be distinguished in this analysis because all measurements were done at harvesting. Preliminary observations suggested that, during ‘normal’ growth, tillers are produced before flowering whereas axillary branches appear after the inflorescence terminating the main culm has been produced.

Data exploration We examined the shape of the frequency histograms for the trait values in each trial and found evidence of two distributions for most trait/trial combinations. In particular, there was a large excess of plants that had no axillary branches, resulting in an extreme positive skew to the data in all trials. As well, in the high density trials, many of the plants had only a single tiller. Previous observations (Doust et al. 2004) indicated that the foxtail millet parent probably lacks axillary meristems and therefore cannot produce axillary branches. Therefore, the large number of plants that do not have axillary branches may be the result of segregation of a gene responsible for this phenotype in the foxtail millet parent. However, lack of axillary branching may also be the effect of crowding (more families had no axillary branching in high density than in low density trials). It is also possible that some replicates in a family could have axillary branches while others lack them, because of incomplete segregation in these F 3 families. Therefore, to analyse quantitative variation in axillary branch number on only the subset of families that we could be confident had the potential to express such variation, we took low density trial 4 as our guide and pruned families that had means for axillary branch number that were less than one. This takes account of the possibility that some families were segregating for the lack of axillary branches, while keeping in the analysis families that had no axillary branches at high density © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

(possibly as the result of crowding) but had axillary branches at low density. After pruning, the data set for quantitative analysis consisted of 79 families. The frequency distributions of the pruned data set approached normality more closely than that of the full data set. Residuals obtained from a univariate analysis of each trait, using trial as the factor, indicated that a log transformation would equalize the variances amongst the four trials for axillary branch number. Tiller number in the low density trials was normally distributed, but we could not find any satisfactory transformation that would convert tiller number in the high density trials to a more normal distribution. Outliers were identified for each trait/trial combination by examining boxplots and searching for extreme values using the Descriptive Statistics– Explore command (SPSS 2005). Analyses were run with and without outliers, and differences noted where present.

Genetic correlations and genotype by environment interactions Genetic correlations were calculated for both the full and reduced data sets for the low density trials, but the small variances and extreme non-normality of trait distributions of tiller number in the high density trials precluded statistical analysis. Genetic correlations of tiller and logtransformed axillary branch number in each trial examine the extent to which the two traits might be under common genetic control (pleiotropy). Genetic correlations were also calculated for each trait between trials (e.g. tiller numbers in trial 3 vs. those in trial 4). These correlations indicate the extent to which the same set of genes underlies variation in that trait amongst trials. Genetic correlations were estimated as cov[trait1, trait2]/σtrait1σtrait2

(eqn 1)

where cov[trait1, trait2] is the covariance of the cross-product of the correlation between the two traits, and σtrait1 and σtrait2 are the square roots of the variance components for each trait (Robertson 1959; Falconer & Mackay 1996). Variance components were estimated using restricted maximum-likelihood estimation (SPSS 2005). Reaction norm plots were also assembled for the four trials, using all 120 families, to highlight both general trends between trials, as well as genotype–environment interactions.

Analysis of variance Two sets of analyses were performed. One set compared log-transformed axillary branch numbers (LOGAXB) between high and low densities and amongst the four trials using mean values for each family and the reduced data set

1338 A . N . D O U S T and E . A . K E L L O G G of 79 families. It was not possible to satisfy conditions for analysis of variance with the two high density tiller trials, so tiller number (TILL) was not analysed in this set of analyses. A nested GLM anova was used to partition variance into sources attributable to density [D], F2:3 family [F], trial nested within density [T(D)], and the interaction between these. The sources of variation in the model can be represented as y = µ + D + T(D) + F + D × F + T(D) × F + E

(eqn 2)

where µ represents the overall mean of the experiment, D is a fixed effect, and T(D), F, D × F, and T(D) × F are random effects. T(D) × F cannot be tested in this experimental design because the use of a single mean value for each family results in too few degrees of freedom. Thus, this source of variation was included in the error, E. Family and trial are treated as random effects because the families available for analysis are only a random selection of all possible genotypes that could be produced from this cross, and because the timing of the four trials was not planned. Because both trial and family are random factors in this analysis, the nested design does not allow an F-test to test the effect of density by itself (Quinn & Keough 2002, p. 318). Instead, a quasi F-test is presented in the nested analysis to test the effect of density, using a combination of the mean squares of all of the interaction terms as a denominator (Quinn & Keough 2002; SPSS 2005). A second set of analyses used data only from the low density trials, where values for individual replicates could be used. Here, the model can be represented as y=µ+T+F+T×F+E

(eqn 3)

Both TILL and LOGAXB were analysed for the low density trials. These two variables were also used as dependent variables in a multivariate anova. Analyses were done using the GLM anova programs in SPSS (SPSS 2005). Variance components were estimated using restricted maximum-likelihood estimation (SPSS 2005), and used to calculate the percentage of variance explained by each model effect. QTL analyses (detailed below) suggested three candidate loci; these were analysed for TILL and LOGAXB in the low density trials using mean values for each family. The three loci replace the family effect, and were treated as fixed effects. A full factorial model was fitted, using trial and the three gene loci and all appropriate interactions between them. LOGAXB was also analysed in both high and low density trials, using a similar nested model as in the first anova above. Appropriate interactions are all those for which all allele combinations were present (that is, to measure an interaction between gene1 and gene2, it is necessary to have present amongst the 79 families all possible

combinations of the three allelic states possible for each gene). After the first analysis, nonsignificant interactions were dropped from the model and the analyses re-run. These analyses had unequal numbers of samples in each group, so error mean squares for each of the terms were weighted to remove this bias (SPSS 2005).

QTL analyses Means for all 120 families were used in the QTL analyses, as we wanted to explore all possible genetic regions (QTL) underlying variation in phenotypic traits. The small number of families used makes it likely that not all QTL regions can be identified and that the size of the effect of the identified QTL regions may be exaggerated (Beavis 1998). Our previous report on tillering and axillary branching incorporated data only from the two low density trials analysed individually (Doust et al. 2004). Here we present new data on the high density trials and a joint analysis of all data on tillering and axillary branching. In trials 1, 2 and 4, the parental and hybrid ranges were examined for evidence of transgressive segregation; parents of the cross were not incorporated into the planting design in trial 3. QTL were detected using composite interval mapping (CIM), as implemented in QTL Cartographer (Basten et al. 2002). Background markers were selected at P = 0.05, and five background parameters were included as cofactors in each CIM model. Tests were made at 2.0-cM intervals, with a window size of 10 cM. QTL in separate trials were considered to be identical if their 1-LOD support intervals overlapped, and the sign of their additive effects was the same. Joint analysis of each trait for all four trials taken together was analysed using the module JZmapqtl (Basten et al. 2002). The joint analysis allows an estimate of GEI effects between the trial values for each trait, and provides a measure of both the main and interaction effects of detected QTL. Significance thresholds for QTL were calculated by 1000 permutations of the original data, using the same parameter settings as for the original analysis (Churchill & Doerge 1994; Doerge & Churchill 1996). We calculated significance levels both genome-wide and chromosome-wide (both at P < 0.05). Identification of QTL based on multiple chromosome-wide significance levels will increase type I error compared to the genome-wide level because nine different tests are being performed (one for each chromosome), but will also increase the probability of identifying more true QTL (Cheverud 2001). We considered occurrence of QTL significant at the chromosome-wide P < 0.05 level in more than one trial as evidence that a real QTL had been detected. The program epistacy (Holland 1998) was used to identify digenic epistatic interactions. There are (119 × 118)/2 © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

E N V I R O N M E N T A L E F F E C T S O N B R A N C H I N G I N M I L L E T S 1339 Fig. 1 Histograms of tiller number and axillary branch number in the four trials (using means from all 120 families). The continuous and dashed arrows show the positions of the means for foxtail and green millet, respectively, in trials 1, 2, and 4. HD, high density; LD, low density.

Fig. 2 Scatter plots of tiller number vs. axillary branch number in each of the four trials (using means from all 120 families). TILL, tiller number; AXB, axillary branch number; HD, high density; LD, low density.

comparisons that are made in this analysis, and a Bonferroni correction to the P < 0.05 experiment-wide significance level gives a per-comparison significance level of 7 × 10−5 (Rieseberg et al. 2003a). To identify possible candidate genes from maize and rice, we used markers mapped on maize and foxtail millet or rice and foxtail millet to define intervals on the maize or rice maps that correspond to QTL regions in foxtail millet. MaizeGDB (Lawrence et al. 2004) and Gramene (Ware et al. 2002a, b) were used to identify genes that had mutant phenotypes or putative functions that might affect some aspect of branching. Several of these were analysed by anova, as detailed above.

Results

millet consistently had more tillers than foxtail millet (Fig. 1). At high density, both green and foxtail millet produced fewer tillers than at low density. Our observations indicate that an increase in branching in either parent always increases the amount of seed produced, even if the inflorescences are somewhat smaller on laterproduced branches. Hybrids in trial 2 for tillering and trials 1, 2 and 4 for axillary branching showed transgressive segregation, where the trait values for the hybrids had values that exceeded the range between the two parents. However, no transgressive segregation was seen for tillering in trial 1 or 4, where the parental values were at the extremes of the range of values found for the hybrid population. In the F2:3 families both tiller number and axillary branch number were much lower at high density than at low density.

Phenotypic distribution of traits Both parents adjusted their number of tillers to planting density, but only green millet modified the number of axillary branches, producing fewer axillary branches at high than at low density (Fig. 1). Conversely, foxtail millet never produced axillary branches. In the three trials in which the two parents were grown (trials 1, 2 and 4), green © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

Genetic correlations The phenotypic relationships between tiller and axillary branch number in each trial indicate that there is a strong relationship in trial 3 but almost no relationship in the other three trials (Fig. 2). Genetic correlations between tiller and axillary branch number were examined for both

1340 A . N . D O U S T and E . A . K E L L O G G

Fig. 3 Scatter plots of relationships for each trait between all possible pairs of trials (using means from all 120 families).

the full data set (FDS) and pruned data set (PDS) in each of the two low density trials (3 and 4), and were significant for trial 3 (FDS r = 0.40, P < 0.001; PDS r = 0.38, P < 0.001) but not for trial 4 (FDS r = 0.10, NS; PDS r = 0.05, NS). This indicates that there is a significant overlap in genetic control of the two traits in one trial but not in the other. Scatter diagrams of each trait for pairs of trials indicate that there is a stronger relationship between trait values in the two low density trials than in the two high density trials (Fig. 3). Genetic correlations were not calculated for the high density trials but the relationship between tiller numbers in trials 3 and 4 was significant (FDS r = 0.72, P < 0.001; PDS r = 0.85, P < 0.001), as was that between logtransformed axillary branch numbers (FDS r = 0.64, P < 0.001; PDS r = 0.60, P < 0.001). This indicates that genetic control of variation in each trait is, to a large extent, similar between the two low density trials. To explore responses to environment, we examined variation in the reaction norms of families amongst trials (Fig. 4). A common trend in reaction norms is evident for tiller number, with low mean values at high density and higher mean values at low density. The reaction norms for

Fig. 4 Plots of reaction norms for means of traits across the four trials. All reaction norms are plotted as light grey lines. As examples of the range of reaction norms observed, several were arbitrarily chosen and highlighted with thicker black lines, together with standard error bars for each trial.

four exemplar families, chosen to show the range of variation in reaction norms, show the same general trend, although they also illustrate changes in both the magnitude of effects and in rank order (crossing of lines). The reaction norms for axillary branch number show a complex pattern of crossing reaction norms, with an especially large response in low density trial 3 (Fig. 4). © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

E N V I R O N M E N T A L E F F E C T S O N B R A N C H I N G I N M I L L E T S 1341

Factor

d.f.

MS

F

P value

VC

%V

D T(D) F D×F E

1 2 78 78 156

10.98† 4.38‡ 0.18§ 0.06¶ 0.047

2.50 92.89 3.19 1.17

NS *** *** NS

0.042 0.055 0.030 0.004 0.047

23.6 30.9 16.9 2.3 26.4

Table 1 Nested anova for density (D), trial (density) [T(D)], and family (F) and their associated interactions for LOGAXB (using family means)

VC, variance component; %V, percentage of total variance explained. Significance levels: *P < 0.05; **P < 0.01; ***P < 0.001. †Error term for MSD = MST(D) + MSD × F – MSE; ‡Error term for MST(D) = MSE; §Error term for MSF = MSD × F; ¶Error term for MSD × F = MSE.

Table 2 anova for trial (T), family (F) and trial by family (T × F) for TILL and LOGAXB (using replicates within families), for low density trials. Significance levels and abbreviations as for Table 1 Factor

d.f.

MS

F

Tiller number T 1 24.94† 5.42 F 78 15.61‡ 3.37 T×F 78 4.63§ 1.62 E 1771 2.85 Axillary branch number (log) T 1 103.77† 216.75 F 78 1.43‡ 2.96 T×F 78 0.48§ 2.32 E 1771 0.21

P value

VC

%V

* *** **

0.02 0.468 0.155 2.85

1.0 13.4 4.4 81.6

*** *** ***

0.11 0.039 0.023 0.208

29.0 10.3 6.1 54.7

†Error term for MST = 0.987 MST × F + 0.013 MSE; ‡Error term for MSF = MST × F; §Error term for MST × F = MSE.

Analysis of variance — density, trial and family A nested anova testing for the effect of density, trial(density) and family, using family means, on LOGAXB found significant effects for trial(density) and family, but no significance for density (Table 1). However, the percentage of the variance explained by density is large, and it is possible that the design of the experiment, with only two density levels (and thus only 1 degree of freedom for the comparison) did not allow for adequate testing of density responses. The two anovas performed on TILL and LOGAXB in the low density trials used all available replicates, giving greater power to separate error variation from variation due to trial, family or their interaction (Table 2). All factors and their interaction had significant effects on both TILL and LOGAXB. However, the percentage of variation explained by the factors is much smaller for TILL than for LOGAXB, with error variation accounting for over 80% of the variation in TILL. A manova was also performed, using TILL and LOGAXB as the dependent variables. This was significant for the multivariate Pillai’s Trace test as well as for the two univariate F-tests (SPSS 2005) (results not shown). © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

QTL analyses In previously published analyses on vegetative branching at low density, we reported QTL as worthy of interest if they were significant at the chromosome-wide level of P < 0.05 and were found in more than one trial (Doust et al. 2004, 2005). In the present study we are interested in both QTL with fixed differences between the species (QTL replicated across environments) and those that differ across environments (GEI). For this reason we have looked at QTL significant at the chromosome-wide P < 0.05 level in both individual trials and across all four trials in joint analyses. Individual analyses of the separate trials found varying numbers of QTL: for tiller number, 5, 10, 6, and 7 QTL were found in trials 1 to 4, respectively, of these 6 were also significant at the genome-wide P < 0.05 level (Table 3). These represent 16 QTL regions, as a number of QTL amongst the trials map to the same genomic region (Table 3, Fig. 5; note that, for clarity, only the stringent genome-wide significance line is shown in the figure). Analyses of the individual trials for log-transformed axillary branch number found 25 QTL significant at the chromosome P < 0.05 level (5, 8, 8, and 4 QTL in trials 1 to 4, respectively), of which 8 were also significant at the genome-wide P < 0.05 level. These represent 16 QTL regions (Table 3, Fig. 5). Segregation distortion is severe in one of these regions, on chromosome VIII, and moderate in another, on chromosome IX (Wang et al. 1998). This distortion may affect our ability to place QTL precisely in these regions. Several QTL regions identified for tillering and axillary branching appear to be controlling both traits, such as those at chromosome positions III-8, V-5, V-10, VI-6, and IX-14. The QTL region on chromosome VI is found for axillary branching in all four trials but is only found in high density trials for tillering. The QTL region on chromosome III is detected for tillering in trials 2 and 4, and for axillary branching in trials 1 and 3. These QTL may reflect multiple closely linked genes, each affecting one trait. Alternatively, these QTL might reflect the action of genes that are generally involved in branching, but that are detected in tillering or axillary branching depending on growing conditions.

1342 A . N . D O U S T and E . A . K E L L O G G Table 3 QTL (chromosome and closest marker to LOD peak) detected at chromosome-wide significance levels (P < 0.05) in individual trials. QTL in bold are also significant at the genome-wide level of P < 0.05, corresponding to loci shown as significant in Fig. 5. HD, high density; LD, low density; spring, planted late spring; summer, planted midsummer; TILL, tiller number; LOGAXB, log transformed axillary branch per tiller number; A , additive effect; D, dominance effect; R2, explained phenotypic effect (%). A positive sign for additive or dominance effects indicates an increase and a negative sign a decrease in the trait value

Trait TILL

Trial 1 HD Chrom A — II-4 — — — IV-7 V-5 — VI-6 VII-1

Spring D R2

+0.15 +0.02

+0.19 +0.14 +0.06 +0.16 −0.15 +0.07 −0.01 +0.12

Trial 2 HD Chrom A

I-10 7.6 — II-11 II-12 III-8

Summer D R2

+0.17 +0.13 +0.19 +0.01 −0.27 +0.07 +0.14 −0.17

12.2 — 8.2 — V-11 12.0 VI-8 4.9 VII-1

−0.07 +0.19 −0.43 +0.72

VII-9

−0.09 +0.38

IX-4 — IX-16

+0.22 +0.11

+0.08 +0.21

— — — — LOGAXB — — — — III-9 — V-5 V-9 — VI-7 — — — — IX-14

+0.27 −0.26

−0.48 +0.12 +0.38 +0.52 −0.55 −0.21

+0.51 +0.23

I-5 — — — 6.4 — — V-4 8.9 9.0 — V-14 16.6 VI-6 VII-4 VIII-2 VIII-11 — 5.1 IX-14

+0.43 −0.26 +0.48 −0.33

+0.38 +0.23

+0.14 −0.86 −0.37 −0.16 −0.37

+0.61 −0.08 +0.01 −0.45 −0.24

+0.54 +0.1

Joint analyses for each trait across the four individual trials found 8 significant QTL for tillering and 5 for axillary branching at the chromosome-wide P < 0.05 level. Of these, one QTL for tillering and one for axillary branching were also significant at the genome-wide level of P < 0.05 (Table 4, Fig. 5). In two cases, QTL for tillering and for axillary branching appear to be in the same genomic region (VI-7/8, IX-14), suggesting that individual QTL may have effects on more than one trait. Approximately half of the joint QTL for each trait had significant GEI effects, indicating that the effect of the QTL is significantly influenced by the environment. However, differences between trials observed in the individual analyses are not always discernable as GEI effects in the joint analyses.

Trial 3 LD Chrom A

7.2 I-9 — 6.5 — 6.9 — 7.6 III-11 IV-7 V-7 4.8 V-10 49.1 — — 5.9 — 14.7 IX-1 7.7 — 23.3 — 12.1 — II-1 II-7 II-9 III-8 — 12.1 V-5 — 8.2 — 26.0 VI-6 4.5 — 5.3 VIII-1 5.9 — IX-10 10.2 —

Spring D R2

+0.22 −0.27

+0.20 +0.56 +0.42 −0.70

−0.48 +0.31 +0.21 +0.24

−0.11 +0.37

+0.11 +0.12 −0.13 +0.13

+0.07 −0.08 −0.13 −0.02

+0.24 −0.01

−0.23 −0.02 −0.07 −0.11 +0.14 −0.07

Trial 4 LD Chrom A

Summer D R2

5.6 I-9 II-4 — — III-9 11.5 10.3 — 9.1 — 28.3 V-11 — —

+0.40 +0.22 +0.02 −0.43

11.6 5.5

+0.51 +0.12

15.1

−0.39 +0.16

9.2

VII-7

+0.31 +0.03

5.4

IX-3 IX-14 —

+0.21 +0.1 +0.46 +0.39

2.7 12.5

−0.22 +0.28

5.4

V-10 +0.32 +0.10 — −0.69 −0.14 17.1 VI-5 — 3.9 — VIII-10 −0.21 −0.14 8.1 — —

7.1

5.6

2.4 7.0 9.2 7.4

— — — — — IV-3 —

23.7

32.2

3.9

Nine QTL for tiller number and nine for axillary branch number were neither replicated between trials nor found in the joint analysis. Some of these may represent false positives but others are likely to represent loci that are affected by the environment, and thus provide a window into environmental plasticity. There is considerable variation in the pattern of LOD peaks between the individual trials (Fig. 5), and the pattern for the joint trials appears to incorporate many of the patterns seen amongst the individual trials. Figure 5 also shows that the trials differ substantially in how many LOD peaks approach significance, with the two high density trials having fewer QTL significant at the genome-wide P < 0.05 level. The amount of phenotypic variation explained by the QTL was generally high for all traits and trials. Individual QTL were © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

E N V I R O N M E N T A L E F F E C T S O N B R A N C H I N G I N M I L L E T S 1343 Table 4 QTL detected at the chromosome-wide level of P < 0.05 by the joint analysis of the four individual trials. QTL in bold are also significant at the genome-wide level of P < 0.05, corresponding to loci shown as significant in Fig. 5. GEI, genotype by environment interactions; *, significant at the chromosome-wide level of P < 0.05; NS, not significant; other abbreviations are as for Table 3 Trait

Chrom. position

Main effect

GEI

A

D

TILL

IV-2 IV-7 V-11 VI-8 VII-2 VIII-4 IX-14 IX-16 II-7 VI-7 VII-1 IX-11 IX-14

* * * * * * * * * * * * *

NS NS * * NS NS * * * * * NS NS

−0.01 +0.14 −0.06 +0.01 −0.02 −0.01 +0.06 +0.09 +0.12 −0.18 −0.09 +0.06 +0.16

−0.06 +0.07 +0.01 +0.07 −0.13 +0.10 −0.01 +0.06 −0.04 +0.04 +0.03 −0.02 +0.01

LOGAXB

only of moderate size, although one QTL for axillary branching on chromosome VI explained between 16% and 32% of the variation. QTL that were replicated across trials explained varying amounts of phenotypic variation in those trials, sometimes varying by a factor of two but sometimes, as for tillering on chromosome V-10/11, varying by as much as a factor of 6. This is not surprising, as the proportion of variation explained depends in part on the total variation and the amount explained by other loci, both of which differ among the trials. The additive effects of QTL for both traits were a mixture of positive and negative effects (Tables 3 and 4). Separate QTL with additive effects of differing sign for a particular trait indicate that each parent contains a mixture of alleles for that trait, some acting to change the phenotype towards that of one parent and others acting to change the phenotype towards the other. The additive effects of these QTL, if all alleles were of the same sign, are not sufficient to explain the range in variation seen in the F2:3 hybrids, indicating that dominance and epistatic interactions are important in generating phenotypic variation.

Trait

Trial

Locus1

Position

TILL TILL TILL TILL LOGAXB LOGAXB LOGAXB LOGAXB LOGAXB LOGAXB

1 1 1 1 1 1 1 3 3 3

RGC147 PSF470 PSM713 PSM713 PSM371 RGR1943 PSM671 RG83.6 PSF163 PSF63.1

II-9 I-4 V-4 V-4 III-12 VI-12 VI-6 VI-1 VIII-10 V-3

QTL?

Y Y

Y

Dominance effects are often as large as additive effects for both tillering and axillary branching, although in other cases the effects observed are below the power of this experimental design to detect unambiguously (Lander & Botstein 1989). The dominance effects for vegetative branching traits are substantially greater than those reported previously for inflorescence branching in this mapping population (Doust et al. 2005). Digenic epistasis was detected between some markers associated with QTL for tillering and axillary branching, but also between markers not associated with QTL (Table 5).

Comparative mapping There is good colinearity between the millet genome and that of other cereals (Devos et al. 1998, 2000), so that comparison of species is useful for suggesting candidate genes. We looked particularly for genes that were in the identified QTL regions and that were known to control branching in other species. One well-studied gene that controls outgrowth of tillers and axillary branches in maize is teosinte branched1 (tb1), which has been associated with QTL controlling vegetative branching in maize × teosinte crosses, and which has been hypothesized to suppress axillary meristem elongation (Doebley et al. 1997; Hubbard et al. 2002). We hybridized a maize cDNA clone of tb1 to the F2 mapping filters for millet, and placed the gene on Setaria chromosome IX, in the same region as several QTL for branching. Another gene which is associated with QTL on chromosome V is barrenstalk1 (ba1). We confirmed this placement by hybridization of a foxtail millet DNA clone of this gene, and placed it between markers psm768 and rgc385. This corresponds to QTL regions found for tillering in both individual trials 2, 3, and 4 as well as the joint analysis. We had initially suspected that phytochrome B would be a candidate for one of our QTL, as this has been shown to control changes in branching and stem elongation that are associated with plant crowding. However, only a single minor QTL was found for tillering in trial 2 in the region on chromosome IX where phytochrome B maps, accounting for only 6% of the variation.

Locus2

Position

RGC950 RGC389 RGC601.2 RGR642 RGR1943 RGR830 PSM713 RGC950 UGT737 PSM671

VIII-12 VII-1 V-12 VIII-3 VI-12 III-13 V-4 VIII-12 VI-4 VI-6

© 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

QTL?

J, Y J J

Y

Y

R2 0.32 0.25 0.26 0.26 0.41 0.34 0.16 0.35 0.24 0.16

Table 5 Digenic interactions calculated using epistacy (Holland 1998). Results only reported if the overall probability of interaction is P < 7 × 10−5 (Bonferroni correction to give experiment-wide error rate of P < 0.05). Locus1 and locus2 are the two markers for which an interaction is being tested; QTL? recognizes those loci that are associated with a QTL for that trait and trial combination; Y, presence of QTL in individual analyses; J, presence of QTL in joint analysis; R2, proportion of variation explained by the digenic interaction. Trait abbreviations as for Table 3

1344 A . N . D O U S T and E . A . K E L L O G G

Fig. 5 Plots of LOD scores at 10-cM intervals across the nine chromosomes of the genome in the four individual analyses and in the joint analysis of tiller number and log-transformed axillary branch number. The plot is constructed as if the chromosomes have been laid end to end. The LOD score curve in the joint analysis composed of closed diamonds represents the main effects while the curve made of open circles represents GEI. The straight dashed line in trials 1 to 4 indicates the genome-wide significance level of P < 0.05, as do the open and compact dashed lines in the joint analyses. In the joint analyses, the significance of the main effect is given by the compact dashed line and that of the GEI effect by the open dashed line. N.B. This is a more stringent significance cut-off than the chromosome-wide significance level applied in Tables 3,4.

The QTL region on chromosome VI, which accounts for up to 30% of the variation in axillary branch number, has no obvious candidate genes from other grass species.

Analysis of variance — density, trial, and the three major QTL We used tb1 on chromosome I, ba1 on chromosome V, and hhu33, a marker closely associated with the large QTL on chromosome VI, to represent the three major QTL for tillering and axillary branching.

A nested anova, testing for the effect of density, trial (density), and the three markers and their interactions on LOGAXB found that density was not a significant factor, but that trial(density) and the three QTL markers were significant. Density × ba1 and trial(density) × hhu33 were also significant. Interestingly, as already noted above, density explains a large proportion of the variation, but has only one degree of freedom and thus fails to achieve significance. Trial(density) and hhu33 also explain large proportions of the variation, suggesting that these factors are important in determining LOGAXB variation (Table 6). © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

E N V I R O N M E N T A L E F F E C T S O N B R A N C H I N G I N M I L L E T S 1345 Table 6 Nested anova for density (D), trial(density) [T(D)], ba1, tb1, hhu33 and selected interactions for LOGAXB (using family means). Significance levels and abbreviations as for Table 1

Table 8 anova for trial (T), ba1, tb1, hhu33 and significant interactions for LOGAXB (using family means), in the low density trials. Significance levels and abbreviations as for Table 1

Factor

d.f.

MS

F

P value

VC

%V

Factor

d.f.

MS

F

P value

VC

%V

D T(D) ba1 tb1 hhu33 D × ba1 T(D) × hhu33 E

1 2 2 2 2 2 6 278

7.59† 2.60‡ 0.48§ 0.64§ 1.88¶ 0.24§ 0.15§ 0.06

3.34 18.67 8.27 10.95 12.57 4.06 2.66

NS ** *** *** ** * *

0.055 0.051 0.004 0.009 0.026 0.004 0.005 0.058

25.9 24.1 1.9 4.2 12.3 1.9 2.4 27.4

T ba1 tb1 hhu33 T × hhu33 E

1 2 2 2 2 136

5.18† 0.67‡ 0.31‡ 0.83§ 0.32‡ 0.05

19.04 13.85 6.38 2.73 6.55

* *** ** NS **

0.094 0.016 0.006 0.017 0.012 0.051

48.0 8.2 3.1 8.7 6.1 26.0

†Error term for MSD = 0.87 MST(D) + 0.13 MSE; ‡Error term for MST(D) = 0.84 MShhu33 × T(D) + 0.16 MSE; §Error term for MSba1, MStb1, MSD × ba1, MShhu33 × T(D) = MSE; ¶Error term for MShhu33 = 0.95 MShhu33 × T(D) + 0.05 MSE.

Table 7 anova for trial (T), ba1, tb1, hhu33 and significant interactions for TILL (using family means), in the low density trials. Significance levels and abbreviations as for Table 1 Factor

d.f.

MS

F

P value

VC

%V

T ba1 tb1 hhu33 ba1 × tb1 E

1 2 2 2 4 134

2.34† 5.21† 1.93† 0.03† 2.23† 0.72

3.27 7.25 2.69 0.05 3.11

NS ** NS NS *

0.022 0.095 0.006 0.000‡ 0.120 0.707

2.3 10.0 6.3 0 12.6 74.4

†Error term = MSE. ‡VC estimate set to zero, because term is redundant.

Analysis of TILL in the low density trials shows significant effects for ba1, and ba1 × tb1, but no significant effects for tb1, hhu33, trial or other interactions (Table 7). Analysis of LOGAXB in the low density trials differs from the analysis including high density trials in that hhu33 is no longer significant, although trial and trial × hhu33 are still significant. Ba1 and tb1 are also still significant (Table 8). Hhu33 is not significant, even though it explains as much variation as ba1 and tb1. This is because of the inclusion of the interaction term between trial and hhu33, which is used as the error term to test the significance of hhu33. Almost 50% of the variation in LOGAXB is explained by variation between the two (low density) trials.

Discussion Phenotype of the weed vs. the domesticate Our data show that the weed green millet is more responsive to environmental differences than its domesticated © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

†Error term MST = 0.84 MST × hhu33 + 0.16 MSE; ‡Error term MSba1, MStb1, MST × hhu33 = MSE; §Error term MShhu33 = 0.95 MST × hhu33 + 0.05 MSE.

descendant foxtail millet. At low density, green millet produces more than three times as many axillary branches as it does at high density, whereas foxtail never produces axillary branches. Although both species produce more tillers at low density, the increased number in green millet is more than double the increase observed in foxtail millet. These results confirm the observations of Dekker (2003) that green millet has highly plastic responses to shading and resource availability, whereas foxtail millet exhibits less plasticity to environmental conditions. Many of the environmentally labile changes seen in green millet should have strong fitness benefits, especially those increasing the duration of flowering and thus the amount of seed produced. In contrast, the lack of sensitivity to environmental changes seen in foxtail millet is likely the result of human selection for ease of cultivation and harvesting. Human selection for a reduction in axillary branching in foxtail millet is evident in the absence of axillary branches at either planting density. This, together with our failure to find any visible axillary branch primordia in the leaf axils of foxtail millet, even with the aid of a microscope, suggests that axillary branch primordia are either initiated but arrested at a very early stage or that they are not initiated at all. These two hypotheses suggest different candidate genes: in maize ba1 controls axillary meristem initiation while tb1 controls the elongation of the axillary branches. The pattern of expression of orthologous copies of these genes in foxtail and green millet is currently under study.

Plasticity and its genetic basis We only calculated genetic correlations in the low density trials (3 and 4), because of the extreme non-normality of tiller distribution at high densities. Tiller and axillary branch number had significant genetic correlations in trial 3 but not in trial 4. This, coupled with the results from the QTL analyses and our observations on the timing of tiller and axillary branch production, suggested that the two

1346 A . N . D O U S T and E . A . K E L L O G G forms of branching are actually different phenomena that might be expected to be under separate genetic control. This has also been reported in pearl millet (Pennisetum glaucum), a species that is more closely related to foxtail millet than either maize or rice (Poncet et al. 2000). Tillers are initiated early in the life cycle of the plant, whereas axillary branches only elongate after the inflorescence terminating the main culm has initiated. We also observed that some genotypes only produced axillary branches on the main culm, whereas other genotypes produced axillary branches on other axillary branches, leading to multiple orders of branching. Although we made no quantitative measurements of overall plant size, there appeared to be no relationship between plant size and number of branches in any particular trial. Rather, genotypes with multiple orders of branching could be either relatively small or large compared to other genotypes with few or no branches. Thus there is potential for considerable phenotypic plasticity in branching, dependent on particular allele combinations. Genetic correlations between tiller numbers in the two low density trials were high, as were those for axillary branch numbers. This indicates that the same suite of genes underlies the phenotypes in at least these two low density trials. This is not surprising for TILL, but is more so for LOGAXB, because of the large differences seen in the reaction norms between these two trials. The QTL analysis provides further insight into the correlations observed. The correlations between trials presumably are reflected in the 13 QTL detected by the joint analysis (Table 4). However, many QTL were found in individual trials but not in the joint analysis, and these may reflect the activity of additional loci that are significantly affected by environment. Environmental effects are also evident in the QTL identified by the joint analysis, as indicated by significant GEI, as well as by their detection in only a subset of trials. These ‘environmentally responsive regions’ of the genome are found on I-9/10, II-4, II-7, III-8/ 9, IV-7, V-4/5/7, V-9/10/11, VI-6/8, VII-1, VIII-1/2, VIII10/11, IX-3/4, IX-14, and IX-16 (Tables 3 and 4). We observed transgressive segregation in trials 1, 2 and 4, in which some hybrid families had mean tiller number and/or axillary branch number exceeding the range seen between the two parents. Not all trait/trial combinations showed transgressive segregation, so that, for example, the green millet parent in trial 4 actually had more tillers than any of the hybrids. Both tillering and axillary branching in all trials had QTL with both positive and negative additive effects, the favourable recombination of which could lead to transgressive segregation. Transgressive segregation is commonly found in QTL studies (Rieseberg et al. 2003b), and our data are no exception. Dominance and epistatic effects also contributed to the observed variation in some trials, indicating that they too varied between environments. Dominance effects in most

cases were substantial and could account for differences in the hybrid populations over and above those attributable to the combination of additive effects. Significant epistatic effects were also found in trial 1 for both tillering and axillary branching, and in trial 3 for axillary branching. These effects were found both between QTL positions identified in the individual and joint analyses, as well as between marker positions outside the QTL regions. This raises the intriguing possibility that some regions of the genome do not by themselves contribute to changes in phenotype, but can do so in specific combinations with other regions. It will be necessary to construct RILs with combinations of these particular regions to test this hypothesis. We show that the F2:3 families have different reaction norms, as predicted for any species likely to disperse into heterogeneous environments (Sultan & Spencer 2002). Each family could have a distinct combination of alleles for each of the 14 ‘environmentally responsive regions’, which would immediately lead to high variability in the magnitude — and possibly the direction — of response. We found that several of the QTL for tillering and axillary branching occurred in the same genomic region, and it may be that these regions contain genes that pleiotropically affect both types of branching. Two of these regions, found in individual and joint analyses for both tillering and axillary branching, are the lower region of chromosome IX (near markers 14 and 16 and the mapped location of teosinte branched1), and the region near the centromere on chromosome V (near markers 10 and 11 and the mapped location of barren stalk1). The anovas provide their own insight into the cause of variability in the hybrid population. We deliberately removed a portion of the families from the anovas in order to improve the tractability of the data, and to remove what appeared to be a subpopulation of families within the trial that never branched. This was presumably because of the inheritance of a homozygous factor from the foxtail millet parent, which itself never has axillary branches. Even so, there was still an extreme positive skew to the distribution of tiller values in the high density trials, which is likely the result of severe suppression of branching in crowded conditions. Difference in planting density is obviously not the only environmental factor at work in this experiment, because differences between trials at the same density were often as large as those between different densities. This is consistent with the unexpected result from the individual QTL trials that the pattern of QTL at either density differs substantially between trials, and that trials at one density may share replicated QTL regions with trials from the other planting density. These observations are supported by the anova results, where the nested analysis showed that axillary branching was affected much more by differences between trials and between genotypes than by density. © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

E N V I R O N M E N T A L E F F E C T S O N B R A N C H I N G I N M I L L E T S 1347 Differences in growing times may help explain differences between trials at the same density. By chance, each density level was grown in early summer one year and in late summer the following year. The different growing periods were characterized by differences in overall temperature and in light intensity, although artificial illumination maintained constant day length. It is not possible to disentangle differences in growing period from other unmeasured differences between trials, because of the lack of replication of densities at each growing time.

Candidate genes It is challenging to identify genes responsible for the phenotypic effects discussed above, especially so in a group of plants that lacks the genetic resources of model systems such as maize, rice, and Arabidopsis (CTC 2003; Phillips 2005). Our approach has been comparative, using the rice genome and abundant developmental mutants of maize (Doust et al. 2004, 2005) to flag likely candidate genes that can then be tested one by one. Any candidate gene in this study had to have a mutant phenotype in a model organism that was similar to some aspect of the phenotypic differences between foxtail and green millet, and had to be localized to a QTL region implicated in control of that phenotype. Two candidate genes that had relevant phenotypes and co-localized with QTL were teosinte branched1 (tb1) and barren stalk1 (ba1). Tb1 maps to the bottom of chromosome IX and is a transcription factor that suppresses axillary meristem elongation in maize relative to its wild ancestor, teosinte (Doebley et al. 1991,1995; Doebley 1995; Lukens & Doebley 1999; Hubbard et al. 2002). Its promoter region has been shown to have been under selection during maize domestication (Wang et al. 1999; Clark et al. 2004), and is a good candidate for involvement in domestication of foxtail millet. One joint and two individual QTL for tillering map to the region of tb1; that in trial 2 explains approximately 20% of the variation while that in trial 4 explains approximately 10%. Two QTL in trials 1 and 2 for axillary branching are also in this region. Another possible candidate gene from maize that may be involved in vegetative branching is barren stalk1 (ba1). The mutant phenotype of ba1 causes reduction or elimination of lateral branches throughout the plant (Neuffer et al. 1997; Gallavotti et al. 2004). Ba1 maps to a position of chromosome V that is covered by QTL for tillering found in three individual trials as well as the joint analysis, and by QTL for axillary branching found in two individual trials. Several replicated QTL fall in regions with no obvious candidate gene from maize, including QTL on chromosome III and VI. The QTL for axillary branching on chromosome VI is found in all individual trials and in the joint analysis, and is one of the largest and most consistently © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

found QTL. The failure to identify candidate genes in this region may be due to lack of synteny between millet and maize in those regions, or to genes underlying these QTL that have not been identified by forward genetics in maize, or to genes that are novel in foxtail millet. This QTL is unusual in that the additive and dominant effects are negative, i.e. they should act to decrease axillary branching. This QTL might be responsible for stopping initiation of axillary branch meristems, and thus acts as a ‘switch’ that is overlain by quantitative variation in branching controlled by other genes. We used the marker hhu33 as representative of this QTL (it is the closest marker to the peak of the QTL likelihood score on that chromosome, and is within 3 cM of the peak). anovas using tb1 and ba1, and hhu33, found contrasting patterns between tillering and axillary branching, as well as between axillary branching analysed over low density or over both high and low density trials. Differences between QTL and anovas can be attributed to the different data sets used for each (with or without plants that lack the ability to produce axillary branches). In the low density trials, anovas on tillering found significant effects for ba1 and ba1 × tb1. These results match emerging evidence from maize that tb1 and ba1 are both involved in control of tiller number (Lukens et al. 1999; Gallavotti et al. 2004). In maize, genetic experiments using NILs that are homozygous for either the QTL containing tb1 (on maize chromosome 1) or that containing ba1 (on maize chromosome 3) have demonstrated epistasis between the two regions, with the combined effect of the two regions having more phenotypic effect than the addition of the separate effects of each region would have (Lukens et al. 1999). Genetic experiments using single and double mutants have also shown that ba1 is completely epistatic to tb1 (Ritter et al. 2002; Gallavotti et al. 2004). This makes sense functionally as ba1 appears to be necessary to establish the population of meristematic cells that will be the new axillary meristem (Gallavotti et al. 2004), whereas tb1 regulates elongation of axillary shoots (Hubbard et al. 2002). Failure of ba1 expression will result in the absence of axillary meristems, and thus no opportunity for tb1 to act. Tb1 is regarded as the most important gene affecting tillering in maize (Doebley et al. 1995; Hubbard et al. 2002), but, due to differences in the initiation of axillary meristems between millet and maize, it appears that ba1 may be as important as tb1 in controlling tiller number in millets. anovas of axillary branch number also found that trial was by far the most important factor explaining variation, accounting for almost 50% of the variation in the low density analysis. Ba1 also has a significant effect, as does, to a lesser extent tb1. Hhu33 accounts for more variation than tb1, but is nonsignificant, because the F-test uses the interaction term between trial and hhu33 as the error term. However, the interaction between trial and hhu33 is

1348 A . N . D O U S T and E . A . K E L L O G G still significant. In the combined high and low density analysis, the three gene loci are significant, along with interactions between density and ba1, and trial(density) and hhu33. The results indicate that gene effects, both singly and in combination with other genes and environmental factors, explain a significant proportion of the variation in branching. The major effect of differences between trials is also clear. To attempt to identify other possible genes, we examined sections of the rice genome that were presumed orthologous to regions of the millet genome covered by QTL. Many transcription factors were identified, but, as the phenotypic effect of these genes is unknown, they cannot as yet provide good candidate genes for the present study. Rice genomic regions colinear to QTL regions on chromosomes V, VI and VII also include a variety of auxin and gibberellin pathway mediators (Doust et al. 2004). This approach offers considerable promise in identifying novel candidate genes, and will be significantly improved by increasingly comprehensive annotation of the rice genome. As this happens, we expect further candidate genes from rice and maize to be identified. Our results indicate that the two components of vegetative branching in millet grasses — tillering and axillary branching — are under largely separate genetic control, and respond differently to density and other environmental effects. QTL analyses revealed that each trait was under polygenic control, and that different trials uncovered different combinations of significant QTL. Several of these QTL regions contain candidate genes whose position, mutant phenotype, and statistical significance suggest that they may be involved in the control of branching in millets. These include tb1 and ba1. Our approach, utilizing quantitative genetic, QTL, and comparative mapping, gives new insight into the control of vegetative branching, a key element of the weediness strategy in grasses such as green millet, and a focus of reverse selection in the conversion of green millet into the domesticated foxtail millet.

Acknowledgements We thank Katrien Devos and James Beales for access to the mapping populations and for help with mapping maize RFLP markers onto the foxtail millet map, Kathy Upton for her patience in growing and maintaining thousands of plants, John Doebley for the clone of tb1, Andrea Gallavotti and Anthony Verboom for unpublished primers for ba1 and kn1, respectively, Bob Marquis for discussions on statistics, Torbert Rocheford and Jim Cherevud for advice throughout this project, and the Kellogg Laboratory Group for helpful comments. This project was supported by National Science Foundation grant MCB-0110809 to E.A.K.

References Basten CJ, Weir BS, Zeng Z-B (2002) QTLCartographer: A Reference Manual and Tutorial for QTL Mapping. North Carolina State University, Raleigh, North Carolina.

Beavis WD (1998) QTL analyses: power, precision, and accuracy. In: Molecular Dissection of Complex Traits (ed. Paterson AH), pp. 145–162. CRC Press, Boca Raton, Florida. Bell AD (1991) Plant Form. Oxford University Press, Oxford, UK. Benabdelmouna A, Abirached-Darmency M, Darmency H (2001) Phylogenetic and genomic relationships in Setaria italica and its close relatives based on the molecular diversity and chromosomal organization of 5S and 18S-5.8S−25S rDNA genes. Theoretical and Applied Genetics, 103, 668–677. Cheverud JM (2001) A simple correction for multiple comparisons in interval mapping genome scans. Heredity, 87, 52–58. Churchill GA, Doerge RW (1994) Empirical threshold values for quantitative trait mapping. Genetics, 138, 963–971. Clark LG, Fisher JB (1987) Vegetative morphology of grasses: shoots and roots. In: Grass Systematics and Evolution (eds Soderstrom TR, Hilu KW, Campbell CS, Barkworth ME), pp. 37–45. Smithsonian Institution Press, Washington, D.C. Clark RM, Linton E, Messing J, Doebley JF (2004) Pattern of diversity in the genomic region near the maize domestication gene tb1. Proceedings of the National Academy of Sciences, USA, 101, 700–707. Complex Trait Consortium (CTC) (2003) The nature and identification of quantitative trait loci: a community’s view. Nature Reviews Genetics, 4, 911–916. Dekker J (2003) The foxtail (Setaria) species-group. Weed Science, 51, 641–656. Devos KM, Pittaway TS, Reynolds A, Gale MD (2000) Comparative mapping reveals a complex relationship between the pearl millet genome and those of foxtail millet and rice. Theoretical and Applied Genetics, 100, 190–198. Devos KM, Wang ZM, Beales J, Sasaki T, Gale MD (1998) Comparative genetic maps of foxtail millet (Setaria italica) and rice (Oryza sativa). Theoretical and Applied Genetics, 96, 63–68. Dharmasiri N, Dharmasiri S, Estelle M (2005) The F-box protein TIR1 is an auxin receptor. Nature, 435, 441–445. Doebley J (1995) Genetics, development, and the morphological evolution of maize. In: Experimental and Molecular Approaches to Plant Biosystematics (eds Hoch PC, Stephenson AG), pp. 57–70. Missouri Botanical Garden, St Louis, Missouri. Doebley J, Stec A (1991) Genetic analysis of the morphological differences between maize and teosinte. Genetics, 129, 285–295. Doebley J, Stec A, Gustus C (1995) Teosinte branched1 and the origin of maize: evidence for epistasis and the evolution of dominance. Genetics, 141, 333–346. Doebley J, Stec A, Hubbard L (1997) The evolution of apical dominance in maize. Nature, 386, 485–488. Doerge RW, Churchill GA (1996) Permutation tests for multiple loci affecting a quantitative character. Genetics, 142, 285 – 294. Doust AN, Devos KM, Gadberry MD, Gale MD, Kellogg EA (2004) Genetic control of branching in foxtail millet. Proceedings of the National Academy of Sciences, USA, 101, 9045–9050. Doust AN, Devos KM, Gadberry MD, Gale MD, Kellogg EA (2005) The genetic basis for inflorescence variation between foxtail and green millet (Poaceae). Genetics, 169, 1659–1672. Doust AN, Kellogg EA (2002) Inflorescence diversification in the panicoid ‘bristle grass’ clade (Paniceae, Poaceae): evidence from molecular phylogenies and developmental morphology. American Journal of Botany, 89, 1203–1222. Falconer DS, Mackay TFC (1996) Introduction to Quantitative Genetics, 4th edn. Pearson Educational Limited, Harlow, Essex, UK. Gallavotti A, Zhao Q, Kyozuka J et al. (2004) The role of barren stalk1 in the architecture of maize. Nature, 432, 630–635. © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

E N V I R O N M E N T A L E F F E C T S O N B R A N C H I N G I N M I L L E T S 1349 Gould FW, Shaw RB (1983) Grass Systematics, 2nd edn. Texas A&M University Press, College Station, Texas. Harlan JR (1992) Crops and Man, 2nd edn. American Society of Agronomy, and Crop Science Society of America, Madison, Wisconsin. Holland JB (1998) epistacy: a SAS program for detecting twolocus epistatic interactions using genetic marker information. Journal of Heredity, 89, 374 –375. Hubbard L, McSteen P, Doebley J, Hake S (2002) Expression patterns and mutant phenotype of teosinte branched1 correlate with growth suppression in maize and teosinte. Genetics, 162, 1927–1935. Juenger T, Purugganan M, Mackay TFC (2000) Quantitative trait loci for floral morphology in Arabidopsis thaliana. Genetics, 156, 1379–1392. Koornneef M, Alonso-Blanco C, Vreugdenhil D (2004) Naturally occurring genetic variation in Arabidopsis thaliana. Annual Review of Plant Biology, 55, 141–172. Lander ES, Botstein D (1989) Mapping Mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics, 121, 185–199. Lawrence CJ, Dong Q, Polacco ML, Seigfried TE, Brendel V (2004) MaizeGDB, the community database for maize genetics and genomics. Nucleic Acids Research, 32, D393–D397. Lukens LN, Doebley J (1999) Epistatic and environmental interactions for quantitative trait loci involved in maize evolution. Genetical Research, Cambridge, 74, 291–302. Mabberley DJ (1987) The Plant-Book. University of Cambridge, Cambridge, UK. McSteen P, Leyser O (2005) Shoot branching. Annual Review of Plant Biology, 56, 353 –374. Neuffer MG, Coe EH, Wessler SR (1997) Mutants of Maize. Cold Spring Harbor Press, Cold Spring Harbor, New York. Perreta MG, Vegetti AC (2004) Structure and development of the branching system in Melica macra (Poaceae). Flora, 199, 36–41. Perrins J, Williamson M, Fitter A (1992) Do annual weeds have predictable characters? Acta Oecologica, 13, 517–533. Phillips PC (2005) Testing hypotheses regarding the genetics of adaptation. Genetica, 123, 15 –24. Pigliucci M (1998) Developmental phenotypic plasticity: where internal programming meets the external environment. Current Opinions in Plant Biology, 1, 87–91. Poncet V, Lamy F, Devos KM et al. (2000) Genetic control of domestication traits in pearl millet (Pennisetum glaucum L., Poaceae). Theoretical and Applied Genetics, 100, 147–159. Quinn GP, Keough MJ (2002) Experimental Design and Data Analysis for Biologists. Cambridge University Press, Cambridge, UK. Rieseberg LH, Raymond O, Rosenthal DM et al. (2003a) Major ecological transitions in wild sunflowers facilitated by hybridization. Science, 301, 1211–1216. Rieseberg LH, Widmer A, Arntz AM, Burke JM (2003b) The genetic architecture necessary for transgressive segregation is common in both natural and domesticated populations. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, 358, 1141–1147. Ritter MK, Padilla CM, Schmidt RJ (2002) The maize mutant barren stalk1 is defective in axillary meristem development. American Journal of Botany, 89, 203–210. Robertson A (1959) The sampling variance of the genetic correlation coefficient. Biometrics, 15, 469– 485. Sawers RJH, Linley PJ, Farmer PR et al. (2002) Elongated mesocotyl1, a phytochrome-deficient mutant of maize. Plant Physiology, 130, 155–163. Sawers RJH, Linley PJ, Gutierrez-Marcos JF et al. (2004) The Elm1 © 2006 Blackwell Publishing Ltd, Molecular Ecology, 15, 1335–1349

(ZmHy2) gene of maize encodes a phytochromobilin synthase. Plant Physiology, 136, 2771–2781. Sawers RJH, Sheehan MJ, Brutnell TP (2005) Cereal phytochromes: targets of selection, targets for manipulation? Trends in Plant Science, 10, 138–143. Schlichting CD (2002) Phenotypic plasticity in plants. Plant Species Biology, 17, 85–88. Schlichting CD (2003) Origins of differentiation via phenotypic plasticity. Evolution and Development, 5, 98–105. Schmitt J, Stinchcombe JR, Heschel MS, Huber H (2003) The adaptive evolution of plasticity: phytochrome-mediated shade avoidance responses. Integrative and Comparative Biology, 43, 459– 469. Sinha N, Hake S (1994) The Knotted leaf blade is a mosaic of blade, sheath, and auricle identities. Developmental Genetics, 15, 401– 414. Sinha NR, Williams RE, Hake S (1993) Overexpression of the maize homeobox gene, KNOTTED-1, causes a switch from determinate to indeterminate cell fates. Genes and Development, 7, 787–795. SPSS (2005) Statistical Package for the Social Sciences (Windows), Version 13. SPSS, Chicago. Sultan SE, Spencer HG (2002) Metapopulation structure favors plasticity. American Naturalist, 160, 271–283. Takano M, Kanegae H, Shinomura T et al. (2001) Isolation and characterization of rice phytochrome A mutants. Plant Cell, 13, 521–534. Thomas P (2000) Trees: Their Natural History. Cambridge University Press, Cambridge, UK. Ungerer MC, Halldorsdottir SS, Purugganan MD, Mackay TFC (2003) Genotype–environment interactions as quantitative trait loci affecting inflorescence development in Arabidopsis thaliana. Genetics, 165, 353–365. Wang ZM, Devos KM, Liu CJ, Wang RQ, Gale MD (1998) Construction of RFLP-based maps of foxtail millet, Setaria italica (L.) P. Beauv. Theoretical and Applied Genetics, 96, 31–36. Wang R-L, Stec A, Hey J, Lukens L, Doebley J (1999) The limits of selection during maize domestication. Nature, 398, 236–239. Wang R-L, Wendel JF, Dekker JH (1995) Weedy adaptation in Setaria spp. I. Isozyme analysis of genetic diversity and population genetic structure in Setaria viridis. American Journal of Botany, 82, 308–317. Ware D, Jaiswal P, Ni J et al. (2002a) Gramene: a resource for comparative grass genomics. Nucleic Acids Research, 30, 103–105. Ware DH, Jaiswal P, Ni J et al. (2002b) Gramene, a tool for grass genomics. Plant Physiology, 130, 1606–1613. Weinig C (2005) Rapid evolutionary responses to selection in heterogeneous environments among agricultural and nonagricultural weeds. International Journal of Plant Sciences, 166, 641–647. Weinig C, Schmitt J (2004) Environmental effects on the expression of quantitative trait loci and implications for phenotypic evolution. Bioscience, 54, 627–635.

Andrew Doust studies the evolution of morphological diversity in plants, focusing on phylogenetic patterns and underlying developmental and genetic processes. The work reported here is part of an intergrated analysis of phylogeny, development and genetics in millet grasses. Elizabeth Kellogg studies the evolution of development in plants. Her work focuses on the cereal grasses and their wild relatives, and integrates data from molecular phylogenetics, developmental morphology, quantitative genetics, and comparative gene expression.