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

Revealing the transitory and local effect of zebularine on development and on proteome dynamics of Salix purpurea

Introduction

DNA methylation plays major roles in the epigenetic regulation of gene expression, transposon and transcriptional silencing, and DNA repair, with implications in developmental processes and phenotypic plasticity. Relevantly for woody species, DNA methylation constitutes a regulative layer in cell wall dynamics associated with xylogenesis. The use of methyltransferase and/or demethylase inhibitors has been proven informative to shed light on the methylome dynamics behind the regulation of these processes.

Methods

The present work employs the cytidine analog zebularine to inhibit DNA methyltransferases and induce DNA hypomethylation in Salix purpurea plantlets grown in vitro and in soil. An integrative approach was adopted to highlight the effects of zebularine on proteomic dynamics, revealing age-specific (3 weeks of in vitro culture and 1 month of growth in soil) and tissue-specific (stem and root) effects.

Results and discussion

After 3 weeks of recovery from zebularine treatment, a decrease of 5-mC levels was observed in different genomic contexts in the roots of explants that were exposed to zebularine, whereas a functionally heterogeneous subset of protein entries was differentially accumulated in stem samples, including entries related to cell wall biosynthesis, tissue morphogenesis, and hormonal regulation. Significant proteomic remodeling was revealed in the development from in vitro to in-soil culture, but no significant changes in 5-mC levels were observed. The identification of tissue-specific proteomic hallmarks in combination with hypomethylating agents provides new insights into the role of DNA methylation and proteome in early plant development in willow species. Proteomic data are available via ProteomeXchange with identifier PXD045653. WGBS data are available under BioProject accession PRJNA889596.

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.

Comparative transcriptomics and eQTL mapping of response to Melampsora americana in selected Salix purpurea F2 progeny

Background

Melampsora spp. rusts are the greatest pathogen threat to shrub willow (Salix spp.) bioenergy crops. Genetic resistance is key to limit the effects of these foliar diseases on host response and biomass yield, however, the genetic basis of host resistance has not been characterized. The addition of new genomic resources for Salix provides greater power to investigate the interaction between S. purpurea and M. americana, species commonly found in the Northeast US. Here, we utilize 3′ RNA-seq to investigate host-pathogen interactions following controlled inoculations of M. americana on resistant and susceptible F₂S. purpurea genotypes identified in a recent QTL mapping study. Differential gene expression, network analysis, and eQTL mapping were used to contrast the response to inoculation and to identify associated candidate genes.

Results

Controlled inoculation in a replicated greenhouse study identified 19 and 105 differentially expressed genes between resistant and susceptible genotypes at 42 and 66 HPI, respectively. Defense response gene networks were activated in both resistant and susceptible genotypes and enriched for many of the same defense response genes, yet the hub genes of these common response modules showed greater mean expression among the resistant plants. Further, eight and six eQTL hotspots were identified at 42 and 66 HPI, respectively. The combined results of three analyses highlight 124 candidate genes in the host for further analysis while analysis of pathogen RNA showed differential expression of 22 genes, two of which are candidate pathogen effectors.

Conclusions

We identified two differentially expressed M. americana transcripts and 124 S. purpurea genes that are good candidates for future studies to confirm their role in conferring resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08254-1.

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.

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.