- Research article
- Open Access
Genetic architecture of spring and autumn phenology in Salix
BMC Plant Biology volume 14, Article number: 31 (2014)
In woody plants from temperate regions, adaptation to the local climate results in annual cycles of growth and dormancy, and optimal regulation of these cycles are critical for growth, long-term survival, and competitive success. In this study we have investigated the genetic background to growth phenology in a Salix pedigree by assessing genetic and phenotypic variation in growth cessation, leaf senescence and bud burst in different years and environments. A previously constructed linkage map using the same pedigree and anchored to the annotated genome of P. trichocarpa was improved in target regions and used for QTL analysis of the traits. The major aims in this study were to map QTLs for phenology traits in Salix, and to identify candidate genes in QTL hot spots through comparative mapping with the closely related Populus trichocarpa.
All traits varied significantly among genotypes and the broad-sense heritabilities ranged between 0.5 and 0.9, with the highest for leaf senescence. In total across experiment and years, 80 QTLs were detected. For individual traits, the QTLs explained together from 21.5 to 56.5% of the variation. Generally each individual QTL explained a low amount of the variation but three QTLs explained above 15% of the variation with one QTL for leaf senescence explaining 34% of the variation. The majority of the QTLs were recurrently identified across traits, years and environments. Two hotspots were identified on linkage group (LG) II and X where narrow QTLs for all traits co-localized.
This study provides the most detailed analysis of QTL detection for phenology in Salix conducted so far. Several hotspot regions were found where QTLs for different traits and QTLs for the same trait but identified during different years co-localised. Many QTLs co-localised with QTLs found in poplar for similar traits that could indicate common pathways for these traits in Salicaceae. This study is an important first step in identifying QTLs and candidate genes for phenology traits in Salix.
Adaptation of the annual cycles of growth and dormancy to the local climate is critical for survival and competitive success of woody plants. Such local adaptation has been described in several species, and it can be seen as clines in phenology traits, including the timing of leaf emergence, leaf senescence, and growth cessation . This pattern reflects the trade-off between frost tolerance and enhanced growth . Cessation and initiation of growth determines the period of active stem elongation, mark the shift between frost resistant and vulnerable phases, and their timing is critical for overall biomass production, fitness and the long-term survival of species [3, 4]. After height growth cessation and until leaves senesce, deciduous shrubs and trees, including willows and poplars, continue to be photosynthetically active and may accumulate considerable biomass [5–7]. In this phase, photosynthesis is critical to cold acclimation and survival in winter . The timing of leaf senescence in autumn has a strong impact on nutrient retranslocation, reserve storage and the next early-season growth [5, 9]. In addition, spring and autumn leaf phenology are evolutionary important traits for herbivore and pathogen resistance [10–13].
Decreasing photoperiod (day-length) is the main environmental cue inducing growth cessation and bud set in many perennial plants , including poplar . This response to photoperiod is under strong genetic control [16–18] and is maintained when trees are moved between latitudes . Temperature, alone or in combination with photoperiod, also induces growth cessation in some tree species [20–24]. Senescence and shedding of leaves are also influenced by photoperiod, often in interaction with temperature . In poplars, leaf senescence is induced by shortening day-lengths but needs to be preceded by bud set. However, it is not clear to what extent senescence and bud set are under independent genetic control [26, 27].
After growth cessation, dormancy is initiated, which is a prerequisite for the development of cold acclimation and freezing tolerance . Dormancy release requires exposure to chilling temperatures . Light does not seem to play a role in this process, but once dormancy is broken, bud burst and growth resumption are regulated by temperature, light and photoperiod. Growth resumption and bud burst depend on accumulation of temperature units over a specific threshold (thermal time) . Temperature is the most important factor regulating bud burst in temperate woody plants , but photoperiod also plays a role in some populations and species [31, 32].
Willows (Salix) belong together with poplars (Populus) to the plant family Salicaceae. Based on the fossil record the divergence of the two genera was dated to approximately 45 mya [33, 34]. Willows and poplars share many characteristics such as dioecy, rapid growth and seed development, and ease with which they can be vegetatively propagated. They typically have a haploid chromosome number of 19 and similar genomes sizes of approximately 500 Mbp. Also there is strong syntheny and colinearity between willow and poplar genomes . The Salix genus shows a remarkable phenotypic diversity ranging from small shrubs to large trees. Salix spp. have a global distribution in temperate and arctic regions and are adapted to a wide range of habitats . Relatively high levels of genetic diversity  and the broad phenotypic diversity make them an excellent model system for studying evolutionary processes such as adaptation. Moreover, willows have generally rapid growth and high biomass yields and these characteristics together with ease of vegetative propagation make them economically attractive as bioenergy crops. Willows have been increasingly used in the last decades for biomass production worldwide and Salix viminalis L. and S. schwerinii E. Wolf and their hybrids are some of the most commonly used willows in the breeding programs in Europe. These two species are dioecious and outcrossing and morphologically very similar. Both are multi-stemmed shrubs with long and slender leaves and are commonly found along rivers and in other wet areas.
Growth cessation in Salix species is marked by the abscission of the shoot apex  and is controlled by photoperiod [38, 39]. Large clonal variation in the time of leaf abscission has been observed in Salix, and delayed leaf abscission was shown to impair leaf nitrogen retranslocation and to increase nitrogen losses . Extensive clonal and species variation in timing of bud burst has been observed in willow, mainly determined by differences in thermal time requirement [6, 40, 41]. Moderate to high heritabilities have been reported for timing of bud burst and growth cessation in different S. viminalis families [40, 42, 43].
Phenology traits have a quantitative genetic background and thus QTL mapping is a powerful method to identify genomic regions controlling phenology traits. With a reference genome one can obtain information on the genomic content of the QTL regions. With the advent of high-throughput genotyping technologies and annotated reference genomes, genetic markers in evenly spaced genes throughout the genome can be developed and genotyped for the purpose of constructing dense genome-wide linkage maps.
Here we studied the phenology of growth, including timing of bud burst, timing of cessation of elongation growth and leaf abscission in willows both in controlled and field conditions during multiple years. Our two aims were to map QTLs associated with phenology in Salix, and to identify candidate genes in QTL hot spots through comparative mapping with Populus trichocarpa. This was achieved by using a dense linkage map anchored to the annotated genome of P. trichocarpa and by constructing denser maps in the QTL hot spots.
Phenotypic variation in phenology
In the S1 pedigree, planted in an experimental field in Pustnäs, south of Uppsala (59°48′ N, 17°39′E, 25 m), the mean date for bud burst was 20th of April (day of the year (DOY) 112, Figure 1), the date of apex abscission (growth cessation) was 24th of September (DOY 268, Figure 1), and about 25% yellow leaves were left on the plants at the end of October (leaf senescence index LSI = 1.5; DOY 304, Figure 1). Year-to-year variation was significant for all traits (Table 1, Figure 1), and genotype ranking significantly changed across years (Table 1, significant genotype × year interaction) indicating different responses to seasonal variation between genotypes. In the indoor experiment, elongation growth ceased on average 2 weeks after progressive reduction of the photoperiod. After nine weeks of artificial winter in the indoor experiment, bud burst occurred after 4 weeks of increased temperature and day length, corresponding to 336 day degrees > 0°C (Figure 1g, h).
In the progeny, all traits varied significantly among genotypes (Table 1, Figure 1). The variation of genotypic means for bud burst was smaller in the field than indoors (11 vs. 17 days, Figure 1b, h). A large part of the phenotypic variation for all the traits and in all experiments was due to genetic factors as shown by relatively high broad-sense heritabilities (0.5–0.9) (Table 2). The highest broad-sense heritability was found for leaf senescence. There were strong block effects (Table 1), but all genotypes were similarly affected in the field (Table 1, no genotype × block interaction) and indoors (among-blocks correlation r=0.5–0.7, p<0.001). Therefore, QTL analyses were performed on unadjusted clonal means, which should be a good estimate of the average behaviour of a genotype.
Bud burst and growth cessation showed a significant correlation only in 2008. Positive but weak correlations between the field and indoor were found both for bud burst and growth cessation. In the field, bud burst, growth cessation, and leaf senescence were each positively correlated between years (Table 3, plots in Additional file 1: Figure S1).
In total across experiments and years, 80 QTLs were detected (Table 4, Figures 2 and 3). QTLs were named by the trait (BB = bud burst, GC = growth cessation, LS = leaf senescence), environment (I = indoors or F = field), year of assessment, linkage group where it was located, and if multiple QTLs were mapped to the same group, by an ordering number. On average about 40% (21.5–56.5% depending on the trait) of the variation among genotypic means was explained by a general model including all the QTLs that remained significant after backward selection. Individual QTLs generally explained a small proportion of the variance of the genotypic means. Eleven QTLs contributed to 7% or more of the trait variation, 3 QTLs explained more than 15%, and 1 QTL explained up to 34% of the variation (Table 4).
QTLs for bud burst
For bud burst, 9 QTLs were identified indoors and 26 in the field during two years of assessments. Each of the QTLs included in the final model explained 1.0 to 14.1% of the variation in genotypic means (Table 4). The final model including all QTLs explained from 27.6% to 53.2% of the observed variance depending on environment and year of assessment (Table 4). No significant QTL by QTL interaction was detected for bud burst.
QTLs for growth cessation and leaf senescence
For growth cessation, 3 QTLs were identified indoors and 18 in the field during two years of assessments. The contribution of each QTL to the phenotypic variance was low to moderate (1.2% - 16.7%) but all QTLs together explained from 21.5 to 49.1% of the trait variation (Table 4). For leaf senescence, 24 QTLs were identified across the two years of assessment. The final model (10 QTLs) explained 56.5% of the trait variation in 2009 and 43.8% in 2010. The majority of QTLs made a low contribution to the total variance (1.0 – 6.0%) with the exception of LSF09.II and LSF10.II, which explained 34.2% and 19.3% respectively (Table 4).
Comparison of QTL positions among years, environments and traits
QTLs were considered to co-localize when their peak positions were less than 10 cM apart. Twenty-eight QTLs were identified only once while at 20 other positions at least 2 QTLs were mapped. Among the 20 regions where co-localizing QTLs were found, 11 included 2 QTLs for the same trait identified in 2 years or both environments, 4 included 2 QTLs for different traits, and 5 included 3 to 6 QTLs affecting different phenology traits (Table 4, Figures 2 and 3).
For bud burst in the field, in five cases QTLs appeared more than once at the same map position comparing the two years (Figures 2 and 3). In contrast, for growth cessation in the field, most of the QTLs mapped to the same position both years (Figures 2 and 3). Six genomic regions were detected for leaf senescence where QTL appeared both years (Figures 2 and 3).
In only one case were QTLs involved in bud burst found at the same map location when comparing indoors and the field conditions (Figure 2). One QTL for growth cessation was only identified indoors while another two co-localized with QTLs identified in the field (Figures 2).
Improved mapping of two QTL hot spots for phenology on LG II and LG X
The QTLs on LG II and X for bud burst, growth cessation and leaf senescence were confirmed with the new denser linkage map. A comparison of the QTL mapping with the linkage map from Berlin et al.  and the new maps is illustrated in Figures 4 and 5. The peak position of the QTLs was changed for all traits on LG II while only to some extent for the growth cessation traits on LG X. The 1.5 LOD intervals of the QTLs were considerable shortened with the new linkage maps (Figures 4 and 5). For the traits BBI and GCF09 the new analysis identified two QTLs on LG II for each trait instead of one as in the original analysis (Figure 4).
Positional information of QTLs and identification of candidate genes
We predicted the genomic interval for each QTL and summed all intervals for the traits, BBF09, BBF10, BBI, GCF08, GCF09, GCI, LSF09 and LSF10. This was possible by anchoring the SNP markers to the P. trichocarpa genome. The total genomic intervals varied between 15.1 and 69.0 Mbp for the traits (interval for each QTL see Additional file 2: Table S1) and the total number of gene models ranged from 565 to 6,604. The largest number of gene models was found for BB10, however, 39% is due to one QTL with a large 1.5 LOD interval on LG I. Since some QTLs covered the same genomic interval, some gene models appeared more than once, we therefore estimated the number of unique gene models for BB, GC and LSF to 9,633, 4,355 and 4,815 respectively (Table 5). The total number of gene models in the intervals as well as the putative candidate genes for growth cessation are presented in Table 5. Several putative candidate genes were identified among these gene models (Figures 2 and 3; Additional file 3: Table S2). Candidate genes in the QTL intervals include photoreceptors as well as several circadian clock genes and downstream components. Among photoreceptors both cryptochrome and phytochrome genes were identified (Figures 2 and 3; Additional file 3: Table S2). The phytochrome gene PHYB2 was located in a narrow range between two QTLs on LG X and should not be ruled out as a potential gene influencing growth control in willows. Gene models within QTLs also included several core circadian clock genes such as LATE HYPOCOTYL (LHY), PSEUDO-RESPONSE REGULATOR 7 (PRR7), LUX ARRHYTHMO (LUX), EARLY FLOWERING 3 (ELF3), ZEITLUPE (ZTL) and GIGANTEA (GI) (Figures 2 and 3; Additional file 3: Table S2) . Interestingly the FT2 gene is located in the proximity of the fine mapped region on LG X (Figure 5). Noteworthy is that PtFT2  and PttLHY  are involved in the control of growth cycle in Populus.
The present study explores the genetic architecture of growth phenology in a pedigree between S. viminalis and S. schwerinii. Using a linkage map based on several hundreds of SNP markers, we identified QTLs for bud burst, growth cessation and leaf senescence in different years and environments. In the field, the QTLs explained together more than 40% of the variation in each trait. Several regions were identified where many QTLs co-localized for different traits or for the same trait across years. Since the SNPs were developed using the P. trichocarpa genome as a template, we obtained positional information for the Salix QTLs, projected them on the poplar genome, and identified the corresponding genomic intervals. The results suggest that some QTLs might be homologous to Populus QTLs. Moreover, in the projected QTL intervals we could identify putative candidate genes for the traits.
All phenotypic traits varied among the progeny. A large part of the variation was imputable to genetic factors, which confirms moderate to high broad sense heritabilities for bud burst and growth cessation in Salix species [42, 43]. Although heritability estimates are known to be environment and population specific, it is now well established that bud phenology is under strong genetic control in Salicaceae [16, 18, 47–51].
The behaviour of the progenies changed between controlled and field conditions. Growth cessation in short days indoors and apex abscission in the field displayed particularly weak correlations. One possible explanation to this discrepancy could be that the variation observed indoors only reflects a photoperiodic response, while the variation observed in the field could reflect effects of other environmental factors. In fact, short days alone can induce apical growth cessation both in seedlings and rooted cuttings of several Salix species, while in field conditions, apical growth cessation does not seem to be regulated by photoperiod, but other factors seem to be involved [38, 52]. In S. viminalis, the coincidence of growth rate decline and tip senescence with development of low leaf water potential in summer suggests an effect by water stress .
The timing of phenological events changed across years, as shown by Weih  for other Salix clones. In addition, the ranking of the clones changed from one year to another, which indicates a plastic behaviour in the family. This might represent a differential response among the progenies to the seasonal differences in rainfall and temperature observed between years (Figure 6). Plasticity is crucial for a species to respond to the demands of a changing environment. An increasing number of studies indicate that temperature and stress factors may variably interact with photoperiod in controlling the timing of phenological events in woody species, including willow and poplar [20, 53, 54]. In S. viminalis, the effects of photoperiod seem to be superimposed on those of water stress in controlling apical growth cessation . A drier and warmer summer, i.e. higher temperature sum, higher maximum temperatures and a greater number of dry days (Figure 6), might indeed be partly responsible for the markedly earlier growth cessation observed in 2009 in this study. Combined effects of photoperiod and temperature could be also responsible for the inter-annual variability observed for leaf senescence in our willow pedigree. In Populus, yellowing of the leaves is initiated by a photoperiodic stimulus, but the progression of senescence is accelerated under low temperature . Moreover in Populus, senescence seems to be faster in trees that have a late onset of senescence, independently on the effect of temperature. There is no data about regulation of seasonal leaf senescence in Salix. However, if we hypothesize a similar regulation as in poplar, the smaller total phenotypic variation (the full range of phenotypic variation in leaf abscission between individual plants was 40 - 100% in 2009, while 10 - 100% in 2010; not shown) and the more advanced stage of senescence observed in 2009 compared to 2010 could both be explained as an effect of lower temperatures (means, minima and maxima) in October 2009 (Figure 6a, b, c).
QTLs for phenology traits
The phenology traits were typically quantitative. Several QTLs were found for each trait of which the majority explained less than 10% of the total variation. This was not surprising and seems to be a common feature of phenological traits among many species in the Salicaceae plant family [18, 50, 51]. Indeed, the regulation of phenology is quite complex and involves different pathways [14, 55–57]. When estimating the total phenotypic variation of all QTLs for each trait we found particularly high values for leaf senescence, e.g. LSF10 explained as much as 56.5% of the total variation. The total phenotypic variation explained by the QTLs for the other traits was less strong but still substantial.
Colocalization of QTLs
At the level of resolution attained, there is considerably coincidence in map positions among a number of the QTLs. The most striking instances are found on LG II and X, where QTLs with narrow confidence intervals for BB, GC and LS consistently mapped across years and environments are located at similar positions. The QTLs on LG II also co-locate with bud burst QTLs found by Tsarouhas et al.  where they used a subset of the same willow pedigree and a different linkage map with some common markers (Figures 2 and 3). Fine mapping with additional markers on LG II and LG X further supported that QTLs for all three traits might represent the same locus. Similarly, on LG V, QTLs for growth cessation and leaf senescence co-localized, on LG VI, QTLs for leaf senescence and bud burst co-localized and on LG XIV QTLs for bud burst and growth cessation were found at a similar genomic position. This suggests that pleiotropic effects of individual QTLs on several traits could result from shared components of the pathways controlling the different traits. This is perhaps not unexpected given that both bud burst and growth cessation are both controlled jointly by photoperiod and temperature. The time of leaf senescence has been shown to be an important trait for seasonal acclimation as well as for biomass production in willows . It generally occurs well after growth cessation, but data from Populus tremula suggest that also the onset of leaf senescence is under photoperiodic control although it might respond to a different photoperiod than bud set . The frequent co-localization of QTLs for growth cessation and leaf senescence suggests that genes controlling growth cessation may to a large extent also affect the timing of leaf senescence. Furthermore, the common QTLs on some of the linkage groups suggest that these locations indeed contain genes with a central role in controlling seasonal growth in willows.
Comparison of QTL regions between poplar and willows
When comparing the willow QTLs with QTLs found in poplar many QTLs were found in similar genomic positions, which suggests common mechanisms controlling these traits in willows and poplar. Rohde et al.  mapped six narrow QTLs for bud set in poplar to LG III, V, VI and two on VIII and XIII, and we find QTLs for growth cessation or leaf senescence at five of these positions (Figures 2 and 3). Fabbrini et al.  identified several QTLs for bud set in P. nigra of which those on LG I and IV overlap with QTLs for growth cessation and leaf senescence in willows. Frewen et al.  mapped three QTLs in poplar for bud set to LG III, VI and X of which those on III and VI was also identified in willows. There are also some major differences as for example the QTLs on LG II and IX in willows are missing in poplars. However the presence or absence of a QTL for a trait depends on several factors such as accuracy of the phenotyping, environmental factors and perhaps most importantly on whether or not the trait or gene is variable in the mapping population under study. In this study we have used one mapping population from a back cross between two species and how general these results are should be validated in another genetic background.
Identification of positional candidate genes
Comparative mapping with poplars represents a first step to further dissect the genetic basis of phenological variation in willows and ultimately identify genes or alleles responsible for the trait variation that we observe. Our understanding of the physiology and biochemistry of the traits of interest very much define the success of this approach as it relies on previous identification of genes potentially involved in the control of the traits as well as on genome conservation between poplars and willows. We have previously shown that overall gene order is conserved between willows and poplars except for few large-scale chromosomal rearrangements , which justify anchoring the Salix QTLs to the physical map of poplar. We started by estimating the number of genes within the QTLs for each trait and found as many as 9,633 genes for bud burst and about half the number of genes for growth cessation and leaf senescence. Functionally characterized genes for phenology are still scarce in poplar apart from two examples. The poplar gene PtFT2 located on LG X that is known to be involved in growth cessation , was located close to the fine mapped region. LHY-genes have recently been shown to be involved in both growth cessation and budburst  and LHY1 was located close to the QTL hot spot region on LG II and LHY2 was located in a cluster of QTLs for bud burst and growth cessation on LG XIV. The fine mapped region on LG II contains a SVP homolog. In perennial species SVP genes have been shown to be involved in the growth cycle typically with high expression levels during bud dormancy [58–61]. The above-mentioned genes are strong candidate genes for growth cessation and bud burst that warrant further investigation in willows.
We identified substantial variation in all traits in the pedigree and all traits were associated with many QTLs that each explained less than 10% of the variation, a typical pattern of quantitative characters. In total, we identified 80 QTLs, of which some were clustered in hotspots where QTLs for the different traits co-localised. Two such hotspots on LG II and X were further investigated by the construction of denser linkage maps in these regions, an effort that greatly reduced the QTL intervals (and number of gene models). Some QTLs appear to co-localize with those found in poplars, which could indicate common pathways for these traits in Salicaceae. This study is an important first step in identifying QTLs and candidate genes for phenology traits in Salix but further work is needed e.g. to confirm the QTLs in other genetic backgrounds, further fine mapping and functional studies, to verify candidate genes.
Plant material, experimental design and phenotyping
The S1 pedigree consists of 463 F1 progeny from the cross between the diploid hybrid male ‘Björn’ (Salix viminalis L. × S. schwerinii E. Wolf) and the diploid female S. viminalis ‘78183’ originating from southern Sweden. The S. schwerinii parent (79069) of Björn originates from Siberia while the S. viminalis parent is an interspecific cross between the male clone 78101 from Western Sweden and the female clone 78195 from southern Sweden. The parental clones of S1 were selected based on variation in phenology traits . The pedigree is conserved in an orchard near to Uppsala (59°49′ N 17°40′ E, central Sweden) where all the plant material used in this study was collected.
Experiment 1 - Growth cessation and bud burst indoors
In spring 2008, growth cessation was assessed in the S1 pedigree in a phytotron experiment under controlled day length and temperature conditions. The two parental genotypes and 294 randomly drawn genotypes were propagated by means of hardwood cuttings and planted in 1.1 litre pots filled with Weibulls ‘Kron Mull’ (organic matter 95%; pH 5·5–6·5; 180 g m-3 N, 110 g m-3 P, 195 g m-3 K, 260 g m-3 Mg, 100 g m-3 S, 2000 g m-3 Ca). All genotypes were assigned to each one of three walk-in growth chambers in a complete randomized block design, where each chamber represented a block with one replicate. Plants were grown for four weeks under 20°C constant temperature, 70% relative humidity and 20 h photoperiod (300 μmol PAR m-2 s-1). The day length was then reduced to 16 h for one week, to 14 h for another week, and then one hour per week down to 10 h. After the first two weeks of growth, the plants were pruned and only the main shoot was preserved for the experiment. Plant height, i.e. the length of the stem from the emerging point on the wooden cutting to the tip of the apex, was measured once a week during week three and four, and three times per week thereafter. A sigmoid curve on the form y=b/(1+a × e(-k × x)), where y is plant height and x is the number of days since the first measurement, was fitted to the data of each plant in order to estimate the end of elongation growth, i.e. the day on which the estimated plant height reached 95% of the final value. The date of growth cessation was expressed as the number of days since the start of the reduction of daylength.
After the simulated autumn, the plants were subjected to an artificial winter at 8°C constant temperature and 9 h photoperiod for 9 weeks. The plants were then cut back, leaving 5 cm stem. Bud burst was forced by keeping a constant temperature of 12°C for six weeks and 11 h photoperiod. During this time the plants were checked for bud burst daily or every second day, depending on the speed of the process. Bud burst was defined as stage 3 according to the phenological scale previously used by Weih . Individual plants were recorded as flushing on the day when at least one bud reached stage 3. The date of bud burst was expressed for each plant as the number of days since the beginning of forcing.
Experiment 2 - Bud burst, growth cessation and leaf senescence in the field
In spring 2008, plants of 463 genotypes were planted in an experimental field in Pustnäs, south of Uppsala (59°48′ N, 17°39′E, 25 m), at a spacing of 130 × 50 cm (i.e., about 20,000 plants ha– 1), according to a randomized complete block design comprising six blocks with one plant per genotype in each. Two border rows were planted around the experiment to reduce marginal effects. The plants were obtained by 5 cm hardwood cuttings rooted in 0.5 L peat pots with Weibulls ‘Kron Mull’ as growing medium. The plants were grown for five weeks in a greenhouse and then transferred outside for hardening. Before planting, the site was appropriately prepared , including ploughing, harrowing and repeated application of a systemic herbicide (Glyphomax, Dow AgroSciences, Indianapolis, IN). The plantation was irrigated in summer 2008 and weed controlled during the whole experimental period. The plants were cut back in winter 2009 and fertilized in spring 2009, 2010 with N P K (21-4-7) corresponding to 80 kg N/ha and year.
Growth cessation as defined by shoot apex abscission was scored using the highest shoot of each plant in 2008 and 2009 from the end of August and once to twice per week depending on the rate of progression. Leaf senescence and abscission was visually estimated on October 31 2009 and 2010 according to the following leaf senescence index (LSI): 0 = no leaves left on the plant (100% abscission); 0.5 = less than 10% brownish leaves (~ 95% abscission); 1= 10 to 20% brownish leaves (~ 85% abscission); 1.5 = 20 to 30% brownish or yellow leaves (~75% abscission); 2 = 30 to 40% yellow leaves (~65% abscission); 2.5 = 40 to 50 yellow and green leaves (~55% abscission); 3 = 50 to 65 green leaves (~40% abscission); 3.5 = 65 to 80% green leaves (~30% abscission); 4=more than 80% green leaves (~10% abscission). Bud burst, defined as in Experiment 1, was assessed twice a week during April and May 2009 and 2010. Date of bud burst and date of apex abscission were expressed as day of the year (DOY), i.e. number of days since January 1. Plants were cut back in January 2009. Therefore spring phenology was assessed on stumps in 2009 and on one-year old shoots in 2010; and autumn phenology was assessed on one-year shoots in 2008 and 2009, and on two-year old shoots in 2010. In order to characterize the weather conditions at the plantation site, temperature sum (cumulated mean temperatures CMT,°C), number of chill days with mean temperature below 5°C (CD<5°C), and the number of days without rain (cumulated dry days, CMDD) were calculated in relevant periods of time in all years of study from the records of a nearby meteorological station (Figure 6).
Statistical analyses and QTL mapping
The complete set of data from each experiment, which included all genotypes that had records for at least three ramets, was analysed with a mixed model ANOVA and variance components analysis to determine the effects of genotype (set as random factor) and block (set as fixed factor), and to estimate the broad-sense heritability, i.e. the ratio between genetic variance and total phenotypic variance. The genotype × block interactions in the field were tested using grouped data from two years. Subsequent analyses were performed on unadjusted mean values among blocks in each experiment. Correlations among genotype means were calculated and plotted (Additional file 1: Figure S1) for all combinations of traits in each experiment and across experiments. The year effect and its interaction with genotype were tested in a separate ANOVA for all the traits assessed in the field. QTL analyses were performed with MapQTL ® 6.0  and the linkage map previously developed for the Salix pedigree S1 was used . First, interval mapping  was applied using 1.0 cM steps across the genome to determine putative QTLs involved in the variation of each trait. Following Churchill and Doerge , the logarithm of the odds (LOD) threshold for QTL significance was empirically estimated from 5,000 permutations of phenotypic data. Two theoretical critical thresholds were considered for detection of a putative QTL: the first corresponding to genome wide error rate of 5% was used to define significant QTLs, and the second corresponding to a type I error of 5% at the chromosome level was used to define suggestive QTLs (‘suggestive linkage’, [66, 67]. Multiple QTL mapping (MQM)  was then performed on the same data: the nearest SNP marker to each putative QTL peak was used as a cofactor to control the genetic background while testing at a position in the genome. Only markers close to QTLs significant at the genome wide level were used as cofactors in MQM. When a cofactor was also a flanking marker of the tested region, it was automatically excluded from the model. The number of cofactors used varied between 1 and 5. A 1.5-unit drop off on either side of the local LOD score peak was used to determine heuristic support intervals for significant QTLs.
A backward selection procedure, at a test level of 0.05, was performed to determine whether a particular QTL could be dropped from the model resulted from the MQM analysis. For each trait, the model involved the genotype at the closest marker to the corresponding putative QTL. Pairwise epistatic interactions between all putative QTLs were tested, via the corresponding marker × marker interaction treated as random effects, by comparing the deviance of a model including all main effect QTLs and the specific epistatic interaction, with the deviance of a model including only the main effect QTLs and no epistatic interaction . The deviance test (p ≤ 0.001) was applied for testing interactions. The contribution of each significant QTL and epistatic interaction was then estimated by maximum-likelihood variance components analysis. Statistical analyses were performed using R  and SPSS v.19.
Improved mapping of two QTL hot spots for phenology on LG II and X
Genomic DNA from 463 genotypes from the S1 pedigree was extracted as in Berlin et al. (2010). Two chromosomal regions on LG II (75–90 cM) and LG X (84–93 cM ) containing QTLs for all traits except bud burst in the field 2010 were selected to increase the marker density and new primer pairs were designed in these regions (Additional file 4: Methods and Additional file 5: Table S3). Gene segments were amplified by PCR and sequenced in each parent as described in Berlin et al. . A total of 96 SNPs were selected for genotyping, 59 from LG II region and 37 from LG X region. Linkage maps were constructed as in Berlin et al. , for details see Additional file 4: Methods.
New QTL analyses for group II and X with the denser linkage maps were conducted. In the linkage map used for the QTL analysis all AFLP markers were removed since only 96 out of 463 individuals were genotyped with AFLPs and thus did not add much information to the analysis. When several markers were located at the same position only the most informative was kept. MQM mapping using MapQTL ® 6.0  was conducted for comparison of QTLs between new and earlier linkage maps. The most informative marker close to the peak position of the QTL was used as cofactor in the MQM analysis.
Comparative mapping and in-silico selection of candidate genes
Physical coordinates of the QTLs were obtained from anchored markers in the P. trichocarpa genome assembly version 3 (http://www.phytozome.net/poplar). For each QTL interval, positions of sequences containing the SNP markers flanking 1.5-LOD on both sides of the LOD score peaks were determined by BLASTN searches. Once the regions were determined, gene models (predicted by the Gnomon gene prediction tool) were downloaded using BioMart that were subsequently annotated using the P. trichocarpa version 3 annotation information. Annotated gene models positioned in QTLs were searched for putative candidate genes. Since growth cessation and bud set are mainly controlled by photoperiod, genes in the photoperiod pathway and the circadian clock were considered as candidate genes for these traits (Additional file 3: Table S2). Genes controlling budburst and leaf senescence are less well known and candidate genes for those traits with functional characterization are scarce or absent. We therefore did not attempt to identify any such candidate genes.
Morgenstern EK: Geographic Variation in Forest Trees: Genetic Basis and Application of Knowledge in Silviculture. 1996, Vancouver, Canada: UBC Press
Levitt J: Growth and survival of plants at extreme temperature – A unified concept. Dormancy and Survival, Symposia of the Society for Experimental Biology, XXIII. Edited by: Woolhouse HW. 1969, Cambridge: Cambridge University Press, 395-448.
Chuine I: Why does phenology drive species distribution?. Phil Trans R Soc B. 2010, 365: 3149-3160. 10.1098/rstb.2010.0142.
Chuine I, Beaubien EG: Phenology is a major determinant of tree species range. Ecol Lett. 2001, 4: 500-510. 10.1046/j.1461-0248.2001.00261.x.
Nelson ND, Isebrands JG: Late-season photosynthesis and photosynthate distribution in an intensively-cultured Populus nigra x laurifolia clone. Photosynthetica. 1983, 17: 537-549.
Weih M: Genetic and environmental variation in spring and autumn phenology of biomass willows (Salix spp.): effects on shoot growth and nitrogen economy. Tree Physiol. 2009, 29: 1479-1490. 10.1093/treephys/tpp081.
Pallardy SG: Physiology of Woody Plants. 2008, London, UK: Elseiver Inc.
Fuchigami LH, Weise CJ, Evert DR: Induction of cold acclimation in Cornus stolonifera Michx. Plant Physiol. 1971, 47: 98-103. 10.1104/pp.47.1.98.
Scarascia-Mugnozza GE, Hinckley TM, Stettler RF, Heilman PE, Isebrands JG: Production physiology and morphology of Populus species and their hybrids grown under short rotation. III. Seasonal carbon allocation patterns from branches. Can J Forest Res. 1999, 29: 1419-1432. 10.1139/cjfr-29-9-1419.
Dodd RS, Huberli D, Mayer W, Harnik TY, Afzal-Rafli Z, Garbelotto M: Evidence for the role of synchronicity between host phenology and pathogen activity in the distribution of sudden oak death canker disease. New Phytol. 2008, 179: 505-514. 10.1111/j.1469-8137.2008.02450.x.
Ghelardini L, Santini A: Avoidance by early flushing: a new perspective on Dutch elm disease. iForest. 2009, 2: 143-153. 10.3832/ifor0508-002.
McKinney LV, Nielsen LR, Hansen JK, Kjaer ED: Presence of natural genetic resistance in Fraxinus excelsior (Oleraceae) to Chalara fraxinea (Ascomycota): an emerging infectious disease. Heredity. 2011, 106: 788-797. 10.1038/hdy.2010.119.
Sinkkonen A, Somerkoski E, Paaso U, Holopainen JK, Rousi M, Mikola J: Genotypic variation in yellow autumn leaf colours explains aphid load in silver birch. New Phytol. 2012, 195: 461-469. 10.1111/j.1469-8137.2012.04156.x.
Lagercrantz U: At the end of the day: a common molecular mechanism for photoperiod responses in plants?. J Exp Bot. 2009, 60: 2501-2515. 10.1093/jxb/erp139.
Howe GT, Gardner G, Hackett WP, Furnier GR: Phytochrome control of short-day-induced bud set in black cottonwood. Physiol Plant. 1996, 97: 95-103. 10.1111/j.1399-3054.1996.tb00484.x.
Bradshaw HD, Stettler RF: Molecular genetics of growth and development in Populus. IV. Mapping QTLs with large effects on growth, form, and phenology traits in a forest tree. Genetics. 1995, 139: 963-973.
Keller SR, Soolanayakanahally RY, Guy RD, Silim SN, Olson MS, Tiffin P: Climate-driven local adaptation of ecophysiology and phenology in balsam poplar, Populus balsamifera L. (Salicaceae). Am J Bot. 2011, 98: 99-108. 10.3732/ajb.1000317.
Rohde A, Storme V, Jorge V, Gaudet M, Vitacolonna N, Fabbrini F, Ruttink T, Zaina G, Marron N, Dillen S, Steenackers M, Sabatti M, Morgante M, Boerjan W, Bastien C: Bud set in poplar - genetic dissection of a complex trait in natural and hybrid populations. New Phytol. 2011, 189: 106-121. 10.1111/j.1469-8137.2010.03469.x.
Weih M: Intensive short rotation forestry in boreal climates: present and future perspectives. Can J Forest Res. 2004, 34: 1369-1378. 10.1139/x04-090.
Tanino KK, Kalcsits L, Silim S, Kendall E, Gray GR: Temperature-driven plasticity in growth cessation and dormancy development in deciduous woody plants: a working hypothesis suggesting how molecular and cellular function is affected by temperature during dormancy induction. Plant Mol Biol. 2010, 73: 49-65. 10.1007/s11103-010-9610-y.
Rohde A, Bastien C, Boerjan W: Temperature signals contribute to the timing of photoperiodic growth cessation and bud set in poplar. Tree Physiol. 2011, 31: 472-482. 10.1093/treephys/tpr038.
Heide OM, Prestud AK: Low temperature, but not photoperiod, controls growth cessation and dormancy induction and relaese in apple and pear. Tree Physiol. 2005, 25: 109-114. 10.1093/treephys/25.1.109.
Heide OM: Interaction of photoperiod and temperature in the control of growth and dormancy of Prunus species. Sci Hortic-Amsterdam. 2008, 115: 309-314. 10.1016/j.scienta.2007.10.005.
Heide OM: Temperature rather than photoperiod controls growth cessation and dormancy in Sorbus species. J Exp Bot. 2011, 62: 5397-5404. 10.1093/jxb/err213.
Kozlowski TT, Pallardy SG: Growth Control in Woody Plants. 1997, San Diego, USA: Academic Press Inc
Fracheboud Y, Luquez V, Bjorken L, Sjodin A, Tuominen H, Jansson S: The control of autumn senescence in European aspen. Plant Physiol. 2009, 149: 1982-1991. 10.1104/pp.108.133249.
Keskitalo J, Bergquist G, Gardestrom P, Jansson S: A cellular timetable of autumn senescence. Plant Physiol. 2005, 139: 1635-1648. 10.1104/pp.105.066845.
van der Schoot C, Rinne PLH: Dormancy cycling at the shoot apical meristem: transitioning between self-organization and self-arrest. Plant Sci. 2011, 180: 120-131. 10.1016/j.plantsci.2010.08.009.
Perry TO: Dormancy of trees in winter. Science. 1971, 171: 29-36. 10.1126/science.171.3966.29.
Polgar CA, Primack RB: Leaf-out phenology of temperate woody plants: from trees to ecosystems. New Phytol. 2011, 191: 926-941. 10.1111/j.1469-8137.2011.03803.x.
Heide OM: Dormancy release in beech buds (Fagus sylvatica) requires both chilling and long days. Physiol Plant. 1993, 89: 187-191. 10.1111/j.1399-3054.1993.tb01804.x.
Sanz-Perez V, Castro-Diez P, Valladares F: Differential and interactive effects of temperature and photoperiod on budburst and carbon reserves in two co-occurring Mediterranean oaks. Plant Biol. 2009, 11: 142-151. 10.1111/j.1438-8677.2008.00119.x.
Boucher LD, Manchester SR, Judd WS: An extinct genus of salicaceae based on twigs with attached flowers fruits, and foliage from the eocene green river formation of Utah and Colorado, USA. Am J Bot. 2003, 90: 1389-1399. 10.3732/ajb.90.9.1389.
Manchester SR, Judd WS, Handley B: Foliage and fruits of early poplars (Salicaceae : Populus) from the eocene of Utah, Colorado, and Wyoming. Int J Plant Sci. 2006, 167: 897-908. 10.1086/503918.
Berlin S, Lagercrantz U, von Arnold S, öst T, Rönnberg-Wästljung AC: High-density linkage mapping and evolution of paralogs and orthologs in Salix and Populus. BMC Genomics. 2010, 11: 129-10.1186/1471-2164-11-129.
Kuzovkina Y, Weih M, Abalos Romero M, Belyaeva I, Charles J, Hurst S, Karp A, Labrecque M, McIvor I, Singh NB, Smart LB, Volk TA: Salix: botany and global horticulture. Hortic Revs. 2008, 34: 447-489.
Berlin S, Fogelqvist J, Lascoux M, Lagercrantz U, Rönnberg-Wästljung AC: Polymorphism and divergence in two willow species, Salix viminalis L. and Salix schwerinii E. Wolf. G3. 2011, 1: 387-400. 2011.
Junttila O: Apical growth cessation and shoot tip abscission in Salix. Physiol Plant. 1976, 38: 278-286. 10.1111/j.1399-3054.1976.tb04004.x.
Junttila O: Effect of photoperiod and temperature on apical growth cessation in two ecotypes of Salix and Betula. Physiol Plant. 1980, 48: 347-352. 10.1111/j.1399-3054.1980.tb03266.x.
Rönnberg-Wästljung AC, Gullberg U: Genetics of breeding characters with possible effects on biomass production in Salix viminalis (L.). Theor Appl Genet. 1999, 98: 531-540. 10.1007/s001220051101.
Lennartsson M, ögren E: Screening for efficient cold hardening in a breeding population of Salix using near infrared reflectance spectroscopy. Ann Forest Sci. 2004, 61: 449-454. 10.1051/forest:2004038.
Tsarouhas V, Gullberg U, Lagercrantz U: Mapping of quantitative trait loci controlling timing of bud flush in Salix. Hereditas. 2003, 138: 172-178. 10.1034/j.1601-5223.2003.01695.x.
Tsarouhas V, Gullberg U, Lagercrantz U: Mapping of quantitative trait loci (QTLs) affecting autumn freezing resistance and phenology in Salix. Theor Appl Genet. 2004, 108: 1335-1342. 10.1007/s00122-003-1544-1.
Pokhilko A, Fernandez AP, Edwards KD, Southern MM, Halliday KJ, Millar AJ: The clock gene circuit in Arabidopsis includes a repressilator with additional feedback loops. Mol Syst Biol. 2012, 8: 574-
Hsu CY, Adams JP, Kim H, No K, Ma C, Strauss SH, Drnevich J, Vandervelde L, Ellis JD, Rice BM, Wickett N, Gunter LE, Tuskan GA, Brunner AM, Page GP, Barakat A, Carlson JE, DePamphilis CW, Luthe DS, Yuceer C: FLOWERING LOCUS T duplication coordinates reproductive and vegetative growth in perennial poplar. Proc Natl Acad Sci USA. 2011, 108: 10756-10761. 10.1073/pnas.1104713108.
Ibanez C, Kozarewa I, Johansson M, Ogren E, Rohde A, Eriksson ME: Circadian clock components regulate entry and affect exit of seasonal dormancy as well as winter hardiness in Populus trees. Plant Physiol. 2010, 153: 1823-1833. 10.1104/pp.110.158220.
Dunlap JM, Stettler RF: Genetic variation and productivity of Populus trichocarpa and its hybrids.IX. Phenology and Melampsora rust incidence of native black cottonwood clones from four river valleys in Washington. Forest Ecol Manag. 1996, 87: 233-256. 10.1016/S0378-1127(96)03774-7.
Frewen BE, Chen TH, Howe GT, Davis J, Rohde A, Boerjan W, Bradshaw HD: Quantitative trait loci and candidate gene mapping of bud set and bud flush in Populus. Genetics. 2000, 154: 837-845.
Howe GT, Saruul P, Davis J, Chen THH: Quantitative genetics of bud phenology, frost damage, and winter survival in an F-2 family of hybrid poplars. Theor Appl Genet. 2000, 101: 632-642. 10.1007/s001220051525.
Fabbrini F, Gaudet M, Bastien C, Zaina G, Harfouche A, Beritognolo I, Marron N, Morgante M, Scarascia-Mugnozza G, Sabatti M: Phenotypic plasticity, QTL mapping and genomic characterization of bud set in black poplar. BMC Plant Biol. 2012, 12: 47-10.1186/1471-2229-12-47.
Marron N, Storme V, Dillen SY, Bastien C, Ricciotti L, Salani F, Sabatti M, Rae AM, Ceulemans R, Boerjan W: Genomic regions involved in productivity of two interspecific poplar families in Europe. 2. Biomass production and its relationships with tree architecture and phenology. Tree Genet Genomes. 2010, 6: 533-554. 10.1007/s11295-010-0270-5.
Barros RS, Neill SJ: Shoot growth in willow (Salix viminalis) in relation to abscisic acid, plant water status and photoperiod. Physiol Plant. 1987, 70: 708-712. 10.1111/j.1399-3054.1987.tb04328.x.
Mølmann JA, Asante DKA, Jensen JB, Krane MN, Ernstsen A, Junttila O, Olsen JE: Low night temperature and inhibition of gibberellin biosynthesis override phytochrome action and induce bud set and cold acclimation, but not dormancy in PHYA overexpressors and wild-type of hybrid aspen. Plant Cell Environ. 2005, 28: 1579-1588. 10.1111/j.1365-3040.2005.01395.x.
Kalcsits LA, Silim S, Tanino K: Warm temperature accelerates short photoperiod-induced growth cessation and dormancy induction in hybrid poplar (Populus x spp.). Trees. 2009, 23: 971-979. 10.1007/s00468-009-0339-7.
Chen THH, Howe GT, Bradshaw HD: Molecular genetic analysis of dormancy-related traits in poplars. Weed Sci. 2002, 50: 232-240. 10.1614/0043-1745(2002)050[0232:MGAODR]2.0.CO;2.
Böhlenius H, Huang T, Charbonnel-Campaa L, Brunner AM, Jansson S, Strauss SH, Nilsson O: CO/FT regulatory module controls timing of flowering and seasonal growth cessation in trees. Science. 2006, 312: 1040-1043. 10.1126/science.1126038.
Ruttink T, Arend M, Morreel K, Storme V, Rombauts S, Fromm J, Bhalerao RP, Boerjan W, Rohde A: A molecular timetable for apical bud formation and dormancy induction in poplar. Plant Cell. 2007, 19: 2370-2390. 10.1105/tpc.107.052811.
Horvath DP, Chao WS, Suttle JC, Thimmapuram J, Anderson JV: Transcriptome analysis identifies novel responses and potential regulatory genes involved in seasonal dormancy transitions of leafy spurge (Euphorbia esula L.). BMC Genomics. 2008, 9: 536-10.1186/1471-2164-9-536.
Yamane H, Ooka T, Jotatsu H, Hosaka Y, Sasaki R, Tao R: Expressional regulation of PpDAM5 and PpDAM6, peach (Prunus persica) dormancy-associated MADS-box genes, by low temperature and dormancy-breaking reagent treatment. J Exp Bot. 2011, 62: 3481-3488. 10.1093/jxb/err028.
Diaz-Riquelme J, Lijavetzky D, Martinez-Zapater JM, Carmona MJ: Genome-wide analysis of MIKCC-type MADS box genes in grapevine. Plant Physiol. 2009, 149: 354-369. 10.1104/pp.108.131052.
Li Z, Reighard GL, Abbott AG, Bielenberg DG: Dormancy-associated MADS genes from the EVG locus of peach [Prunus persica (L.) Batsch] have distinct seasonal and photoperiodic expression patterns. J Exp Bot. 2009, 60: 3521-3530. 10.1093/jxb/erp195.
Danfors B, Ledin S, Rosenqvist H: Short-rotation Willow Coppice. Growers’ Manual. 1997, Uppsala, Sweden: Swedish Institute of Agricultural and Environmental Engineering
Van Ooijen JW: MapQTL ® 6, Software for the Mapping of Quantitative Trait Loci in Experimental Populations. 2009, Wageningen, Netherlands: Kyazama B.V.
Lander ES, Botstein D: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121: 185-199.
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138: 963-971.
Lander E, Kruglyak L: Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat Genet. 1995, 11: 241-247. 10.1038/ng1195-241.
Van Ooijen JW: LOD significance thresholds for QTL analysis in experimental populations of diploid species. Heredity. 1999, 83: 613-624. 10.1038/sj.hdy.6886230.
Jansen RC, Stam P: High resolution of quantitative traits into multiple loci via interval mapping. Genetics. 1994, 136: 1447-1455.
Snijders TAB, Bosker RJ: Multilevel Analysis: An Introduction to Basic and Advanced Multilevel Modeling. 1999, London, UK: Sage Publications Ltd
R Development Core Team: R: A Language and Environment for Statistical Computing. 2010, Vienna, Austria: R Foundation for Statistical Computing
We thank Nils-Erik Nordh, Richard Childs and Lorenzo Bonosi for conducting phenotypic mesurements in the field; Marika Pettersson and Andreas Fritz for phenotypic measurements indoors. We also thank Pascal Pucholt for constructing the pearl script for the primer design and Francois Besnier for suggestions on the QTL analyses. This project was supported by funds from the Swedish Energy Agency and the NL-faculty at the Swedish University of Agricultural Sciences.
The authors declare that they have no competing interests.
LG performed phenotypic measurements, conducted and interpreted the phenotypic and QTL analyses and wrote the manuscript. SB created the new linkage map, performed the genome analysis and wrote the manuscript. MW read and commented on the manuscript and provided funding for measurements in the field experiment. NG and UL were involved in phenotypic measurements, selected candidate genes, read and commented on the manuscript. ACRW conceived the study, performed phenotypic measurements, conducted QTL analysis with the new linkage map, wrote the manuscript and provided funding for the project. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Pairwise plots of phenology traits across years and environments based on mean values for each individual in the mapping population S1. Red line show the linear fit of points. (PPTX 2 MB)
Additional file 4: Methods: Improved linkage map in two QTL hot spots for phenology on LG II and X. (DOCX 103 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Ghelardini, L., Berlin, S., Weih, M. et al. Genetic architecture of spring and autumn phenology in Salix. BMC Plant Biol 14, 31 (2014). https://doi.org/10.1186/1471-2229-14-31
- Candidate genes