|Home | About | Journals | Submit | Contact Us | Français|
When plants become shaded by neighbouring plants, they perceive a decrease in the red/far-red (R/FR) ratio of the light environment, which provides an early and unambiguous warning of the presence of competing vegetation. The mechanistic bases of the natural genetic variation in response to shade signals remain largely unknown. This study demonstrates that a wide range of genetic variation for hypocotyl elongation in response to an FR pulse at the end of day (EOD), a light signal that simulates natural shade, exists between Arabidopsis accessions. A quantitative trait locus (QTL) mapping analysis was done in the Bayreuth×Shahdara recombinant inbred line population. EODINDEX1 is the most significant QTL identified in response to EOD. The Shahdara alleles at EODINDEX1 caused a reduced response to shade as a consequence of an impaired hypocotyl inhibition under white light, and an accelerated leaf movement rhythm, which correlated positively with the pattern of circadian expression of clock genes such as PRR7 and PRR9. Genetic and quantitative complementation analyses demonstrated that ELF3 is the most likely candidate gene underlying natural variation at EODINDEX1. In conclusion, ELF3 is proposed as a component of the shade avoidance signalling pathway responsible for the phenotypic differences between Arabidopsis populations in relation to adaptation in a changing light environment.
Plants growing in stands with increasing vegetation are exposed to a reduction in the red to far-red (R/FR) ratio of the incoming light, caused by the enrichment of FR light reflected by neighbouring plants. With increasing plant density, plants also experience a significant reduction of photosynthetically active radiation (PAR), resulting from the depletion of red and blue light absorbed by leaves. Changes in the R/FR ratio are finely sensed primarily by phyB and secondarily by phyD and phyE, and the changes in blue light photons of PAR are sensed by cryptochromes and phototropins (Devlin et al., 1998, 1999; Pierik et al., 2004). All these photoreceptors act coordinately to trigger important adaptive responses known collectively as the shade avoidance syndrome (SAS). The SAS includes increased growth rate of stems and petioles and an acceleration of flowering that enables plants to anticipate light competition by neighbouring plants (Franklin, 2008).
A significant progress in the understanding of SAS molecular mechanisms has been achieved in the last decade through the analysis of mutants with impaired shade avoidance responses. In plants subject to a low R/FR ratio, the active form of phyB (Pfr) is converted to the inactive form (Pr) and exits the nucleus. The reduction of phyB Pfr increases the nuclear abundance of PIF4 and PIF5, two proteins that interact with the Pfr form of the phytochromes, enhancing the expression of genes involved in elongation growth (Nozue et al., 2007; Lorrain et al., 2008). Furthermore, recent studies suggest that gibberellin (GA) and phytochrome signalling pathways converge to regulate elongation growth (de Lucas et al., 2008; Feng et al., 2008). In plants exposed to a low R/FR ratio or to a depletion of blue light, the promotion of hypocotyl and petiole growth is partially accounted for by the degradation of DELLAs, a family of growth-repressing proteins acting in the GA signalling pathway (Djakovic-Petrovic et al., 2007).
Although induced mutagenesis continues to be an important approach for gene discovery (Ostergaard and Yanofsky, 2004), the study of the molecular bases of natural variation has an enormous potential to overcome several limitations of mutagenesis, such as no obvious phenotypic variation, unstable phenotypes, narrow genetic background, and a limited range of phenotypic variation (Tonsor et al., 2004). The analysis of important life history traits such as seed dormancy and flowering time led to the cloning of some genes responsible for the natural variation of these traits (Johanson et al., 2000; Michaels and Amasino, 2000; Sheldon et al., 2000; Bentsink et al., 2006). Albeit phyB is the primary phytochrome mediating SAS (Nagatani et al., 1991; Somers et al., 1991), and numerous allelic polymorphisms exist in its amino acid sequence (Filiault et al., 2008), the molecular bases underlying phenotypic variation in shade avoidance responses in natural populations remain largely ignored.
A huge natural genetic variation has been reported in SAS. Studying the hypocotyl growth and flowering time responses to shade in >100 accessions of Arabidopis thaliana, Botto and Smith (2002) found that both physiological responses induced by the same FR signal are not correlated, suggesting that shade avoidance signalling pathways have evolved separately. More recently, Adams et al. (2009) demonstrated that the insensitivity of flowering time in the Bla-6 accession to a low R/FR ratio is due to the high expression level of the floral repressor FLOWERING LOCUS C (FLC), conferred by a combination of functional FLC and FRIGIDA (FRI) alleles, with a weak FY allele regulating the expression of the floral integrator FLOWERING LOCUS T (FT). By quantitative trait locus (QTL) mapping analysis, it has been shown that the timing of flowering is regulated by intricate interactions between environmental signals and genetic factors (Alonso-Blanco et al., 2009). When plants of three recombinant inbred populations were exposed to a combination of seasonal and vegetation shade environments, it was demonstrated that a major proportion of loci were mapped in specific environments, and only a reduced proportion of them were pleiotropic to different conditions (Botto and Coluccio, 2007).
In this study QTLs controlling hypocotyl growth under white light (WL) and WL+EOD (end of day far-red light) were mapped. The EODINDEX1 QTL was identified on chromosome 2, which explained >30% of the relative growth variation to EOD segregating in recombinant inbred lines (RILs) from the Bayreuth (Bay) and Shahdara (Sha) accessions of A. thaliana. Sha alleles at EODINDEX1 were shown to impair growth inhibition under WL, and accelerate the leaf movement circadian rhythms. A fine mapping and complementation analysis strongly suggest that ELF3 is the most likely candidate gene for EODINDEX1. It is hypothesized that ELF3, a component of the circadian clock involved in the shade avoidance signalling, is a gene subject to selective pressure that provides an adaptive advantage in shade environments.
The core population of 164 RILs derived from crosses between Bay and Sha (Loudet et al., 2002), originated in opposite altitudinal places, was used. This population is available at the Arabidopsis Biological Resource Center (ABRC; Ohio State University, Columbus, OH, USA). The RIL population was phenotyped in two different light treatments to map QTLs affecting hypocotyl elongation. Due to failure in germination, the final number of lines tested was 141. Twenty seeds of each genotype were sown in clear plastic boxes (20cm long×10cm wide×1.5cm tall) containing 100ml of 0.8% agar and incubated in darkness at 5°C for 4d. Chilled seeds were exposed to a 1h red (R) pulse, and transferred to darkness for 24h to induce homogeneous seed germination at 23°C. Then the boxes with seedlings were placed in a WL (70μmol m−2 s−1) chamber under short-day conditions (8/16h light/dark; SD) at 23°C. Two light treatments were established: WL and WL+EOD. The EOD consisted of a 30min pulse of FR light (40μmol m−2 s−1), provided by incandescent lamps in combination with an RG9 filter (Schott, Germany). The experiments lasted 4d. For each box (a replicate), the 10 longest seedlings were measured with a ruler. At least three replicates were measured for each genotype and treatment, unless otherwise indicated in the text.
Accessions used for natural genetic variation for shade avoidance response on hypocotyl growth experiments were Est-0, St-0, Ct-1, Ba-1, Ll-2, Mr-0, Lm-2, Pog-0, Oy-0, Gr-1, Mh-0, Di-1, Ms-0, Kas-1, Ra-0, Kin-0, Van-0, Pi-0, Ge-1, Mt-0, Bur-0, Bay-0, Sha, Col-0, Ler, Cvi, Co-4, Bla-1, Bla-6, Sf-2, and mutant lines phyB-1 and phyB-9, all obtained from the ABRC.
Marker segregation data for the Bay×Sha RIL population were obtained from Loudet et al. (2002). Thirty-eight markers that cover the five chromosomes with an average genetic distance of 10.8cM between markers were used. MAPMAKER/EXP 3.0 was used to construct the linkage map (Lander et al., 1987). Linkage groups were verified with a minimum LOD=3 and a maximum distance=50cM (Kosambi function). Both the linkage map data and the phenotypic data were then imported to QTL Cartographer version 2.0 obtained from http://statgen.ncsu.edu/qtlcart/WQTLCart.htm (Wang et al., 2004). A final number of 141 RILs was used in the QTL analysis. The likelihood, location, additive effect, and percentage of variance explained by each QTL were calculated using model 6 based on the composite interval mapping (CIM) method (Zeng, 1994). QTL cofactors were initially selected by using forward–backward stepwise multiple regression. Mapping was conducted with a walking speed of 0.5cM and a window size of 3cM. For precise determination of significant QTLs, the thresholds of LOD for each linkage group were calculated by a permutation test method (Doerge and Churchill, 1996) with 1000 permutations at the permutation significance level (P <0.05). The support interval of each QTL was constructed using the 2–LOD rule with the confidence intervals being defined by all those values falling within a 2–LOD score of the maximum value (Lynch and Walsh, 1998). To perform EODINDEX1 QTL fine mapping, a new QTL analysis was conducted using the phenotypic data, an additional set of 29 single-feature polymorphisms markers mapping on chromosome 2 (West et al., 2006), and also a marker at the ELF3 locus position that distinguishes the polyglutamine region between both accessions (Tajima et al., 2007). The re-analysis was done as described earlier using the markers corresponding to chromosome 2.
Many RILs are still individually segregating for one limited genomic region. Taking advantage of this, near isogenic lines were generated in the form of heterogeneous inbred families (HIFs) as previously described (Tuinstra, 1997). RIL84 and RIL163 are still heterozygous for markers MSAT2.41 and MSAT 2.7; by screening the F7 seeds with these markers, the region was fixed for each parental allele, generating in this way two lines that differed only at the region of interest (HIF84 and HIF163). Seeds originating from two independent plants for each genotype were used for the experiments. The selected HIFs were screened for the ELF3 marker to confirm that they were not segregating at this particular locus. Conditions for the PCR amplifications were obtained from Loudet et al. (2002).
For quantitative complementation analysis of ELF3 as a candidate gene, Bay and Sha parental lines were crossed with elf3-1 and Col (wild-type control). The F1 and the original genotypes were evaluated for hypocotyl growth under WL and WL+EOD following the experimental protocols previously described.
ELF3 gene sequence was analysed in both parental lines (Bay and Sha). DNA sequence analysis was performed using a dye terminator sequencing system (Big Dye Terminator, Applied Biosystems) on the Genetic Analyzer 3600 (Applied Biosystems). The polymorphisms found were confirmed by those available now at www.arabidopsis.org, and the number of glutamine repeats in Sha was confirmed by Tajima et al. (2007). The polymorphisms found are shown in Supplementary Fig. S1 available at JXB online.
For leaf movement analysis, plants were entrained under SD conditions, until the first pair of leaves emerged, and were then transferred to continuous 40μmol m−2 s−1 white fluorescent light (LL). The position of the first pair of leaves was recorded every hour for 6d (Hicks et al., 1996). Period estimates were calculated with Brass 3.0 software (Biological Rhythms Analysis Software System; available from http://www.amillar.org) (Millar et al., 1995; Plautz et al., 1997).
For qRT-PCR expression analysis, Arabidopsis seedlings were entrained under SD for 2 weeks and then they were released into LL and harvested every 4h. Total RNA was extracted with a Total RNA Extraction Kit-Plant (RBC Real Genomics) following the manufacturer's protocols and it was subjected to a DNase treatment with RQ1 RNase-free DNase (Promega). cDNA derived from this RNA was synthesized using Invitrogen SuperScript III and an oligo(dT) primer. The synthesized cDNAs were amplified with FastStart Universal SYBR Green Master (Roche) using the 7500 Real Time PCR System (Applied Biosystems) cycler. The Protein Phosphatase 2A Subunit A3 (PP2A) gene was used as the normalization control (Czechowski et al., 2005). A list of primer sequences is provided in Supplementary Table S1 at JXB online.
The variation of hypocotyl length in response to shade avoidance signals was studied in seedlings grown in WL or WL+EOD in SD conditions (8h light+16h dark). At the end of the day, half of the seedlings received a pulse of FR (WL+EOD), a light treatment that mimics the shade avoidance responses displayed by plants exposed to vegetation shade in nature (Franklin, 2008). A wide genetic variation for hypocotyl growth was observed in WL and WL+EOD between 30 accessions of Arabidopsis (between 2.21mm and 5.08mm in WL, and between 4.55mm and 8.65mm in WL+EOD). The most common laboratory accessions displayed a very strong inhibition of hypocotyl growth under WL and a strong promotion under WL+EOD. phyB mutant seedlings showed a null response to shade principally due to an impaired inhibition of hypocotyl elongation under WL (Fig. 1A, phyB-1 and phyB-9). Bay and Sha, two accessions broadly used in natural genetic studies because a RIL population was developed, showed different hypocotyl elongation patterns under conditions that simulated open (WL) or shade (WL+EOD) environments (Fig. 1B). Indeed, Sha seedlings displayed 1.7-fold longer hypocotyls than Bay in WL and a similar hypocotyl length to Bay in response to WL+EOD. Therefore, Sha shows a constitutive shade avoidance phenotype similar to that observed in phyB mutants. In the following experiments hypocotyl growth in response to EOD signals among these accessions was explored by QTL analysis.
The frequency distributions of 141 RILs originated from the cross between Bay and Sha (Loudet et al., 2002) in response to WL and WL+EOD are shown in Fig. 2A. A huge variation within the RIL population was detected in WL (between 2.33mm and 10.70mm). Furthermore, phenotypic variation between lines was also observed in WL+EOD (from 4.43mm to 11.53mm). An EOD index was calculated as hypocotyl length in (WL+EOD–WL)/WL+EOD to obtain a relative magnitude of the shade avoidance response for each line. The Sha accession showed a reduced response to EOD in contrast to the Bay accession. The EOD index ranged from 0.047 for RILs less responsive to shade to 0.60 for the most responsive genotypes (Fig. 2B, Supplementary Table S2 at JXB online). Three, two, and four QTLs were mapped in WL, WL+EOD, and the EOD index, respectively. The three QTLs identified in WL had been previously reported (Loudet et al., 2008) when growing seedlings of the same RIL population in continuous WL supplemented with incandescent lamps at 22°C and 26°C. The most significant QTL for the EOD index, which we named EODINDEX1, was mapped on chromosome 2 and co-localized with a QTL that appeared in WL (Fig. 3, Supplementary Table S3 at JXB online). The EODINDEX1 QTL explains 31% of the total phenotypic variation, and Bay alleles have a positive contribution to the promotion of hypocotyl length in response to simulated shade, evaluated through the EOD index trait. As expected, the Bay allelic effects of EODINDEX1 in WL showed a negative additive contribution (Supplementary Table S3 at JXB online).
Additional genetic evidence confirmed the allelic differences detected at the EODINDEX1 QTL. Two independent families of polymorphic lines at the EODINDEX1 QTL were developed. HIF163 and HIF84 carrying Bay alleles at the EODINDEX1 QTL showed a substantially higher EOD index for hypocotyl growth than those carrying Sha alleles (Fig. 4).
The EODINDEX1 QTL contains several genes within the confidence interval defined by MSAT2.36 and MSAT2.7, including the ELF3 gene but excluding PHYB. To reduce the confidence interval of the EODINDEX1 QTL, a fine mapping QTL analysis using 30 additional single-feature polymorphisms including an ELF3 marker was performed. The analysis reduced the confidence interval of the EODINDEX1 QTL to a region of 1399kb that contained potentially two LOD peaks with their confidence intervals (579kb and 820kb, respectively). The highest additive effect was associated with the ELF3 marker, suggesting that it was the most firm candidate gene for the EODINDEX1 QTL (Fig. 5A). In addition, the constitutive shade avoidance phenotypes of Sha and HIF163-Sha resemble the elf-3 phenotype under WL (Figs 1B, 4B, 5C). Then, a quantitative complementation analysis was performed to confirm that the ELF3 gene was responsible for the EODINDEX1 QTL. The phenotype of F1 seedlings coming from a cross between Bay and Sha with elf3-1 and Col (wild-type) exposed to WL and WL+EOD showed that the Sha-ELF3 allele did not fully complement the phenotype of elf3-1. This result suggests that ELF3 is indeed the most likely candidate gene responsible for the EODINDEX1 QTL (Fig. 5B).
In order to find the molecular basis of the phenotypic differences between accessions in response to simulated shade, a study of the genomic sequence of Bay-ELF3 and Sha-ELF3 was performed. Previous evidence demonstrated that ELF3 encodes a protein with a variable number of glutamine repeats, which vary between different accessions (Tajima et al., 2007). However, correlation between the EOD index phenotype and the glutamine number was not found in 22 accessions of Arabidopsis, indicating that this polymorphism is probably not responsible for the differences observed in EOD index between Bay and Sha accessions (Fig. 6A). In addition, eight nucleotide changes were found, only one of them resulting in a non-synonymous amino acid substitution (Supplementary Fig. S1 at JXB online). This substitution corresponds to an alanine at position 362, conserved for all accessions with the exception of Sha, which carries a valine at this position (Fig. 6B). This suggests that this particular amino acid change could be responsible for the molecular basis of the phenotypic variation observed on hypocotyl growth in response to shade.
ELF3 is a component of the circadian clock and has also been implicated in the light control of hypocotyl growth (Hicks et al., 1996; Reed et al., 2000; Liu et al., 2001). The circadian rhythm of leaf movement in the parental accessions Bay and Sha, as well as in polymorphic lines at the EODINDEX1 QTL, was evaluated. Seedlings were entrained under SD conditions for 10d and then transferred to continuous WL for several days (LL, free running conditions). Albeit Bay and Sha accessions showed similar periods for the circadian rhythm of leaf movement, two independent lines of HIFs (84/163) carrying Sha alleles at the EODINDEX1 QTL exhibited an acceleration in the leaf movement rhythm (i.e. a shorter period) compared with those carrying Bay alleles at this position (Fig. 7). Expression analysis by qRT-PCR demonstrated that the differential behaviour of polymorphic lines at the EODINDEX1 QTL was correlated with the expression pattern of six core clock genes and one clock output gene that mediates the clock regulation of hypocotyl elongation. HIF84-Sha showed a small but significant acceleration in the expression rhythm of PRR7 and PRR9 transcripts, already observed during the second day in continuous WL. A similar acceleration was detected for CCA1, LHY, TOC1, ELF3, and TZP transcripts during the third day under free-running conditions when HIF84-Sha is compared with HIF84-Bay (Fig. 8). These results demonstrate that the EODINDEX1 QTL is involved in the regulation of circadian rhythms, as expected if ELF3 is indeed the gene responsible for the EODINDEX1 QTL.
Shade avoidance responses have been shown to confer a significant adaptive advantage in plants growing under natural conditions (Schmitt et al., 2003). At the same time, our understanding of the molecular mechanisms underlying shade avoidance responses has increased dramatically over the last decade (Franklin, 2008). In spite of these advances, the genetic bases underlying phenotypic diversity in responses to shade signals in natural populations is just beginning to be understood. Indeed, whilst some of the molecular mechanisms underlying differential flowering time responses to a low R/FR ratio have been recently uncovered (Adams et al., 2009), the genetic bases responsible for the phenotypic diversity observed in growth responses triggered by shade are completely ignored.
Here attention was focused on the natural genetic variation underlying differential growth responses to simulated shade in the Bay×Sha RIL population. In the present laboratory conditions, the Sha accession displayed an impaired shade avoidance response due to an enhanced promotion of hypocotyl length under WL, and a lack of further growth promotion in response to EOD. EODINDEX1 was identified as the most important QTL in the EOD index response. The generation of polymorphic lines at EODINDEX1 and a quantitative complementation analysis suggest that ELF3, a component of the circadian clock, is the gene candidate for EODINDEX1. An amino acid substitution of a conservative alanine for a ‘rare’ valine amino acid at position 362 of the Sha-ELF3 protein is likely to be the cause of the constitutive shade avoidance phenotype in the Sha accession. In addition, physiological and RNA expression analysis comparing the rhythm of leaf movement between genotypes demonstrate that Sha alleles at EODINDEX1 shortened the period of circadian rhythms, these effects being dependent of the genetic background. Null elf3 plants are arrhythmic (Hicks et al., 1996), but plants overexpressing ELF3 display a long circadian period phenotype (Covington et al., 2001), suggesting that the ELF3 allele present in Sha is hypomorphic.
ELF3 is a highly conserved plant-specific nuclear protein, which regulates flowering time, photomorphogenic development, and circadian rhythms in Arabidopsis (Hicks et al., 1996; Zagotta et al., 1996). The precise biochemical and molecular mechanisms through which ELF3 acts to regulate the above physiological processes are largely unknown. Several years ago a role for ELF3 in mediating or modulating the activity of phyB, the main SAS photoreceptor, was proposed based on the interaction of these two proteins in the yeast two-hybrid system (Liu et al., 2001). At the same time, the physiological characterization of single and double mutants of these genes revealed both phyB-dependent and -independent effects of ELF3 on hypocotyl growth (Reed et al., 2000; Liu et al., 2001). In relation to its circadian function, ELF3 was originally thought to allow circadian rhythm progression modulating the response to light cues that set the circadian clock to the correct time of day (McWatters et al., 2000; Covington et al., 2001). More recently, however, ELF3 has been proposed to act as a component of the core oscillator to allow clock progression from day to night in a light-independent fashion (Thines and Harmon, 2010). Thus, both light-dependent and light-independent effects of ELF3 on different physiological processes have been reported, suggesting that ELF3 may have multiple molecular functions.
The present data support the idea that ELF3 is a component of the light signalling pathway regulating SAS, in agreement with previous evidence (Reed et al., 2000), and this function appears to be independent of ELF3’s role as a component of the circadian clock. EODINDEX1 mediated the acceleration of circadian rhythms, at least in part, by regulating the expression of some core clock genes. This observation is in agreement with the result of the quantitative complementation approach suggesting that the ELF3 locus is responsible for the impaired SAS phenotypes of Sha. Interestingly, the effect of the Sha EODINDEX1 allele on clock function was evident in the two different HIF backgrounds evaluated, but was not present when the parental line backgrounds were compared. In contrast, Sha EODINDEX1 effects on hypocotyl growth responses to shade were impaired both in the parental line background and in the HIF backgrounds. Therefore, although the long hypocotyl phenotype of elf3 has been proposed to result from its defects in clock function, the lack of correlation between circadian and SAS alterations caused by EODINDEX1 in different genetic backgrounds supports the idea of a direct role for ELF3 in phyB signalling, in addition to its well-established role as a circadian clock component (Thines and Harmon, 2010). The recent observation that ELF3 may act as an adaptor protein facilitating the ability of the E3-ubiquitin ligase COP1 of destabilizing GIGANTEA (Yu et al., 2008) suggests that ELF3 may affect clock function and light signalling independently regulating the stability of different proteins.
In addition to identifying ELF3 as the first candidate gene to mediate phenotypic variation in SAS, the present work also expanded previous knowledge on the phenotypic variation in different shade avoidance responses among A. thaliana accessions (Botto and Smith, 2002). Specifically in this work, a wide range of variation for the control of hypocotyl growth under WL and WL+EOD in 30 accessions of Arabidopsis was shown. The most common laboratory accessions displayed reduced hypocotyl elongation in WL and a clear response to EOD that mimics a natural shade avoidance response. However, a few accessions such as CT-1 from Catania (Italy) and No-0 from Nossen (Germany) already showed large hypocotyl phenotypes similar to those observed in phyB mutants under WL. Although this study explores the molecular bases of the natural variation in SAS, the genetic architecture of this central feature of plant development is far from being understood. The recent availability of many segregating populations (Simon et al., 2008; Balasubramani et al., 2009), coupled with the study of other shade avoidance responses, will surely provide very valuable new alleles and genes involved in the SAS. For example, QTL mapping using RILs originated from the cross of Ct-1 or No-0 with other contrasting parental genotypes will be a fruitful avenue to find new loci involved in the SAS. Furthermore, the comparison between different RIL populations sharing a parent, such as the Ct-1 accession, should contribute to improve the understanding of the function of loci involved in the SAS in different genetic environments.
Finally, although the natural adaptive advantage of shade responses using mutagenic approaches has been demonstrated (Schmitt et al., 1995; Ballaré and Scopel, 1997), the potential adaptive advantage of natural genetic variants in plastic responses to shade remain to be explored (Mitchell-Olds and Schmitt, 2006; Shindo et al., 2007). To determine whether natural variation in SAS indeed resulted from local adaptations among Arabidopsis accessions, a substantial effort in ecological analysis is required. The adaptive contribution of EODINDEX1 to different density environments will have to be assessed through the generation and evaluation of near isogenic lines in field conditions under different levels of shading, to estimate the relative fitness of Bay and Sha alleles. The shade avoidance hypothesis would predict that the fitness of Sha-EODINDEX1 plants relative to Bay-EODINDEX1 plants would decline at low density. In summary, the difficult ecological question of whether the Sha-allele at EODINDEX1 is a deleterious variant that has become fortuitously fixed in the population, or if this allele indeed provides a fitness advantage in its local environment will soon be answerable.
Supplementary data are available at JXB online.
Table S1: List of primer sequences used for qRT-PCR expression analysis.
Table S2: Summary statistics for hypocotyl length under WL, WL+EOD and EOD index.
Table S3: QTL mapping for hypocotyl length under WL, WL+EOD and EOD index in Bay × Sha RIL population.
Figure S1: Polymorphisms found in ELF3 sequence of Bay and Sha accessions.
We thank Romina Sellaro for her invaluable technical support. This research was supported by Préstamo BID PICT 32137 and UBACyT grant G041 to JFB, and Préstamo BID PICT 1026 to MJY. The work was supported by doctoral fellowships from CONICET (MPC and SES) and ANPCyT (LK).