Clear effects on root system architecture of winter wheat cultivars (Triticum aestivum L.) from cultivation environment and practices.

Roots play a pivotal role in the adaption of a plant to its environment, with different root traits adapting the plant to different stresses. The environment affects the Root System Architecture (RSA), but the genetic factors determine to what extent, and whether stress brought about by extreme environmental conditions is detrimental to a specific crop. This study aimed to identify differences in winter wheat RSA caused by cultivation region and practice, in the form of preceding crop (precrop), and to identify if modern cultivars used in Sweden differ in their reaction to these environments. This was undertaken using high-throughput phenotyping to assess the RSA. Clear differences in the RSA were observed between the Swedish cultivation regions, precrop treatments, and interaction of these conditions with each other and the genetics. Julius showed a large difference between cultivars, with 9.3-17.1% fewer and 12-20% narrower seminal roots. Standardized yield decreased when grown after wheat, 23% less compared to oilseed rape (OSR), and when grown in the Southern region, 14% less than the Central region. Additionally, correlations were shown between the root number, angle, and grain yield, with different root types being correlated depending on the precrop. Cultivars on the Swedish market show differences that can be adapted to the region-precrop combinations. The differences in precrop effect on RSA between regions show global implications and a need for further assessment. Correlations between RSA and yield, based on root-type × precrop, indicate different needs of the RSA depending on the management practices and show the potential for improving crop yield through targeting genotypic and environmental conditions in a holistic manner. Understanding this RSA variance, and the mechanisms of conditional response, will allow targeted cultivar breeding for specific environments, increasing plant health and food security.

Tree diversity enhances predation by birds but not by arthropods across climate gradients.

Tree diversity can promote both predator abundance and diversity. However, whether this translates into increased predation and top-down control of herbivores across predator taxonomic groups and contrasting environmental conditions remains unresolved. We used a global network of tree diversity experiments (TreeDivNet) spread across three continents and three biomes to test the effects of tree species richness on predation across varying climatic conditions of temperature and precipitation. We recorded bird and arthropod predation attempts on plasticine caterpillars in monocultures and tree species mixtures. Both tree species richness and temperature increased predation by birds but not by arthropods. Furthermore, the effects of tree species richness on predation were consistent across the studied climatic gradient. Our findings provide evidence that tree diversity strengthens top-down control of insect herbivores by birds, underscoring the need to implement conservation strategies that safeguard tree diversity to sustain ecosystem services provided by natural enemies in forests.

Large effect of phosphate-solubilizing bacteria on the growth and gene expression of Salix spp. at low phosphorus levels.

Phosphorus is one of the most important nutrients required for plant growth and development. However, owing to its low availability in the soil, phosphorus is also one of the most difficult elements for plants to acquire. Phosphorus released into the soil from bedrock quickly becomes unavailable to plants, forming poorly soluble complexes. Phosphate-solubilizing bacteria (PSB) can solubilize unavailable phosphorus-containing compounds into forms in which phosphorus is readily available, thus promoting plant growth. In this study, two willow species, Salix dasyclados cv. Loden and Salix schwerinii × Salix viminalis cv. Tora, were inoculated with two selected bacterial strains, Pantoea agglomerans and Paenibacillus spp., to evaluate the plant growth parameters and changes in gene expression in the presence of different concentrations of tricalcium phosphate: 0 mM (NP), 1 mM (LP), and 2 mM (HP). Inoculation with PSB increased root, shoot and leaf biomass, and for the HP treatment, significant changes in growth patterns were observed. However, the growth responses to plant treatments tested depended on the willow species. Analysis of the leaf transcriptomes of the phosphate-solubilizing bacterium-inoculated plants showed a large variation in gene expression between the two willow species. For the Tora willow species, upregulation of genes was observed, particularly for those involved in pathways related to photosynthesis, and this effect was strongly influenced by bacterial phosphate solubilization. The Loden willow species was characterized by a general downregulation of genes involved in pathway activity that included ion transport, transcription regulation and chromosomes. The results obtained in this study provide an improved understanding of the dynamics of Salix growth and gene expression under the influence of PSB, contributing to an increase in yield and phosphorus-use efficiency.

Disentangling drivers of litter decomposition in a multi-continent network of tree diversity experiments.

Litter decomposition is a key ecosystem function in forests and varies in response to a range of climatic, edaphic, and local stand characteristics. Disentangling the relative contribution of these factors is challenging, especially along large environmental gradients. In particular, knowledge of the effect of management options, such as tree planting density and species composition, on litter decomposition would be highly valuable in forestry. In this study, we made use of 15 tree diversity experiments spread over eight countries and three continents within the global TreeDivNet network. We evaluated the effects of overstory composition (tree identity, species/mixture composition and species richness), plantation conditions (density and age), and climate (temperature and precipitation) on mass loss (after 3 months and 1 year) of two standardized litters: high-quality green tea and low-quality rooibos tea. Across continents, we found that early-stage decomposition of the low-quality rooibos tea was influenced locally by overstory tree identity. Mass loss of rooibos litter was higher under young gymnosperm overstories compared to angiosperm overstories, but this trend reversed with age of the experiment. Tree species richness did not influence decomposition and explained almost no variation in our multi-continent dataset. Hence, in the young plantations of our study, overstory composition effects on decomposition were mainly driven by tree species identity on decomposer communities and forest microclimates. After 12 months of incubation, mass loss of the high-quality green tea litter was mainly influenced by temperature whereas the low-quality rooibos tea litter decomposition showed stronger relationships with overstory composition and stand age. Our findings highlight that decomposition dynamics are not only affected by climate but also by management options, via litter quality of the identity of planted trees but also by overstory composition and structure.

Mixed growth of Salix species can promote phosphate-solubilizing bacteria in the roots and rhizosphere.

Phosphorus (P) is an essential plant nutrient that can limit plant growth due to low availability in the soil. P-solubilizing bacteria in the roots and rhizosphere increase the P use efficiency of plants. This study addressed the impact of plant species, the level of plant association with bacteria (rhizosphere or root endophyte) and environmental factors (e.g., seasons, soil properties) on the abundance and diversity of P-solubilizing bacteria in short-rotation coppices (SRC) of willows (Salix spp.) for biomass production. Two willow species (S. dasyclados cv. Loden and S. schwerinii × S. viminalis cv. Tora) grown in mono-and mixed culture plots were examined for the abundance and diversity of bacteria in the root endosphere and rhizosphere during two seasons (fall and spring) in central Sweden and northern Germany. Soil properties, such as pH and available P and N, had a significant effect on the structure of the bacterial community. Microbiome analysis and culture-based methods revealed a higher diversity of rhizospheric bacteria than endophytic bacteria. The P-solubilizing bacterial isolates belonged mainly to Proteobacteria (85%), Actinobacteria (6%) and Firmicutes (9%). Pseudomonas was the most frequently isolated cultivable bacterial genus from both the root endosphere and the rhizosphere. The remaining cultivable bacterial isolates belonged to the phyla Actinobacteria and Firmicutes. In conclusion, site-specific soil conditions and the level of plant association with bacteria were the main factors shaping the bacterial communities in the willow SRCs. In particular, the concentration of available P along with the total nitrogen in the soil controlled the total bacterial diversity in willow SRCs. A lower number of endophytic and rhizospheric bacteria was observed in Loden willow species compared to that of Tora and the mix of the two, indicating that mixed growth of Salix species promotes P-solubilizing bacterial diversity and abundance. Therefore, a mixed plant design was presented as a management option to increase the P availability for Salix in SRCs. This design should be tested for further species mixtures.

Temporal complementarity between roots and mycorrhizal fungi drives wheat nitrogen use efficiency.

Improving nitrogen (N) use efficiency (NUE) to reduce the application of N fertilisers in a way that benefits the environment and reduces farmers' costs is an ongoing objective for sustainable wheat production. However, whether and how arbuscular mycorrhizal fungi (AMF) affect NUE in wheat is still not well explored. Three independent but complementary experiments were conducted to decipher the contribution of roots and AMF to the N uptake and utilisation efficiency in wheat. We show a temporal complementarity pattern between roots and AMF in shaping NUE of wheat. Pre-anthesis N uptake efficiency mainly depends on root functional traits, but the efficiency to utilise the N taken up during pre-anthesis for producing grains (EN,g ) is strongly affected by AMF, which might increase the uptake of phosphorus and thereby improve photosynthetic carbon assimilation. Root association with AMF reduced the N remobilisation efficiency in varieties with high EN,g ; whilst the overall grain N concentration increased, due to a large improvement in post-anthesis N uptake supported by AMF and/or other microbes. The findings provide evidence for the importance of managing AMF in agroecosystems, and an opportunity to tackle the contradiction between maximising grain yield and protein concentration in wheat breeding.

Nutrient Accumulation Pattern in Mixtures of Wheat and Faba Bean Is Strongly Influenced by Cultivar Choice and Co-Existing Weeds.

Cereal-legume mixtures are often associated with higher yields than the components grown as sole crops, but the underlying mechanisms are unclear. The study aims to evaluate how different cultivars in a two-species wheat-faba bean mixture influence above- and below-ground nitrogen (N) accumulation in the plant biomass, whether crop mixing affected the accumulation of other nutrients relative to the accumulation of N and phosphorus (P), and how the nutrient accumulation pattern in sole crops and mixtures is influenced by weed competition. Using a growth container experiment, we investigate nutrient accumulation patterns on specific wheat and faba bean cultivars grown as sole crops and mixtures, and with and without weed competition. We found that cereals in the mixture accumulated more N than in the sole crops, and the cultivar used influenced biomass accumulation in the legumes. Competition from weeds reduced the amount of plant N pools accumulated in the crop plant biomass. Based on stoichiometric scaling exponents, the plant neighbor affected the accumulation of other nutrients relative to the accumulation of N and P. These results are relevant for species and cultivar selection, all of which are important prerequisites for maximizing mixture performance.

Supply Chain Perspectives on Breeding for Legume-Cereal Intercrops.

Compared to sole crops, intercropping-especially of legumes and cereals-has great potential to improve crop yield and resource use efficiency, and can provide many other ecosystem services. However, the beneficial effects of intercrops are often greatly dependent on the end use as well as the specific species and genotypes being co-cultivated. In addition, intercropping imposes added complexity at different levels of the supply chain. While the need for developing crop genotypes for intercropping has long been recognized, most cultivars on the market are optimized for sole cropping and may not necessarily perform well in intercrops. This paper aims to place breeding targets for intercrop-adapted genotypes in a supply chain perspective. Three case studies of legumes and cereals intercropped for human consumption are used to identify desirable intercrop traits for actors across the supply chains, many of which are not targeted by traditional breeding for sole crops, including certain seed attributes, and some of which do not fit traditional breeding schemes, such as breeding for synchronized maturity and species synergies. Incorporating these traits into intercrop breeding could significantly reduce complexity along the supply chain. It is concluded that the widespread adoption and integration of intercrops will only be successful through the inclusion and collaboration of all supply chain actors, the application of breeding approaches that take into account the complexity of intercrop supply chains, and the implementation of diversification strategies in every process from field to fork.

Application of Crop Growth Models to Assist Breeding for Intercropping: Opportunities and Challenges.

Intercropping of two or more species on the same piece of land can enhance biodiversity and resource use efficiency in agriculture. Traditionally, intercropping systems have been developed and improved by empirical methods within a specific local context. To support the development of promising intercropping systems, the individual species that are part of an intercrop can be subjected to breeding. Breeding for intercropping aims at resource foraging traits of the admixed species to maximize niche complementarity, niche facilitation, and intercrop performance. The breeding process can be facilitated by modeling tools that simulate the outcome of the combination of different species' (or genotypes') traits for growth and yield development, reducing the need of extensive field testing. Here, we revisit the challenges associated with breeding for intercropping, and give an outlook on applying crop growth models to assist breeding for intercropping. We conclude that crop growth models can assist breeding for intercropping, provided that (i) they incorporate the relevant plant features and mechanisms driving interspecific plant-plant interactions; (ii) they are based on model parameters that are closely linked to the traits that breeders would select for; and (iii) model calibration and validation is done with field data measured in intercrops. Minimalist crop growth models are more likely to incorporate the above elements than comprehensive but parameter-intensive crop growth models. Their lower complexity and reduced parameter requirement facilitate the exploration of mechanisms at play and fulfil the model requirements for calibration of the appropriate crop growth models.

Effects of soil compaction on grain yield of wheat depend on weather conditions.

The use of heavy farm machinery has resulted in widespread soil compaction in many regions of the world. Compacted soil limits the access of crops to soil water and nutrients and is expected to reduce crop productivity, but the influence of weather conditions on the interactions between compacted soil and crop productivity is unclear. Furthermore, early vigor has been regarded as a promising trait for improving the yield of crops grown under edaphic stress such as soil compaction. We aimed to assess the combined effects of soil compaction and contrasting weather conditions on growth and grain yield of spring wheat, and to evaluate the association between early vigor and grain yield under temporal variations of the soil physical conditions. Nine spring wheat genotypes were grown on compacted and non-compacted soils during two cropping seasons with contrasting weather conditions in Central Sweden. Compared to the non-compacted treatment, soil compaction increased the relative growth rate of shoot biomass from sowing to stem elongation, and from stem elongation to flowering in the drier year (2018), but decreased the same traits in the wetter year (2019). The contrasting effects of soil compaction on shoot growth in the two years could be explained by soil moisture and penetration resistance associated with the interactive effects of soil compaction and weather condition. Higher early vigor, here indicated by higher relative growth rate from sowing to stem elongation, was associated with reduced grain yield under the progressively drying and hardening soil conditions during the entire cropping season of both years. We conclude that the interactive effects of soil physical and weather conditions need to be considered when evaluating the impact of soil compaction on crop growth and productivity. The potential of early vigor to increase grain yield is strongly influenced by the temporal dynamics of soil physical conditions.

The Effects of Host Plant Genotype and Environmental Conditions on Fungal Community Composition and Phosphorus Solubilization in Willow Short Rotation Coppice.

Phosphorus (P) is an essential plant nutrient. Low availability of P in soil is mainly caused by high content of Fe2O3 in the clay fraction that binds to P making it unavailable. Beneficial microbes, such as P solubilizing microorganisms can increase the available P in soil and improve plant growth and productivity. In this study, we evaluated the effects of environmental conditions (climate, soil parameters), plant genotype, and level of plant association (rhizosphere or endophytic root organism) on the abundance and diversity of phosphorus solubilizing microorganisms in a Salix production system. We hypothesized that a lower number of endophytic fungi may possess the ability to solubilize P compared to the number of rhizosphere fungi with the same ability. We also expect that the plant genotype and the experimental site with its environmental conditions will influence fungal diversity. Two Salix genotypes grown in pure and mixed cultures were investigated for their fungal microbiome community and diversity in the rhizosphere and endosphere during two growing seasons. We found that the rhizosphere fungal community was more diverse. A general dominance of Ascomycota (Dothideomycetes) and Basidiomycota (Tremellomycetes) was observed. The classes Agaricomycetes and Pezizomycetes were more frequent in the endosphere, while Tremellomycetes and Mortierellomycetes were more abundant in the rhizosphere. Plot-specific soil properties (pH, total organic carbon, and nitrogen) significantly influenced the fungal community structure. Among the culturable fungal diversities, 10 strains of phosphate solubilizing fungi (PSFs) from roots and 12 strains from rhizosphere soil were identified using selective media supplemented with di-calcium and tri-calcium phosphates. The fungal density and the number of PSF were much higher in the rhizosphere than in the endosphere. Penicillium was the dominant genus of PSF isolated from both sites; other less frequent genera of PSFs were Alternaria, Cladosporium, and Clonostachys. Overall the main factors controlling the fungal communities (endophytic vs. rhizosphere fungi) were the soil properties and level of plant association, while no significant influence of growing season was observed. Differences between Salix genotypes were observed for culturable fungal diversity, while in metagenomic data analysis, only the class Dothideomycetes showed a significant effect from the plant genotype.

Evidence for magnesium-phosphorus synergism and co-limitation of grain yield in wheat agriculture.

Modern crop production is characterized by high nitrogen (N) application rates, which can influence the co-limitation of harvested yield by other nutrients. Using a multidimensional niche volume concept and scaling exponents frequently applied in plant ecological research, we report that increased N and phosphorus (P) uptake in a growing wheat crop along with enhanced grain biomass is associated with more than proportional increase of other nutrients. Furthermore, N conversion efficiency and grain yield are strongly affected by the magnesium (Mg) to P ratio in the growing crop. We analyzed a field trial in Central Sweden including nine wheat varieties grown during two years with contrasting weather, and found evidence for Mg co-limitation at lower grain yields and P co-limitation at higher yields. We argue that critical concentrations of single nutrients, which are often applied in agronomy, should be replaced by nutrient ratios. In addition, links between plant P and Mg contents and root traits were found; high root number enhanced the P:N ratio, whilst steep root angle, indicating deep roots, increased the Mg:N ratio. The results have significant implications on the management and breeding targets of agriculturally grown wheat, which is one of the most important food crops worldwide.

Shoot and Root Traits Underlying Genotypic Variation in Early Vigor and Nutrient Accumulation in Spring Wheat Grown in High-Latitude Light Conditions.

Plants with improved nutrient use efficiency are needed to maintain and enhance future crop plant production. The aim of this study was to explore candidate traits for pre-breeding to improve nutrient accumulation and early vigor of spring wheat grown at high latitudes. We quantified shoot and root traits together with nutrient accumulation in nine contrasting spring wheat genotypes grown in rhizoboxes for 20 days in a greenhouse. Whole-plant relative growth rate was here correlated with leaf area productivity and plant nitrogen productivity, but not leaf area ratio. Furthermore, the total leaf area was correlated with the accumulation of six macronutrients, and could be suggested as a candidate trait for the pre-breeding towards improved nutrient accumulation and early vigor in wheat to be grown in high-latitude environments. Depending on the nutrient of interest, different root system traits were identified as relevant for their accumulation. Accumulation of nitrogen, potassium, sulfur and calcium was correlated with lateral root length, whilst accumulation of phosphorus and magnesium was correlated with main root length. Therefore, special attention needs to be paid to specific root system traits in the breeding of wheat towards improved nutrient accumulation to counteract the suboptimal uptake of some nutrient elements.

Site-Dependent Relationships Between Fungal Community Composition, Plant Genotypic Diversity and Environmental Drivers in a Salix Biomass System.

Soil fungi are strongly affected by plant species or genotypes since plants modify their surrounding environment, but the effects of plant genotype diversity on fungal diversity and function have not been extensively studied. The interactive responses of fungal community composition to plant genotypic diversity and environmental drivers were investigated in Salix biomass systems, posing questions about: (1) How fungal diversity varies as a function of plant genotype diversity; (2) If plant genotype identity is a strong driver of fungal community composition also in plant mixtures; (3) How the fungal communities change through time (seasonally and interannually)?; and (4) Will the proportion of ECM fungi increase over the rotation? Soil samples were collected over 4 years, starting preplanting from two Salix field trials, including four genotypes with contrasting phenology and functional traits, and genotypes were grown in all possible combinations (four genotypes in Uppsala, Sweden, two in Rostock, Germany). Fungal communities were identified, using Pacific Biosciences sequencing of fungal ITS2 amplicons. We found some site-dependent relationships between fungal community composition and genotype or diversity level, and site accounted for the largest part of the variation in fungal community composition. Rostock had a more homogenous community structure, with significant effects of genotype, diversity level, and the presence of one genotype ("Loden") on fungal community composition. Soil properties and plant and litter traits contributed to explaining the variation in fungal species composition. The within-season variation in composition was of a similar magnitude to the year-to-year variation. The proportion of ECM fungi increased over time irrespective of plant genotype diversity, and, in Uppsala, the 4-mixture showed a weaker response than other combinations. Species richness was generally higher in Uppsala compared with that in Rostock and increased over time, but did not increase with plant genotype diversity. This significant site-specificity underlines the need for consideration of diverse sites to draw general conclusions of temporal variations and functioning of fungal communities. A significant increase in ECM colonization of soil under the pioneer tree Salix on agricultural soils was evident and points to changed litter decomposition and soil carbon dynamics during Salix growth.

The Arabidopsis Transcription Factor CDF3 Is Involved in Nitrogen Responses and Improves Nitrogen Use Efficiency in Tomato.

Nitrate is an essential macronutrient and a signal molecule that regulates the expression of multiple genes involved in plant growth and development. Here, we describe the participation of Arabidopsis DNA binding with one finger (DOF) transcription factor CDF3 in nitrate responses and shows that CDF3 gene is induced under nitrate starvation. Moreover, knockout cdf3 mutant plants exhibit nitrate-dependent lateral and primary root modifications, whereas CDF3 overexpression plants show increased biomass and enhanced root development under both nitrogen poor and rich conditions. Expression analyses of 35S::CDF3 lines reveled that CDF3 regulates the expression of an important set of nitrate responsive genes including, glutamine synthetase-1, glutamate synthase-2, nitrate reductase-1, and nitrate transporters NRT2.1, NRT2.4, and NRT2.5 as well as carbon assimilation genes like PK1 and PEPC1 in response to N availability. Consistently, metabolite profiling disclosed that the total amount of key N metabolites like glutamate, glutamine, and asparagine were higher in CDF3-overexpressing plants, but lower in cdf3-1 in N limiting conditions. Moreover, overexpression of CDF3 in tomato increased N accumulation and yield efficiency under both optimum and limiting N supply. These results highlight CDF3 as an important regulatory factor for the nitrate response, and its potential for improving N use efficiency in crops.

Quantitative genetic architecture of adaptive phenology traits in the deciduous tree, Populus trichocarpa (Torr. and Gray).

In a warming climate, the ability to accurately predict and track shifting environmental conditions will be fundamental for plant survival. Environmental cues define the transitions between growth and dormancy as plants synchronise development with favourable environmental conditions, however these cues are predicted to change under future climate projections which may have profound impacts on tree survival and growth. Here, we use a quantitative genetic approach to estimate the genetic basis of spring and autumn phenology in Populus trichocarpa to determine this species capacity for climate adaptation. We measured bud burst, leaf coloration, and leaf senescence traits across two years (2017-2018) and combine these observations with measures of lifetime growth to determine how genetic correlations between phenology and growth may facilitate or constrain adaptation. Timing of transitions differed between years, although we found strong cross year genetic correlations in all traits, suggesting that genotypes respond in consistent ways to seasonal cues. Spring and autumn phenology were correlated with lifetime growth, where genotypes that burst leaves early and shed them late had the highest lifetime growth. We also identified substantial heritable variation in the timing of all phenological transitions (h2 = 0.5-0.8) and in lifetime growth (h2 = 0.8). The combination of additive variation and favourable genetic correlations in phenology traits suggests that populations of cultivated varieties of P. Trichocarpa may have the capability to adapt their phenology to climatic changes without negative impacts on growth.

Corrigendum: Multi-Dimensional Plant Element Stoichiometry-Looking Beyond Carbon, Nitrogen, and Phosphorus.

[This corrects the article DOI: 10.3389/fpls.2020.00023.].

Multi-Dimensional Plant Element Stoichiometry-Looking Beyond Carbon, Nitrogen, and Phosphorus.

Nutrient elements are important for plant growth. Element stoichiometry considers the balance between different nutrients and how this balance is affected by the environment. So far, focus of plant stoichiometry has mainly been on the three elements carbon (C), nitrogen (N), and phosphorus (P), but many additional elements are essential for proper plant growth. Our overall aim is to test the scaling relations of various additional elements (K, Ca, Mg, S, Cu, Zn, Fe, Mn), by using ten data sets from a range of plant functional types and environmental conditions. To simultaneously handle more than one element, we define a stoichiometric niche volume as the volume of an abstract multidimensional shape in n dimensions, with the n sides of this shape defined by the plant properties in question, here their element concentrations. Thus, a stoichiometric niche volume is here defined as the product of element concentrations. The volumes of N and P (VNP ) are used as the basis, and we investigate how the volume of other elements (VOth ) scales with respect to VNP¸ with the intention to explore if the concentrations of other elements increase faster (scaling exponent > 1) or slower (<1) than the concentrations of N and P. For example, scaling exponents >1 suggest that favorable conditions for plant growth, i.e., environments rich in N and P, may require proportionally higher uptake of other essential elements than poor conditions. We show that the scaling exponent is rather insensitive to environmental conditions or plant species, and ranges from 0.900 to 2.479 (average 1.58) in nine out of ten data sets. For single elements, Mg has the smallest scaling exponent (0.031) and Mn the largest (2.147). Comparison between laboratory determined stoichiometric relations and field observations suggest that element uptake in field conditions often exceeds the minimal physiological requirements. The results provide evidence for the view that the scaling relations previously reported for N and P can be extended to other elements; and that N and P are the driving elements in plant stoichiometric relations. The stoichiometric niche volumes defined here could be used to predict plant performances in different environments.

Genome Wide Associations of Growth, Phenology, and Plasticity Traits in Willow [Salix viminalis (L.)].

The short rotation biomass crop willow (Salix genera) has been of interest for bioenergy but recently also for biofuel production. For a faster development of new varieties molecular markers could be used as selection tool in an early stage of the breeding cycle. To identify markers associated with growth traits, genome-wide association mapping was conducted using a population of 291 Salix viminalis accessions collected across Europe and Russia and a large set of genotyping-by-sequencing markers. The accessions were vegetatively propagated and planted in replicated field experiments, one in Southern Sweden and one in Central Sweden. Phenology data, including bud burst and leaf senescence, as well as different growth traits were collected and measured repeatedly between 2010 and 2017 at both field environments. A value of the plasticity for each accession was calculated for all traits that were measured the same year in both environments as the normalized accession value in one environment subtracted by the corresponding value in the other environment. Broad-sense accession heritabilities and narrow-sense chip heritabilities ranged from 0.68 to 0.95 and 0.45 to 0.99, respectively for phenology traits and from 0.56 to 0.85 and 0.24 to 0.97 for growth traits indicating a considerable genetic component for most traits. Population structure and kinship between accessions were taken into account in the association analyses. In total, 39 marker-trait associations were found where four were specifically connected to plasticity and interestingly one particular marker was associated to several different plasticity growth traits. Otherwise association consistency was poor, possibly due to accession by environment interactions which were demonstrated by the low structure adjusted accession correlations across environments (ranging from 0.40 to 0.58). However, one marker association with biomass fresh weight was repeatedly observed in the same environment over two harvest years. For some traits where several associations were found, the markers jointly explained over 20% of the accession variation. The result from this study using a population of unrelated accessions has given useful information about marker-trait associations especially highlighting marker-plasticity associations and genotype-by-environment interactions as important factors to take account of in future strategies of Salix breeding.