Evidence of Asexual Overwintering of Melampsora paradoxa and Mapping of Stem Rust Host Resistance in Salix.

Melampsora rust is a devastating disease of shrub willow in North America. Previous work has identified Melampsora paradoxa as one of two identified rust species in New York State that infect Salix purpurea and other important Salix host species, however little is known about the population of this rust species in this region. Genotyping-by-sequencing was used to identify single nucleotide polymorphisms (SNPs) and assess population diversity of M. paradoxa isolates collected from three Salix breeding populations in Geneva, NY between 2015 and 2020. Statistical analyses of SNP revealed that all isolates collected were clonally derived even though they were collected across years. In 2020, isolates were collected from stem infections where uredospore pustules were observed, and these isolates were also identical to M. paradoxa collected in previous seasons. These data suggest that M. paradoxa sampled across multiple years overwintered and reproduced asexually and that stem infection is a possible mechanism for overwintering, both of which are novel findings for this rust species. Additionally, field disease ratings were conducted on a S. purpurea × S. suchowensis F1 breeding population with high disease severity, enabling the discovery of QTL for resistance on chromosomes 1 and 19. Lastly, Colletotrichum salicis was frequently associated with stem rust and may play a role in M. paradoxa stem infection. Together, this work is the first substantial exploration into M. paradoxa population biology, stem infection, and host resistance in Salix.

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 F2 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.

The Melampsora americana Population on Salix purpurea in the Great Lakes Region Is Highly Diverse with a Contributory Influence of Clonality.

Shrub willows (Salix spp.) are emerging as a viable lignocellulosic, second-generation bioenergy crop with many growth characteristics favorable for marginal lands in New York State and surrounding areas. Willow rust, caused by members of the genus Melampsora, is the most limiting disease of shrub willow in this region and remains extremely understudied. In this study, genetic diversity, genetic structure, and pathogen clonality were examined in Melampsora americana over two growing seasons via genotyping-by-sequencing to identify single-nucleotide polymorphism markers. In conjunction with this project, a reference genome of rust isolate R15-033-03 was generated to aid in variant discovery. Sampling between years allowed regional and site-specific investigation into population dynamics, in the context of both wild and cultivated hosts within high-density plantings. This work revealed that this pathogen is largely panmictic over the sampled areas, with few sites showing moderate genetic differentiation. These data support the hypothesis of sexual recombination between growing seasons because no genotype persisted across the two years of sampling. Additionally, clonality was determined as a driver of pathogen populations within cultivated fields and single shrubs; however, there is also evidence of high genetic diversity of rust isolates in all settings. This work provides a framework for M. americana population structure in the Great Lakes region, providing crucial information that can aid in future resistance breeding efforts.

Sexual dimorphism in the dioecious willow Salix purpurea.

PREMISE: The evolution of sex chromosomes is driven by sexual dimorphism, yet it can be challenging to document sexually dimorphic traits in dioecious plant species. At the genetic level, sexual dimorphism can be identified through sequence variation between females and males associated with sexually antagonistic traits and different fitness optima. This study aims to examine sexual dimorphism for 26 traits in three populations of Salix purpurea (a diversity panel and F1 and F2 populations) and determine the effect of the traits on biomass yield, a key trait in Salix bioenergy crops across multiple years, locations, and under manipulated growth conditions. METHODS: Sexual dimorphism was evaluated for morphological, phenological, physiological, and wood composition traits in a diversity panel of unrelated S. purpurea accessions and in full-sib F1 and F2 families produced through controlled cross pollinations and grown in replicated field trials. RESULTS: We observed sexual dimorphism in the timing of development for several traits that were highly predictive of biomass yield across three populations of S. purpurea. Across all populations and years surveyed, males had significantly shallower branching angle. Male plants highly predictive of biomass yield across three populations of S. purpurea also accumulated more nitrogen under fertilizer amendment as measured by SPAD in the diversity panel and had greater susceptibility to the rust fungus Melampsora americana in the F2 family. Allometric modelling of biomass yield showed an effect of sex and of location on the interaction between yield and stem height. CONCLUSIONS: These results provide evidence of sexual dimorphism for certain traits in S. purpurea that may be involved in sex chromosome evolution.

Differential Susceptibility of Diverse Salix spp. to Melampsora americana and Melampsora paradoxa.

Melampsora spp. willow rust is the most serious disease of shrub willow bioenergy production in the northeastern United States. Recent phylogenetic studies have identified several Melampsora spp. present on willow in the Northeast; however, in-depth understanding of Melampsora spp. host susceptibility remain unresolved. In this study, a panel of 82 rust isolates collected from the northeastern United States were genotyped via ribosomal DNA sequencing and a subset of these isolates were assayed for host susceptibility. This work revealed that Melampsora americana is the most prevalent species in the sampled geographic region and that there is potential for rust resistance breeding using the Salix spp. taxa assayed. Additionally, leaf morphology traits of these Salix spp. hosts were quantified for correlation analysis, revealing that trichome density and stomata density are possible contributors to resistance. This work provides foundational rust pathology information, which is crucial for M. americana resistance breeding.

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.

Mapping polycomb response elements at the Drosophilla melanogaster giant locus.

Polycomb-group (PcG) proteins are highly conserved epigenetic transcriptional regulators. They are capable of either maintaining the transcriptional silence of target genes through many cell cycles or enabling a dynamic regulation of gene expression in stem cells. In Drosophila melanogaster, recruitment of PcG proteins to targets requires the presence of at least one polycomb response element (PRE). Although the sequence requirements for PREs are not well-defined, the presence of Pho, a PRE-binding PcG protein, is a very good PRE indicator. In this study, we identify two PRE-containing regions at the PcG target gene, giant, one at the promoter, and another approximately 6 kb upstream. PRE-containing fragments, which coincide with localized presence of Pho in chromatin immunoprecipitations, were shown to maintain restricted expression of a lacZ reporter gene in embryos and to cause pairing-sensitive silencing of the mini-white gene in eyes. Our results also reinforce previous observations that although PRE maintenance and pairing-sensitive silencing activities are closely linked, the sequence requirements for these functions are not identical.