QTL analyses of drought tolerance and growth for a Salix dasyclados x Salix viminalis hybrid in contrasting water regimes

Mapping the sex determination region in the Salix F1 hybrid common parent population confirms a ZW system in six diverse species

Within the genus Salix, there are approximately 350 species native primarily to the northern hemisphere and adapted to a wide range of habitats. This diversity can be exploited to mine novel alleles conferring variation important for production as a bioenergy crop, but also to identify evolutionarily important genes, such as those involved in sex determination. To leverage this diversity, we created a mapping population by crossing 6 Salix species (Salix viminalis, Salix suchowensis, Salix integra, Salix koriyanagi, Salix udensis, and Salix alberti) to common male and female Salix purpurea parents. Each family was genotyped via genotyping-by-sequencing and assessed for kinship and population structure as well as the construction of 16 backcross linkage maps to be used as a genetic resource for breeding and selection. Analyses of population structure resolved both the parents and F₁ progeny to their respective phylogenetic section and indicated that the S. alberti parent was misidentified and was most likely S.suchowensis. Sex determining regions were identified on Salix chromosome 15 in the female-informative maps for seven of the eight families indicating that these species share a common female heterogametic ZW sex system. The eighth family, S. integra × S. purpurea, was entirely female and had a truncated chromosome 15. Beyond sex determination, the Salix F₁ hybrid common parent population (Salix F₁ HCP) introduced here will be useful in characterizing genetic factors underlying complex traits, aid in marker-assisted selection, and support genome assemblies for this promising bioenergy crop.

QTL mapping of drought-related traits in the hybrids of Populus deltoides 'Danhong'×Populus simonii 'Tongliao1'

BACKGROUND

Poplar trees provide a large amount of wood material, but many parts of the world are arid or semi-arid areas because of insufficient annual precipitation, which seriously affects the growth of poplar trees. Populus simonii 'Tongliao1' shows strong tolerance to stress environments, and Populus deltoides 'Danhong' shows a stronger growth rate in a suitable environment. To identify drought tolerance-related QTLs and genes, an F₁ population derived from the cross between the 'Danhong' and 'Tongliao 1' Populus was assessed under drought stress.

RESULTS

We measured drought-related traits such as the relative height growth, relative diameter growth, leaf senescence number, specific leaf area, and leaf relative water content in the population under control and drought environments. The results showed that drought stress reduced the plant height relative growth, ground diameter relative growth, specific leaf area and leaf relative water content and increased the number of leaf drops. A total of 208 QTLs were identified by QTL mapping analysis, and they consisted of 92, 63 and 53 QTLs under control, drought stress treatment and drought index conditions, respectively. A molecular identification marker for drought tolerance, np2841, which was associated with a QTL (qDLRWC-LG10-1) for relative leaf water content, was initially developed. We mined 187 candidate genes for QTL regions of five traits under a drought environment. The reference genome annotation for Populus trichocarpa and a homologous gene analysis of Arabidopsis thaliana identified two candidate genes, Potri.003G171300 and Potri.012G123900, with significant functions in response to drought stress. We identified five key regulatory genes (Potri.006G273500, Potri.007G111500, Potri.007G111600, Potri.007G111700, and Potri.007G111800) related to drought tolerance through the poplar coexpression network.

CONCLUSION

In this study, our results indicate that the QTLs can effectively enhance the drought tolerance of poplar. It is a step closer towards unravelling the genetic basis of poplar drought tolerance-related traits, and to providing validated candidate genes and molecular markers for future genetic improvement.

Time-sequential detection of quantitative trait loci and candidate genes underlying the dynamic growth of Salix suchowensis

Elucidating the genetic factors underlying long-term biological processes remains challenging since the relevant genes and their effects may vary across different developmental stages. In this study, we carried out a large-scale field trial of the progeny of an F1 full-sib pedigree of Salix suchowensis and measured plant height and ground diameter periodically over a time course of 240 days. With the obtained data, we characterized plant growth rhythms and performed time-sequential analyses of quantitative trait loci underlying the dynamic growth of the plants. The dynamic mapping of quantitative trait loci revealed that stem height and ground diameter were under the control of four quantitative trait loci, and the effects of these quantitative trait loci varied greatly throughout the growth process, in which two quantitative trait loci were found to exert a pleiotropic effect determining the correlation between stem height and ground diameter. The analysis of candidate genes in the target genetic intervals showed that the pleiotropic effect of the two quantitative trait loci arises from the colocalization of genes with independent effects on stem height and ground diameter. Further examination of the expression patterns of the candidate genes indicated that height and circumference growth involve different activities of leaf and cambium tissues. This study provides unprecedented information to help us understand the dynamic growth of plants and presents an applicable strategy for elucidating the genetic mechanism underlying a long-term biological process by using plant growth as an example.

QTL mapping for growth-related traits by constructing the first genetic linkage map in Simao pine

BACKGROUND

Simao pine is one of the primary economic tree species for resin and timber production in southwest China. The exploitation and utilization of Simao pine are constrained by the relatively lacking of genetic information. Construction a fine genetic linkage map and detecting quantitative trait locis (QTLs) for growth-related traits is a prerequisite section of Simao Pine's molecular breeding program.

RESULTS

In our study, a high-resolution Simao pine genetic map employed specific locus amplified fragment sequencing (SLAF-seq) technology and based on an F₁ pseudo-testcross population has been constructed. There were 11,544 SNPs assigned to 12 linkage groups (LGs), and the total length of the map was 2,062.85 cM with a mean distance of 0.37 cM between markers. According to the phenotypic variation analysis for three consecutive years, a total of seventeen QTLs for four traits were detected. Among 17 QTLs, there were six for plant height (Dh.16.1, Dh16.2, Dh17.1, Dh18.1-3), five for basal diameter (Dbd.17.1-5), four for needle length (Dnl17.1-3, Dnl18.1) and two for needle diameter (Dnd17.1 and Dnd18.1) respectively. These QTLs individually explained phenotypic variance from 11.0-16.3%, and the logarithm of odds (LOD) value ranged from 2.52 to 3.87.

CONCLUSIONS

In our study, a fine genetic map of Simao pine applied the technology of SLAF-seq has been constructed for the first time. Based on the map, a total of 17 QTLs for four growth-related traits were identified. It provides helpful information for genomic studies and marker-assisted selection (MAS) in Simao pine.

Genome-wide association study of soybean seed germination under drought stress

Drought stress, which is increasing with climate change, is a serious threat to agricultural sustainability worldwide. Seed germination is an essential growth phase that ensures the successful establishment and productivity of soybean, which can lose substantial productivity in soils with water deficits. However, only limited genetic information is available about how germinating soybean seeds may exert drought tolerance. In this study, we examined the germinating seed drought-tolerance phenotypes and genotypes of a panel of 259 released Chinese soybean cultivars panel. Based on 4616 Single-Nucleotide Polymorphisms (SNPs), we conducted a mixed-linear model GWAS that identified a total of 15 SNPs associated with at least one drought-tolerance index. Notably, three of these SNPs were commonly associated with two drought-tolerance indices. Two of these SNPs are positioned upstream of genes, and 11 of them are located in or near regions where QTLs have been previously mapped by linkage analysis, five of which are drought-related. The SNPs detected in this study can both drive hypothesis-driven research to deepen our understanding of genetic basis of soybean drought tolerance at the germination stage and provide useful genetic resources that can facilitate the selection of drought stress traits via genomic-assisted selection.

Marginal Lands to Grow Novel Bio-Based Crops: A Plant Breeding Perspective

The biomass demand to fuel a growing global bio-based economy is expected to tremendously increase over the next decades, and projections indicate that dedicated biomass crops will satisfy a large portion of it. The establishment of dedicated biomass crops raises huge concerns, as they can subtract land that is required for food production, undermining food security. In this context, perennial biomass crops suitable for cultivation on marginal lands (MALs) raise attraction, as they could supply biomass without competing for land with food supply. While these crops withstand marginal conditions well, their biomass yield and quality do not ensure acceptable economic returns to farmers and cost-effective biomass conversion into bio-based products, claiming genetic improvement. However, this is constrained by the lack of genetic resources for most of these crops. Here we first review the advantages of cultivating novel perennial biomass crops on MALs, highlighting management practices to enhance the environmental and economic sustainability of these agro-systems. Subsequently, we discuss the preeminent breeding targets to improve the yield and quality of the biomass obtainable from these crops, as well as the stability of biomass production under MALs conditions. These targets include crop architecture and phenology, efficiency in the use of resources, lignocellulose composition in relation to bio-based applications, and tolerance to abiotic stresses. For each target trait, we outline optimal ideotypes, discuss the available breeding resources in the context of (orphan) biomass crops, and provide meaningful examples of genetic improvement. Finally, we discuss the available tools to breed novel perennial biomass crops. These comprise conventional breeding methods (recurrent selection and hybridization), molecular techniques to dissect the genetics of complex traits, speed up selection, and perform transgenic modification (genetic mapping, QTL and GWAS analysis, marker-assisted selection, genomic selection, transformation protocols), and novel high-throughput phenotyping platforms. Furthermore, novel tools to transfer genetic knowledge from model to orphan crops (i.e., universal markers) are also conceptualized, with the belief that their development will enhance the efficiency of plant breeding in orphan biomass crops, enabling a sustainable use of MALs for biomass provision.

Joint linkage and association mapping of complex traits in shrub willow (Salix purpurea L.)

BACKGROUND AND AIMS

Increasing energy demands and the necessity to reduce greenhouse gas emissions are key motivating factors driving the development of lignocellulosic crops as an alternative to non-renewable energy sources. The effects of global climate change will require a better understanding of the genetic basis of complex adaptive traits to breed more resilient bioenergy feedstocks, like willow (Salix spp.). Shrub willow is a sustainable and dedicated bioenergy crop, bred to be fast-growing and high-yielding on marginal land without competing with food crops. In a rapidly changing climate, genomic advances will be vital for the sustained improvement of willow and other non-model bioenergy crops. Here, joint genetic mapping was used to exploit genetic variation garnered from both recent and historical recombination events in S. purpurea.

METHODS

A panel of North American naturalized S. purpurea accessions and full-sib F2S. purpurea population were genotyped and phenotyped for a suite of morphological, physiological, pest and disease resistance, and wood chemical composition traits, collected from multi-environment and multi-year replicated field trials. Controlling for population stratification and kinship in the association panel and spatial variation in the F2, a comprehensive mixed model analysis was used to dissect the complex genetic architecture and plasticity of these important traits.

KEY RESULTS

Individually, genome-wide association (GWAS) models differed in terms of power, but the combined approach, which corrects for yearly and environmental co-factors across datasets, improved the overall detection and resolution of associated loci. Although there were few significant GWAS hits located within support intervals of QTL for corresponding traits in the F2, many large-effect QTL were identified, as well as QTL hotspots.

CONCLUSIONS

This study provides the first comparison of linkage analysis and linkage disequilibrium mapping approaches in Salix, and highlights the complementarity and limits of these two methods for elucidating the genetic architecture of complex bioenergy-related traits of a woody perennial breeding programme.

Optimizing the bioenergy water footprint by selecting SRC willow canopy phenotypes: regional scenario simulations

BACKGROUND AND AIMS

Bioenergy is central for the future energy mix to mitigate climate change impacts; however, its intricate link with the water cycle calls for an evaluation of the carbon-water nexus in biomass production. The great challenge is to optimize trade-offs between carbon harvest and water use by choosing cultivars that combine low water use with high productivity.

METHODS

Regional scenarios were simulated over a range of willow genotype × environment interactions for the major UK soil × climate variations with the process-based model LUCASS. Soil available water capacity (SAWC) ranged from 51 to 251 mm and weather represented the north-west (wet, cool), north-east (dry, cool), south-west (wet, warm) and south-east (dry, warm) of the UK. Scenario simulations were evaluated for small/open narrow-leaf (NL) versus large/closed broad-leaf (BL) willow canopy phenotypes using baseline (1965-89) and warmer recent (1990-2014) weather data.

KEY RESULTS

The low productivity under baseline climate in the north could be compensated by choosing BL cultivars (e.g. 'Endurance'). Recent warmer climate increased average productivity by 0.5-2.5 t ha-1, especially in the north. The modern NL cultivar 'Resolution' had the smallest and most efficient water use. On marginal soils (SAWC <100 mm), yields remained below an economic threshold of 9 t ha-1 more frequently under baseline than recent climate. In the drought-prone south-east, 'Endurance' yielded less than 'Resolution', which consumed on average 17 mm year-1 less water. Assuming a planting area of 10 000 ha, in droughty years between 1.3 and 4.5 × 106 m3 of water could be saved, with a small yield penalty, for 'Resolution'.

CONCLUSIONS

With an increase in air temperature and occasional water scarcities expected with climate change, high-yielding NL cultivars should be the preferred choice for sustainable use of marginal lands and reduced competition with agricultural food crops.

Genome-wide transcriptomic analysis of a desert willow, Salix psammophila, reveals the function of hub genes SpMDP1 and SpWRKY33 in drought tolerance

BACKGROUND

Drought is a major environmental constraint to plant growth, development and productivity. Compared with most willows that are generally susceptible to drought, the desert willow Salix psammophila has extraordinary adaptation to drought stress. However, its molecular basis of drought tolerance is still largely unknown.

RESULTS

During polyethylene glycol 6000 (PEG 6000)-simulated drought stress, we found that the osmotic adjustment substances were accumulated and the antioxidant enzyme activities were enhanced in S. psammophila roots. A total of 8172 differentially expressed genes were identified in roots of S. psammophila through RNA-Sequencing. Based on K-means clustering, their expression patterns were classified into nine clusters, which were enriched in several stress-related processes including transcriptional regulation, response to various stresses, cell death, etc. Moreover, 672 transcription factors from 45 gene families were differentially expressed under drought stress. Furthermore, a weighted gene co-expression network was constructed, and eight genes were identified as hub genes. We demonstrated the function of two hub genes, magnesium-dependent phosphatase 1 (SpMDP1) and SpWRKY33, through overexpression in Arabidopsis thaliana. Overexpression of the two hub genes enhanced the drought tolerance in transgenic plants, suggesting that the identification of candidate drought tolerance genes in this study was highly efficient and credible.

CONCLUSIONS

Our study analyzed the physiological and molecular responses to drought stress in S. psammophila, and these results contribute to dissect the mechanism of drought tolerance of S. psammophila and facilitate identification of critical genes involved in drought tolerance for willow breeding.

Mapping QTLs for water-use efficiency reveals the potential candidate genes involved in regulating the trait in apple under drought stress

BACKGROUND

Improvement of water-use efficiency (WUE) can effectively reduce production losses caused by drought stress. A better understanding of the genetic determination of WUE in crops under drought stress has great potential value for developing cultivars adapted to arid regions. To identify the genetic loci associated with WUE and reveal genes responsible for the trait in apple, we aim to map the quantitative trait loci (QTLs) for carbon isotope composition, the proxy for WUE, applying two contrasting irrigating regimes over the two-year experiment and search for the candidate genes encompassed in the mapped QTLs.

RESULTS

We constructed a high-density genetic linkage map with 10,172 markers of apple, using single nucleotide polymorphism (SNP) markers obtained through restriction site-associated DNA sequencing (RADseq) and a final segregating population of 350 seedlings from the cross of Honeycrisp and Qinguan. In total, 33 QTLs were identified for carbon isotope composition in apple under both well-watered and drought-stressed conditions. Three QTLs were stable over 2 years under drought stress on linkage groups LG8, LG15 and LG16, as validated by Kompetitive Allele-Specific PCR (KASP) assays. In those validated QTLs, 258 genes were screened according to their Gene Ontology functional annotations. Among them, 28 genes were identified, which exhibited significant responses to drought stress in 'Honeycrisp' and/or 'Qinguan'. These genes are involved in signaling, photosynthesis, response to stresses, carbohydrate metabolism, protein metabolism and modification, hormone metabolism and transport, transport, respiration, transcriptional regulation, and development regulation. They, especially those for photoprotection and relevant signal transduction, are potential candidate genes connected with WUE regulation in drought-stressed apple.

CONCLUSIONS

We detected three stable QTLs for carbon isotope composition in apple under drought stress over 2 years, and validated them by KASP assay. Twenty-eight candidate genes encompassed in these QTLs were identified. These stable genetic loci and series of genes provided here serve as a foundation for further studies on marker-assisted selection of high WUE and regulatory mechanism of WUE in apple exposed to drought conditions, respectively.

Genetic architecture of aerial and root traits in field-grown grafted grapevines is largely independent

QTLs were identified for traits assessed on field-grown grafted grapevines. Root number and section had the largest phenotypic variance explained. Genetic control of root and aerial traits was independent. Breeding new rootstocks for perennial crops remains challenging, mainly because of the number of desirable traits which have to be combined, these traits include good rooting ability and root development. Consequently, the present study analyzes the genetic architecture of root traits in grapevine. A segregating progeny of 138 F1 genotypes issued from an inter-specific cross between Vitis vinifera cv. Cabernet-Sauvignon × V. riparia cv. Gloire de Montpellier, used as rootstock, was phenotyped in grafted plants grown for 2 years in the field. Seven traits, related to aerial and root development, were quantified. Heritability ranged between 0.44 for aerial biomass to 0.7 for root number. Total root number was related to the number of fine roots, while root biomass was related to the number of coarse roots. Significant quantitative trait loci (QTLs) were identified for all the traits studied with some of them explaining approximately 20% of phenotypic variance. Only a single QTL co-localized for root and aerial biomass. Identified QTLs for aerial-to-root biomass ratio suggest that aerial and root traits are controlled independently. Genes known to be involved in auxin signaling pathways and phosphorus nutrition, whose orthologues were previously shown to regulate root development in Arabidopsis, were located in the confidence intervals of several QTLs. This study opens new perspectives for breeding rootstocks with improved root development capacities.

The genetic architecture of growth traits in Salix matsudana under salt stress

Willow (Salix) is one of the most important ornamental tree species in landscape plants. One species, Salix matsudana, is widely used as a shade tree and border tree because of its soft branches and plump crown. Some varieties of S. matsudana were salt tolerant and could grow normally in coastal regions. However, the molecular mechanisms of salt tolerance for S. matsudana have been less clear. Here, we addressed this issue by performing a mapping experiment containing 195 intraspecific F₁ progeny of S. matsudana, derived from salt-sensitive 'yanjiang' and salt-tolerant '9901', grown by cuttings in a 100 mM NaCl solution. Growth performance of these progeny under salt stress was investigated, displaying marked genotypic variability with the coefficients of variance of 28.64-86.11% in shoot and root growth traits. We further mapped specific QTLs contributing to these differences to the Salix genome. Of the 204 QTLs identified, a few were detected to explain a remarkably larger portion of the phenotypic variation than many others. Many detected QTLs were found to reside in the region of candidate genes of known biological function. The discovery of growth QTLs expressed under salt stress provides important information for marker-assisted breeding of salt tolerant Salix varieties and founds the basis for the application of S. matsudana in coastal afforestation.

Identification of Quantitative Trait Loci Conditioning the Main Biomass Yield Components and Resistance to Melampsora spp. in Salix viminalis × Salix schwerinii Hybrids

The biomass of Salix viminalis is the most highly valued source of green energy, followed by S. schwerinii, S. dasyclados and other species. Significant variability in productivity and leaf rust resistance are noted both within and among willow species, which creates new opportunities for improving willow yield parameters through selection of desirable recombinants supported with molecular markers. The aim of this study was to identify quantitative trait loci (QTLs) linked with biomass yield-related traits and the resistance/susceptibility of Salix mapping population to leaf rust. The experimental material comprised a mapping population developed based on S. viminalis × S. schwerinii hybrids. Phenotyping was performed on plants grown in a field experiment that had a balanced incomplete block design with 10 replications. Based on a genetic map, 11 QTLs were identified for plant height, 9 for shoot diameter, 3 for number of shoots and 11 for resistance/susceptibility to leaf rust. The QTLs identified in our study explained 3%–16% of variability in the analyzed traits. Our findings make significant contributions to the development of willow breeding programs and research into shrubby willow crops grown for energy.

Genome-wide transcriptional and physiological responses to drought stress in leaves and roots of two willow genotypes

BACKGROUND

Drought is a major environmental stress that can have severe impacts on plant productivity and survival. Understanding molecular mechanisms of drought responses is crucial in order to breed for drought adapted plant cultivars. The aim of the present study was to investigate phenotypic and transcriptional drought responses in two willow genotypes (520 and 592) originating from an experimental cross between S. viminalis × (S. viminalis × S. schwerinii). Willows are woody perennials in the Salicaceae plant family that are grown as bioenergy crops worldwide.

METHODS

An experiment was conducted where plants were exposed to drought and different eco-physiological parameters were assessed. RNA-seq data was furthermore generated with the Illumina technology from root tips and leaves from plants grown in drought and well-watered (WW) conditions. The RNA-seq data was assembled de novo with the Trinity assembler to create a reference gene set to which the reads were mapped in order to obtain differentially expressed genes (DEGs) between the drought and WW conditions. To investigate molecular mechanisms involved in the drought response, GO enrichment analyses were conducted. Candidate genes with a putative function in the drought response were also identified.

RESULTS

A total of 52,599 gene models were obtained and after filtering on gene expression (FPKM ≥ 1), 35,733 gene models remained, of which 24,421 contained open reading frames. A total of 5,112 unique DEGs were identified between drought and WW conditions, of which the majority were found in the root tips. Phenotypically, genotype 592 displayed less growth reduction in response to drought compared to genotype 520. At the transcriptional level, genotype 520 displayed a greater response in the leaves as more DEGs were found in genotype 520 compared to genotype 592. In contrast, the transcriptional responses in the root tips were rather similar between the two genotypes. A core set of candidate genes encoding proteins with a putative function in drought response was identified, for example MYBs and bZIPs as well as chlorophyll a/b binding proteins.

DISCUSSION

We found substantial differences in drought responses between the genotypes, both at the phenotypic and transcriptional levels. In addition to the genotypic variation in several traits, we also found indications for genotypic variation in trait plasticity, which could play a role in drought adaptation. Furthermore, the two genotypes displayed overall similar transcriptional responses in the root tips, but more variation in the leaves. It is thus possible that the observed phenotypic differences could be a result of transcriptional differences mostly at the leaf level.

CONCLUSIONS

This study has contributed to a better general understanding of drought responses in woody plants, specifically in willows, and has implications for breeding research towards more drought adapted plants.

Genetic and morphological evidence for introgression between three species of willows

Background

Hybridization and introgression are said to occur relatively frequently in plants, and in particular among different species of willows. However, data on the actual frequency of natural hybridization and introgression is rare. Here, we report the first fine-scale genetic analysis of a contact zone shared between the three basket willow species, Salix dasyclados, S. schwerinii and S. viminalis in the vicinity of the Lake Baikal in Southern Siberia. Individuals were sampled in fourteen populations and classified as pure species or hybrids based on a set of morphological characters. They were then genotyped at 384 nuclear SNP and four chloroplast SSR loci. The STRUCTURE and NewHybrids softwares were used to estimate the frequency and direction of hybridization using genotypic data at the nuclear SNP loci.

Results

As many as 19 % of the genotyped individuals were classified as introgressed individuals and these were mainly encountered in the centre of the contact zone. All introgressed individuals were backcrosses to S. viminalis or S. schwerinii and no F1 or F2 hybrids were found. The rest of the genotyped individuals were classified as pure species and formed two clusters, one with S. schwerinii individuals and the other with S. viminalis and S. dasyclados individuals. The two clusters were significantly genetically differentiated, with F_ST = 0.333 (0.282–0.382, p < 0.001). In contrast, for the chloroplast haplotypes, no genetic differentiation was observed as they were completely shared between the species. Based on morphological classification only 5 % of the individuals were classified as introgressed individuals, which was much less than what was detected using genotypic data.

Conclusions

We have discovered a new willow hybrid zone with relatively high frequency of introgressed individuals. The low frequency of F1 hybrids indicates that ongoing hybridization is limited, which could be because of the presence of reproductive barriers or simply because the conditions are not favorable for hybridization. We further conclude that in order to get a complete picture of the species composition of a hybrid zone it is necessary to use a combination of morphological characters and genetic data from both nuclear and chloroplast markers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0461-7) contains supplementary material, which is available to authorized users.

Multispectral airborne imagery in the field reveals genetic determinisms of morphological and transpiration traits of an apple tree hybrid population in response to water deficit

Genetic studies of response to water deficit in adult trees are limited by low throughput of the usual phenotyping methods in the field. Here, we aimed at overcoming this bottleneck, applying a new methodology using airborne multispectral imagery and in planta measurements to compare a high number of individuals.An apple tree population, grafted on the same rootstock, was submitted to contrasting summer water regimes over two years. Aerial images acquired in visible, near- and thermal-infrared at three dates each year allowed calculation of vegetation and water stress indices. Tree vigour and fruit production were also assessed. Linear mixed models were built accounting for date and year effects on several variables and including the differential response of genotypes between control and drought conditions.Broad-sense heritability of most variables was high and 18 quantitative trait loci (QTLs) independent of the dates were detected on nine linkage groups of the consensus apple genetic map. For vegetation and stress indices, QTLs were related to the means, the intra-crown heterogeneity, and differences induced by water regimes. Most QTLs explained 15-20% of variance.Airborne multispectral imaging proved relevant to acquire simultaneous information on a whole tree population and to decipher genetic determinisms involved in response to water deficit.

Genetic architecture of carbon isotope composition and growth in Eucalyptus across multiple environments

In the context of climate change, the water-use efficiency (WUE) of highly productive tree varieties, such as eucalypts, has become a major issue for breeding programmes. This study set out to dissect the genetic architecture of carbon isotope composition (δ(13) C), a proxy of WUE, across several environments. A family of Eucalyptus urophylla × E. grandis was planted in three trials and phenotyped for δ(13) C and growth traits. High-resolution genetic maps enabled us to target genomic regions underlying δ(13) C quantitative trait loci (QTLs) on the E. grandis genome. Of the 15 QTLs identified for δ(13) C, nine were stable across the environments and three displayed significant QTL-by-environment interaction, suggesting medium to high genetic determinism for this trait. Only one colocalization was found between growth and δ(13) C. Gene ontology (GO) term enrichment analysis suggested candidate genes related to foliar δ(13) C, including two involved in the regulation of stomatal movements. This study provides the first report of the genetic architecture of δ(13) C and its relation to growth in Eucalyptus. The low correlations found between the two traits at phenotypic and genetic levels suggest the possibility of improving the WUE of Eucalyptus varieties without having an impact on breeding for growth.

Genetic strategies for dissecting complex traits in biomass willows (Salix spp.)

Willows are highly diverse catkin-bearing trees and shrubs of the genus Salix. They occur in many growth forms, from tall trees to creeping alpines, and successfully occupy a wide variety of ecological niches. Shrubby willows (sub-genus Vetrix) have many characteristics that render them suited to cultivation in much faster growth cycles than conventional forestry. They respond well to coppicing, can be propagated vegetatively as cuttings and achieve rapid growth with low fertilizer inputs. As a result, willows grown as short rotation coppice are now among the leading commercially grown biomass crops in temperate regions. However, although willows have a long history of cultivation for traditional uses, their industrial use is relatively recent and, compared with major arable crops, they are largely undomesticated. Breeding programmes initiated to improve willow as a biomass crop achieved a doubling of yields within a period of <15 years. These advances were made by selecting for stem characteristics (height and diameter) and coppicing response (shoot number and shoot vigour), as well as resistance to pests, diseases and environmental stress, with little or no knowledge of the genetic basis of these traits. Genetics and genomics, combined with extensive phenotyping, have substantially improved our understanding of the basis of biomass traits in willow for more targeted breeding via marker-assisted selection. Here, we present the strategy we have adopted in which a genetic-based approach was used to dissect complex traits into more defined components for molecular breeding and gene discovery.

Molecular and physiological responses to abiotic stress in forest trees and their relevance to tree improvement

Abiotic stresses, such as drought, salinity and cold, are the major environmental stresses that adversely affect tree growth and, thus, forest productivity, and play a major role in determining the geographic distribution of tree species. Tree responses and tolerance to abiotic stress are complex biological processes that are best analyzed at a systems level using genetic, genomic, metabolomic and phenomic approaches. This will expedite the dissection of stress-sensing and signaling networks to further support efficient genetic improvement programs. Enormous genetic diversity for stress tolerance exists within some forest-tree species, and due to advances in sequencing technologies the molecular genetic basis for this diversity has been rapidly unfolding in recent years. In addition, the use of emerging phenotyping technologies extends the suite of traits that can be measured and will provide us with a better understanding of stress tolerance. The elucidation of abiotic stress-tolerance mechanisms will allow for effective pyramiding of multiple tolerances in a single tree through genetic engineering. Here we review recent progress in the dissection of the molecular basis of abiotic stress tolerance in forest trees, with special emphasis on Populus, Pinus, Picea, Eucalyptus and Quercus spp. We also outline practices that will enable the deployment of trees engineered for abiotic stress tolerance to land owners. Finally, recommendations for future work are discussed.

The genetics of water-use efficiency and its relation to growth in maritime pine

To meet the increasing demand of wood biomass worldwide in the context of climate change, developing improved forest tree varieties for high productivity in water-limited conditions is becoming a major issue. This involves breeding for genotypes combining high growth and moderate water loss and thus high water-use efficiency (WUE). The present work provides original data about the genetics of intrinsic WUE (the ratio between net CO2 assimilation rate and stomatal conductance, also estimated by carbon isotope composition of plant material; δ(13)C) and its relation to growth in Pinus pinaster Ait. First, heritability for δ(13)C was estimated (0.29) using a 15-year-old progeny trial (Landes provenance), with no significant differences among three sites contrasting in water availability. High intersite correlations (0.63-0.91) and significant but low genotype-environment interactions were detected. Secondly, the genetic architectures of δ(13)C and growth were studied in a three-generation inbred pedigree, introducing the genetic background of a more-drought-adapted parent (Corsican provenance), at ages of 2 years (greenhouse) and 9 years (plantation). One of the quantitative trait loci (QTLs) identified in the field experiment, explaining 67% of the phenotypic variance, was also found among the QTLs detected in the greenhouse experiment, where it colocalized with QTLs for intrinsic WUE and stomatal conductance. This work was able to show that higher WUE was not genetically linked to less growth, allowing thus genetic improvement of water use. As far as is known, the heritability and QTL effects estimated here are based on the highest number of genotypes measured to date.

The first genetic linkage map of Eucommia ulmoides

In accordance with pseudo-testcross strategy, the first genetic linkage map of Eucommia ulmoides Oliv. was constructed by an F1 population of 122 plants using amplified fragment length polymorphism (AFLP) markers. A total of 22 AFLP primer combinations generated 363 polymorphic markers. We selected 289 markers segregating as 1:1 and used them for constructing the parent-specific linkage maps. Among the candidate markers, 127 markers were placed on the maternal map LF and 108 markers on the paternal map Q1. The maternal map LF spanned 1116.1 cM in 14 linkage groups with a mean map distance of 8.78 cM; the paternal map Q1 spanned 929.6 cM in 12 linkage groups with an average spacing of 8.61 cM. The estimated coverage of the genome through two methods was 78.5 and 73.9% for LF, and 76.8 and 71.2% for Q1, respectively. This map is the first linkage map of E. ulmoides and provides a basis for mapping quantitative-trait loci and breeding applications.

A molecular genetic linkage map of Eucommia ulmoides and quantitative trait loci (QTL) analysis for growth traits

Eucommia ulmoides is an economically important tree species for both herbal medicine and organic chemical industry. Effort to breed varieties with improved yield and quality is limited by the lack of knowledge on the genetic basis of the traits. A genetic linkage map of E. ulmoides was constructed from a full-sib family using sequence-related amplified polymorphism, amplified fragment length polymorphism, inter-simple sequence repeat and simple sequence repeat markers. In total, 706 markers were mapped in 25 linkage groups covering 2133 cM. The genetic linkage map covered approximately 89% of the estimated E. ulmoides genome with an average of 3.1 cM between adjacent markers. The present genetic linkage map was used to identify quantitative trait loci (QTL) affecting growth-related traits. Eighteen QTLs were found to explain 12.4%–33.3% of the phenotypic variance. This genetic linkage map provides a tool for marker-assisted selection and for studies of genome in E. ulmoides.

Plasticity of primary and secondary growth dynamics in Eucalyptus hybrids: a quantitative genetics and QTL mapping perspective

BACKGROUND

The genetic basis of growth traits has been widely studied in forest trees. Quantitative trait locus (QTL) studies have highlighted the presence of both stable and unstable genomic regions accounting for biomass production with respect to tree age and genetic background, but results remain scarce regarding the interplay between QTLs and the environment. In this study, our main objective was to dissect the genetic architecture of the growth trajectory with emphasis on genotype x environment interaction by measuring primary and secondary growth covering intervals connected with environmental variations.

RESULTS

Three different trials with the same family of Eucalyptus urophylla x E. grandis hybrids (with different genotypes) were planted in the Republic of Congo, corresponding to two QTL mapping experiments and one clonal test. Height and radial growths were monitored at regular intervals from the seedling stage to five years old. The correlation between growth increments and an aridity index revealed that growth before two years old (r = 0.5; 0.69) was more responsive to changes in water availability than late growth (r = 0.39; 0.42) for both height and circumference. We found a regular increase in heritability with time for cumulative growth for both height [0.06 - 0.33] and circumference [0.06 - 0.38]. Heritabilities for incremental growth were more heterogeneous over time even if ranges of variation were similar (height [0-0.31]; circumference [0.19 to 0.48]). Within the trials, QTL analysis revealed collocations between primary and secondary growth QTLs as well as between early growth increments and final growth QTLs. Between trials, few common QTLs were detected highlighting a strong environmental effect on the genetic architecture of growth, validated by significant QTL x E interactions.

CONCLUSION

These results suggest that early growth responses to water availability determine the genetic architecture of total growth at the mature stage and highlight the importance of considering growth as a composite trait (such as yields for annual plants) for a better understanding of its genetic bases.

Rootstock control of scion transpiration and its acclimation to water deficit are controlled by different genes

The stomatal control of transpiration is one of the major strategies by which plants cope with water stress. Here, we investigated the genetic architecture of the rootstock control of scion transpiration-related traits over a period of 3 yr. The rootstocks studied were full sibs from a controlled interspecific cross (Vitis vinifera cv. Cabernet Sauvignon × Vitis riparia cv. Gloire de Montpellier), onto which we grafted a single scion genotype. After 10 d without stress, the water supply was progressively limited over a period of 10 d, and a stable water deficit was then applied for 15 d. Transpiration rate was estimated daily and a mathematical curve was fitted to its response to water deficit intensity. We also determined δ(13) C values in leaves, transpiration efficiency and water extraction capacity. These traits were then analysed in a multienvironment (year and water status) quantitative trait locus (QTL) analysis. Quantitative trait loci, independent of year and water status, were detected for each trait. One genomic region was specifically implicated in the acclimation of scion transpiration induced by the rootstock. The QTLs identified colocalized with genes involved in water deficit responses, such as those relating to ABA and hydraulic regulation. Scion transpiration rate and its acclimation to water deficit are thus controlled genetically by the rootstock, through different genetic architectures.

Optimizing nitrogen economy under drought: increased leaf nitrogen is an acclimation to water stress in willow (Salix spp.)

BACKGROUND AND AIMS

The major objective was to identify plant traits functionally important for optimization of shoot growth and nitrogen (N) economy under drought. Although increased leaf N content (area basis) has been observed in dry environments and theory predicts increased leaf N to be an acclimation to drought, experimental evidence for the prediction is rare.

METHODS

A pedigree of 200 full-sibling hybrid willows was pot-grown in a glasshouse in three replicate blocks and exposed to two water regimes for 3 weeks. Drought conditions were simulated as repeated periods of water shortage. The total leaf mass and area, leaf area efficiency (shoot growth per unit leaf area, E(A)), area-based leaf N content (N(A)), total leaf N pool (N(L)) and leaf N efficiency (shoot growth per unit leaf N, E(N)) were assessed.

KEY RESULTS

In the water-stress treatment, shoot biomass growth was N limited in the genotypes with low N(L), but increasingly limited by other factors in the genotypes with greatest N(L). The N(A) was increased by drought, and drought-induced shift in N(A) varied between genotypes (significant G × E). Judged from the E(A)-N(A) relationship, optimal N(A) was 16 % higher in the water-stress compared with the well-watered treatment. Biomass allocation to leaves and shoots varied between treatments, but the treatment response of the leaf : shoot ratio was similar across all genotypes.

CONCLUSIONS

It is concluded that N-uptake efficiency and leaf N efficiency are important traits to improve growth under drought. Increased leaf N content (area basis) is an acclimation to optimize N economy under drought. The leaf N content is an interesting trait for breeding of willow bioenergy crops in a climate change future. In contrast, leaf biomass allocation is a less interesting breeding target to improve yield under drought.

Polymorphism and Divergence in Two Willow Species, Salix viminalis L. and Salix schwerinii E. Wolf

We investigated species divergence, present and past gene flow, levels of nucleotide polymorphism, and linkage disequilibrium in two willows from the plant genus Salix. Salix belongs together with Populus to the Salicaceae family; however, most population genetic studies of Salicaceae have been performed in Populus, the model genus in forest biology. Here we present a study on two closely related willow species Salix viminalis and S. schwerinii, in which we have resequenced 33 and 32 nuclear gene segments representing parts of 18 nuclear loci in 24 individuals for each species. We used coalescent simulations and estimated the split time to around 600,000 years ago and found that there is currently limited gene flow between the species. Mean intronic nucleotide diversity across gene segments was slightly higher in S. schwerinii (π_i = 0.00849) than in S. viminalis (π_i = 0.00655). Compared with other angiosperm trees, the two willows harbor intermediate levels of silent polymorphisms. The decay of linkage disequilibrium was slower in S. viminalis compared with S. schwerinii, and we speculate that this is due to different demographic histories as S. viminalis has been partly domesticated in Europe.

Genetic improvement of willow for bioenergy and biofuels

Willows (Salix spp.) are a very diverse group of catkin-bearing trees and shrubs that are widely distributed across temperate regions of the globe. Some species respond well to being grown in short rotation coppice (SRC) cycles, which are much shorter than conventional forestry. Coppicing reinvigorates growth and the biomass rapidly accumulated can be used as a source of renewable carbon for bioenergy and biofuels. As SRC willows re-distribute nutrients during the perennial cycle they require only minimal nitrogen fertilizer for growth. This results in fuel chains with potentially high greenhouse gas reductions. To exploit their potential for renewable energy, willows need to be kept free of pests and diseases and yields need to be improved without significantly increasing the requirements for fertilizers and water. The biomass composition needs to be optimized for different end-uses. Yields also need to be sustainable on land less productive for food crops to reduce conflicts over land use. Advances in understanding the physiology and growth of willow, and in the identification of genes underlying key traits, are now at the stage where they can start to be used in breeding programs to help achieve these goals.

High-density linkage mapping and evolution of paralogs and orthologs in Salix and Populus

BACKGROUND

Salix (willow) and Populus (poplar) are members of the Salicaceae family and they share many ecological as well as genetic and genomic characteristics. The interest of using willow for biomass production is growing, which has resulted in increased pressure on breeding of high yielding and resistant clones adapted to different environments. The main purpose of this work was to develop dense genetic linkage maps for mapping of traits related to yield and resistance in willow. We used the Populus trichocarpa genome to extract evenly spaced markers and mapped the orthologous loci in the willow genome. The marker positions in the two genomes were used to study genome evolution since the divergence of the two lineages some 45 mya.

RESULTS

We constructed two linkage maps covering the 19 linkage groups in willow. The most detailed consensus map, S1, contains 495 markers with a total genetic distance of 2477 cM and an average distance of 5.0 cM between the markers. The S3 consensus map contains 221 markers and has a total genetic distance of 1793 cM and an average distance of 8.1 cM between the markers. We found high degree of synteny and gene order conservation between willow and poplar. There is however evidence for two major interchromosomal rearrangements involving poplar LG I and XVI and willow LG Ib, suggesting a fission or a fusion in one of the lineages, as well as five intrachromosomal inversions. The number of silent substitutions were three times lower (median: 0.12) between orthologs than between paralogs (median: 0.37 - 0.41).

CONCLUSIONS

The relatively slow rates of genomic change between willow and poplar mean that the genomic resources in poplar will be most useful in genomic research in willow, such as identifying genes underlying QTLs of important traits. Our data suggest that the whole-genome duplication occurred long before the divergence of the two genera, events which have until now been regarded as contemporary. Estimated silent substitution rates were 1.28 x 10-9 and 1.68 x 10-9 per site and year, which are close to rates found in other perennials but much lower than rates in annuals.

Genetic and environmental variation in spring and autumn phenology of biomass willows (Salix spp.): effects on shoot growth and nitrogen economy

Six commercial willow (Salix spp.) varieties were examined to investigate the effects of genotype and environment on spring and autumn phenology and the relationships between phenology, shoot growth and leaf nitrogen (N) retranslocation. The willows were field-grown under different irrigation and fertilization in central Sweden. Two independent data sets of bud-burst, leaf unfolding duration, growth cessation and the timing of leaf abscission were assessed, and the biomass and leaf N data from the end of the first cutting cycle were used. Specific hypotheses were that (1) spring phenology has a greater effect on the shoot biomass production than autumn phenology; (2) later bud-burst is associated with more rapid leaf unfolding; (3) the timing of leaf abscission has a greater effect on the shoot biomass production than height growth cessation; and (4) later leaf fall is associated with poorer leaf N retranslocation. Bud-burst date varied by 19 and 39 days in the 2 years and leaf unfolding duration varied by 13 and 38 days. Growth cessation varied by 2.5 weeks and completion of leaf abscission (> 90% of leaves shed) by more than 3 weeks between the genotypes and treatments. Bud-burst date was inversely correlated with leaf unfolding duration (R(2) = 0.96). Significant effects of the duration of leafy period (bud-burst to leaf abscission) and bud-burst date on shoot growth were found. Delayed growth cessation and leaf abscission were generally associated with a greater biomass production, but especially the relationship between growth cessation and biomass was weak. The results show that the timing of bud-burst and leaf abscission is more important for willow biomass production than growth cessation. Delayed leaf abscission has a negative effect on leaf N retranslocation and increases the N losses. The results have implications for the breeding of perennial energy crops.

Mapping QTLs for seed yield and drought susceptibility index in soybean (Glycine max L.) across different environments

Drought stress has long been a major constraint in maintaining yield stability of soybean (Glycine max (L.) Merr.) in rainfed ecosystems. The identification of consistent quantitative trait loci (QTL) involving seed yield per plant (YP) and drought susceptibility index (DSI) in a population across different environments would therefore be important in molecular marker-assisted breeding of soybean cultivars suitable for rainfed regions. The YP of a recombinant line population of 184 F(2:7:11) lines from a cross of Kefeng1 and Nannong1138-2 was studied under water-stressed (WS) and well-watered (WW) conditions in field (F) and greenhouse (G) trials, and DSI for yield was calculated in two trials. Nineteen QTLs associated with YP-WS and YP-WW, and 10 QTLs associated with DSI, were identified. Comparison of these QTL locations with previous findings showed that the majority of these regions control one or more traits related to yield and other agronomic traits. One QTL on molecular linkage group (MLG) K for YP-F, and two QTLs on MLG C2 for YP-G, remained constant across different water regimes. The regions on MLG C2 for YP-WW-F and MLG H for YP-WS-F had a pleiotropic effect on DSI-F, and MLG A1 for YP-WS-G had a pleiotropic effect on DSI-G. The identification of consistent QTLs for YP and DSI across different environments will significantly improve the efficiency of selecting for drought tolerance in soybean.

Bioenergy from plants and the sustainable yield challenge

Bioenergy from plants, particularly from perennial grasses and trees, could make a substantial contribution to alleviation of global problems in climate change and energy security if high yields can be sustained. Here, yield traits in a range of key bioenergy crops are reviewed, from which several targets for future improvement can be identified. Some are already the focus of genetically modified (GM) and non-GM approaches. However, the efficient growth strategies of perennial bioenergy crops rely on newly assimilated and recycled carbon and remobilized nitrogen in a continually shifting balance between sources and sinks. This balance is affected by biotic (e.g. pest, disease) and abiotic (e.g. drought) stresses. Future research should focus on three main challenges: changing (photo)thermal time sensitivity to lengthen the growing season without risking frost damage or limiting remobilization of nutritional elements following senescence; increasing aboveground biomass without depleting belowground reserves required for next year's growth and thus without increasing the requirement for nutrient applications; and increasing aboveground biomass without increasing water use.

Genetic basis of phenotypic correlations among growth traits in hybrid willow (Salix dasycladosxS. viminalis) grown under two water regimes

Phenotypic correlations and quantitative trait loci (QTL) for important growth traits and a surrogate of intrinsic water-use efficiency (leaf delta(13)C) were analysed in a willow pedigree of 92 full-sibling clones grown under two water regimes. The major objective was to examine the genetic basis of the phenotypic correlations. Cuttings of Salix were glasshouse-grown during one growing season. The relative growth rate (RGR) and underlying traits were assessed. QTL analysis was conducted based on an available linkage map for Salix. Leaf area productivity and leaf nitrogen productivity were more important in determining RGR than leaf area ratio and specific leaf area. However, phenotypic correlations among growth traits partly varied between the two environments. QTL were detected for most growth traits, among them many common QTL for different traits. The QTL pattern reflected the phenotypic correlation pattern. None of the QTL for the complex traits was consistent across the different environments. The results demonstrate a genetic basis for phenotypic correlations among growth traits in Salix, and provide evidence for the existence of 'master switches' regulating some of the traits.

Forest-tree population genomics and adaptive evolution

Forest trees have gained much attention in recent years as nonclassical model eukaryotes for population, evolutionary and ecological genomic studies. Because of low domestication, large open-pollinated native populations, and high levels of both genetic and phenotypic variation, they are ideal organisms to unveil the molecular basis of population adaptive divergence in nature. Population genomics, in its broad-sense definition, is an emerging discipline that combines genome-wide sampling with traditional population genetic approaches to understanding evolution. Here we briefly review traditional methods of studying adaptive genetic variation in forest trees, and describe a new, integrated population genomics approach. First, alleles (haplotypes) at candidate genes for adaptive traits and their effects on phenotypes need to be characterized via sequencing and association mapping. At this stage, functional genomics can assist in understanding gene action and regulation by providing detailed transcriptional profiles. Second, frequencies of alleles in native populations for causative single-nucleotide polymorphisms are estimated to identify patterns of adaptive variation across heterogeneous environments. Population genomics, through deciphering allelic effects on phenotypes and identifying patterns of adaptive variation at the landscape level, will in the future constitute a useful tool, if cost-effective, to design conservation strategies for forest trees.

Mapping of QTL associated with chilling tolerance during reproductive growth in soybean

Low temperatures in summer bring about drastic reduction in seed yield of soybean [Glycine max (L.) Merr.]. To identify quantitative trait loci (QTL) associated with chilling tolerance during the reproductive growth in soybean, a recombinant inbred line (RIL) population consisting of 104 F(6)-derived lines was created from a cross between two cultivars, chilling-tolerant Hayahikari and chilling-sensitive Toyomusume. The RIL were genotyped with 181 molecular and phenotypic markers and were scored with regard to chilling tolerance, which was evaluated by comparison of seed-yielding abilities in two artificial climatic environments at chilling and usual temperatures. Three QTL were detected for chilling tolerance in seed-yielding ability. Two of them, qCTTSW 1 and qCTTSW 2, were mapped near QTL for flowering time, and the latter had an epistatic interaction with a marker locus located near another QTL for flowering time, where no significant QTL for chilling tolerance was detected. The analysis of an F(2) population derived from the cross between Hayahikari and an RIL of the Hayahikari genotype at all QTL for flowering time confirmed the effect of the third QTL, qCTTSW 3, on chilling tolerance and suggested that qCTTSW 1 was basically independent of the QTL for flowering time. The findings and QTL found in this study may provide useful information for marker-assisted selection (MAS) and further genetic studies on soybean chilling tolerance.