Comparative genomic analysis of the Growth-Regulating Factors-Interacting Factors (GIFs) in six Salicaceae species and functional analysis of PeGIF3 reveals their regulatory role in Populus heteromorphic leaves.

BACKGROUND: The growth-regulating factor-interacting factor (GIF) gene family plays a vital role in regulating plant growth and development, particularly in controlling leaf, seed, and root meristem homeostasis. However, the regulatory mechanism of heteromorphic leaves by GIF genes in Populus euphratica as an important adaptative trait of heteromorphic leaves in response to desert environment remains unknown. RESULTS: This study aimed to identify and characterize the GIF genes in P. euphratica and other five Salicaceae species to investigate their role in regulating heteromorphic leaf development. A total of 27 GIF genes were identified and characterized across six Salicaceae species (P. euphratica, Populus pruinose, Populus deltoides, Populus trichocarpa, Salix sinopurpurea, and Salix suchowensis) at the genome-wide level. Comparative genomic analysis among these species suggested that the expansion of GIFs may be derived from the specific Salicaceae whole-genome duplication event after their divergence from Arabidopsis thaliana. Furthermore, the expression data of PeGIFs in heteromorphic leaves, combined with functional information on GIF genes in Arabidopsis, indicated the role of PeGIFs in regulating the leaf development of P. euphratica, especially PeGIFs containing several cis-acting elements associated with plant growth and development. By heterologous expression of the PeGIF3 gene in wild-type plants (Col-0) and atgif1 mutant of A. thaliana, a significant difference in leaf expansion along the medial-lateral axis, and an increased number of leaf cells, were observed between the overexpressed plants and the wild type. CONCLUSION: PeGIF3 enhances leaf cell proliferation, thereby resulting in the expansion of the central-lateral region of the leaf. The findings not only provide global insights into the evolutionary features of Salicaceae GIFs but also reveal the regulatory mechanism of PeGIF3 in heteromorphic leaves of P. euphratica.

Comprehensive bioinformation analysis of homeodomain-leucine zipper gene family and expression pattern of HD-Zip I under abiotic stress in Salix suchowensis.

BACKGROUND: Homeodomain-leucine zipper (HD-Zip) transcription factors are plant-specific and play important roles in plant defense against environmental stresses. Identification and functional studies have been carried out in model plants such as rice, Arabidopsis thaliana, and poplar, but comprehensive analysis on the HD-Zip family of Salix suchowensis have not been reported. RESULTS: A total of 55 HD-Zip genes were identified in the willow genome, unevenly distributed on 18 chromosomes except for chromosome 19. And segmental duplication events containing SsHD-Zip were detected on all chromosomes except chromosomes 13 and 19. The SsHD-Zip were classified into 4 subfamilies subfamilies (I-IV) according to the evolutionary analysis, and members of each subfamily shared similar domain structure and gene structure. The combination of GO annotation and promoter analysis showed that SsHD-Zip genes responded to multiple abiotic stresses. Furthermore, the results of qPCR analysis showed that the SsHD-Zip I gene exhibited different degrees of expression under salt stress, PEG treatment and heat treatment. Moreover, there was a synergistic effect between SsHD-Zip I genes under stress conditions based on coregulatory networks analysis. CONCLUSIONS: In this study, HD-Zip transcription factors were systematically identified and analyzed at the whole genome level. These results preliminarily clarified the structural characteristics and related functions of willow HD-Zip family members, and it was found that SsHox34, SsHox36 and SsHox51 genes were significantly involved in the response to various stresses. Together, these findings laid the foundation for further research on the resistance functions of willow HD-Zip genes.

Differential salt stress resistance in male and female Salix linearistipularis plants: insights from transcriptome profiling and the identification of the 4-hydroxy-tetrahydrodipicolinate synthase gene.

MAIN CONCLUSION: Lysine plays an essential role in the growth differences between male and female S. linearistipularis plants under salt stress. Furthermore, SlDHDPS is identified as a vital gene contributing to the differences in saline-alkali tolerance between male and female plants of S. linearistipularis. Soil salinization is a significant problem that severely restricts agricultural production worldwide. High salinity and low nutrient concentrations consequently prevent the growth of most plant species. Salix linearistipularis is the only woody plant (shrub) naturally distributed in the saline-alkali lands of the Songnen Plain in Northeast China, and it is one of the few plants capable of thriving in soils with extremely high salt and alkaline pH (>9.0) levels. However, insufficient attention has been given to the interplay between salt and nitrogen in the growth and development of S. linearistipularis. Here, the male and female plants of S. linearistipularis were subjected to salt stress with nitrogen-starvation or nitrogen-supplement treatments, and it was found that nitrogen significantly affects the difference in salt tolerance between male and female plants, with nitrogen-starvation significantly enhancing the salt stress tolerance of female plants compared to male plants. Transcriptional analyses showed 66 differentially expressed nitrogen-responsive genes in female and male roots, with most of them showing sexual differences in expression patterns under salinity stress. RNA-seq and RT-qPCR analysis demonstrated that six genes had an opposite salt-induced expression pattern in female and male roots. The expression of the 4-hydroxy-tetrahydrodipicolinate synthase encoding gene (SlDHDPS) in female roots was higher than that in male roots. Further treatment with exogenous lysine could significantly alleviate the inhibitory effect of salt stress on the growth of female and male plants. These results indicate that the SlDHDPS in the nitrogen metabolism pathway is involved in the resistance of S. linearistipularis to salt stress, which lays a foundation for further exploring the mechanism of nitrogen on salt tolerance of S. linearistipularis, and has a significant reference value for saline-alkali land management and sustainable agricultural development.

Gap-free X and Y chromosome assemblies of Salix arbutifolia reveal an evolutionary change from male to female heterogamety in willows, without a change in the position of the sex-determining locus.

In the Vetrix clade of Salix, a genus of woody flowering plants, sex determination involves chromosome 15, but an XY system has changed to a ZW system. We studied the detailed genetic changes involved. We used genome sequencing, with chromosome conformation capture (Hi-C) and PacBio HiFi reads to assemble chromosome level gap-free X and Y of Salix arbutifolia, and distinguished the haplotypes in the 15X- and 15Y-linked regions, to study the evolutionary history of the sex-linked regions (SLRs). Our sequencing revealed heteromorphism of the X and Y haplotypes of the SLR, with the X-linked region being considerably larger than the corresponding Y region, mainly due to accumulated repetitive sequences and gene duplications. The phylogenies of single-copy orthogroups within the SLRs indicate that S. arbutifolia and Salix purpurea share an ancestral SLR within a repeat-rich region near the chromosome 15 centromere. During the change in heterogamety, the X-linked region changed to a W-linked one, while the Z was derived from the Y.

Climatic variation allows montane willows to escape an adaptive trade-off.

Adaptive responses to climate change, based on heritable variation in stress tolerance, may be important for plant population persistence. It is unclear which populations will mount the strongest future adaptive responses. It may be fruitful to identify populations that have escaped trade-offs among performance traits, which can hinder adaptation. Barring strong genetic constraints, the extent of trade-offs may depend on spatial relationships among climate variables shaping different traits. Here, we test for climate-driven ecotypic variation and trade-offs among drought and freezing sensitivity, and growth, for Lemmon's willow (Salix lemmonii) in a common garden study of 90 genotypes from 38 sites in the Sierra Nevada, USA. Salix lemmonii exhibits ecotypic variation in leaf turgor loss point, a measure of drought sensitivity, from -0.95 to -0.74 MPa along a gradient of spring snowpack. We also find variation in spring freezing sensitivity with minimum May temperature. However, we find no trade-off, as the climatic gradients shaping these traits are spatially uncorrelated in our study region, despite being negatively correlated across the Sierra Nevada. Species may escape adaptive trade-offs in geographic regions where climate variables are spatially decoupled. These regions may represent valuable reservoirs of heritable adaptive phenotypic variation.

Phylogenetic insights into the Salicaceae: The evolution of willows and beyond.

The Salicaceae includes approximately 54 genera and over 1,400 species with a cosmopolitan distribution. Members of the family are well-known for their diverse secondary plant metabolites, and they play crucial roles in tropical and temperate forest ecosystems. Phylogenetic reconstruction of the Salicaceae has been historically challenging due to the limitations of molecular markers and the extensive history of hybridization and polyploidy within the family. Our study employs whole-genome sequencing of 74 species to generate an extensive phylogeny of the Salicaceae. We generated two RAD-Seq enriched whole-genome sequence datasets and extracted two additional gene sets corresponding to the universal Angiosperms353 and Salicaceae-specific targeted-capture arrays. We reconstructed maximum likelihood-based molecular phylogenies using supermatrix and coalescent-based supertree approaches. Our fossil-calibrated phylogeny estimates that the Salicaceae originated around 128 million years ago and unravels the complex taxonomic relationships within the family. Our findings confirm the non-monophyly of the subgenus Salix s.l. and further support the merging of subgenera Chamaetia and Vetrix, both of which exhibit intricate patterns within and among different sections. Overall, our study not only enhances our understanding of the evolution of the Salicaceae, but also provides valuable insights into the complex relationships within the family.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives.

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

Challenge accepted: Evolutionary lineages versus taxonomic classification of North American shrub willows (Salix).

PREMISE: The huge diversity of Salix subgenus Chamaetia/Vetrix clade in North America and the lack of phylogenetic resolution within this clade has presented a difficult but fascinating challenge for taxonomists to resolve. Here we tested the existing taxonomic classification with molecular tools. METHODS: In this study, 132 samples representing 46 species from 22 described sections of shrub willows from the United States and Canada were analyzed and combined with 67 samples from Eurasia. The ploidy levels of the samples were determined using flow cytometry and nQuire. Sequences were produced using a RAD sequencing approach and subsequently analyzed with ipyrad, then used for phylogenetic reconstructions (RAxML, SplitsTree), dating analyses (BEAST, SNAPPER), and character evolution analyses of 14 selected morphological traits (Mesquite). RESULTS: The RAD sequencing approach allowed the production of a well-resolved phylogeny of shrub willows. The resulting tree showed an exclusively North American (NA) clade in sister position to a Eurasian clade, which included some North American endemics. The NA clade began to diversify in the Miocene. Polyploid species appeared in each observed clade. Character evolution analyses revealed that adaptive traits such as habit and adaxial nectaries evolved multiple times independently. CONCLUSIONS: The diversity in shrub willows was shaped by an evolutionary radiation in North America. Most species were monophyletic, but the existing sectional classification could not be supported by molecular data. Nevertheless, monophyletic lineages share several morphological characters, which might be useful in the revision of the taxonomic classification of shrub willows.

Impacts of in-vitro zebularine treatment on genome-wide DNA methylation and transcriptomic profiles in Salix purpurea L.

Renewable energy resources such as biomass are crucial for a sustainable global society. Trees are a major source of lignocellulosic biomass, which can vary in response to different environmental factors owing to epigenetic regulation, such as DNA C-methylation. To investigate the effects of DNA methylation on plant development and wood formation, and its impacts on gene expression, with a focus on secondary cell wall (SCW)-associated genes, Salix purpurea plantlets were cloned from buds derived from a single hybrid tree for both treatment and control conditions. For the treatment condition, buds were exposed to 50 µM zebularine in vitro and a combined strategy of whole-genome bisulfite sequencing (WGBS) and RNA-seq was employed to examine the methylome and transcriptome profiles of different tissues collected at various time points under both conditions. Transcriptomic and methylome data revealed that most of the promoter and gene body demethylation had no marked effects on the expression profiles of genes. Nevertheless, gene expression tended to decrease with the increased methylation levels of genes with highly methylated promoters. Results indicated that demethylation is less evident in centromeric regions and sex chromosomes. Promoters of secondary cell wall-associated genes, such as 4-coumarate-CoA ligase-like and Rac-like GTP-binding protein RHO, were differentially methylated in the secondary xylem samples collected from two-month potted treated plants compared to control samples. Our results provide novel insights into DNA methylation and gene expression landscapes and a basis for investigating the epigenetic regulation of wood formation in S. purpurea as a model plant for bioenergy species.

Curtobacterium salicis sp. nov., isolated from willow tree stems in Washington state.

Curtobacterium sp. strain WW7 is a Gram-positive, non-motile, orange rod-shaped bacterium isolated from branches of wild willow (Salix sitchensis) trees. The WW7T strain has optimum growth in the temperature range between 25 and 30 °C, a pH range of 6-7.7, and tolerates up to 5.5% (w/v) of NaCl. The genome sequencing of strain WW7T revealed a genome size of approximately 3.8 Mbp and a G + C content of 71.3 mol%. The phylogenomic analyses support the WW7T affiliation to a novel Curtobacterium lineage, with Curtobacterium herbarum being the closest type-strain. Chemotaxonomic analysis indicates that the carbon sources assimilation profile of strain WW7T was similar to the type strains, i.e. Curtobacterium luteum, Curtobacterium albidum, and Curtobacterium flaccumfaciens, while no assimilation of the organic acids succinate, alpha-ketobutyrate, mono methyl-succinate, and lactate was observed. Finally, fatty acid methyl ester analysis identifies anteiso-C15:0 and anteiso-C17:0 as major cellular fatty acids which is a common feature for members of the Curtobacterium genus. Based on the results of phylogenomic and chemotaxonomic analyses, strain WW7T represents a novel Curtobacterium lineage, for which the name Curtobacterium salicis sp. nov. is proposed. The type strain is WW7T (DSM 34805T-NRRL B-68078T).

Impact of Heavy Metals on Cold Acclimation of Salix viminalis Roots.

In nature, plants are exposed to a range of climatic conditions. Those negatively impacting plant growth and survival are called abiotic stresses. Although abiotic stresses have been extensively studied separately, little is known about their interactions. Here, we investigate the impact of long-term mild metal exposure on the cold acclimation of Salix viminalis roots using physiological, transcriptomic, and proteomic approaches. We found that, while metal exposure significantly affected plant morphology and physiology, it did not impede cold acclimation. Cold acclimation alone increased glutathione content and glutathione reductase activity. It also resulted in the increase in transcripts and proteins belonging to the heat-shock proteins and related to the energy metabolism. Exposure to metals decreased antioxidant capacity but increased catalase and superoxide dismutase activity. It also resulted in the overexpression of transcripts and proteins related to metal homeostasis, protein folding, and the antioxidant machinery. The simultaneous exposure to both stressors resulted in effects that were not the simple addition of the effects of both stressors taken separately. At the antioxidant level, the response to both stressors was like the response to metals alone. While this should have led to a reduction of frost tolerance, this was not observed. The impact of the simultaneous exposure to metals and cold acclimation on the transcriptome was unique, while at the proteomic level the cold acclimation component seemed to be dominant. Some genes and proteins displayed positive interaction patterns. These genes and proteins were related to the mitigation and reparation of oxidative damage, sugar catabolism, and the production of lignans, trehalose, and raffinose. Interestingly, none of these genes and proteins belonged to the traditional ROS homeostasis system. These results highlight the importance of the under-studied role of lignans and the ROS damage repair and removal system in plants simultaneously exposed to multiple stressors.

Genome-wide identification and expression analysis of xyloglucan endotransglucosylase/hydrolase genes family in Salicaceae during grafting.

BACKGROUND: Poplar (Populus cathayana)and willow (Salix rehderiana) are important fast-growing trees in China. Grafting plays an important role in improving plant stress resistance and construction of ornamental plants. It is found that willow scions grafted onto poplar rootstocks can form ornamental plants. However, this grafted combination has a low survival rate. Many studies have reported that the xyloglucan endotransglucosylase/hydrolase (XTH) family plays an important role in the healing process of grafts. RESULTS: A total of 38 PtrXTHs and 32 SpuXTHs were identified in poplar and willow respectively, and were classified into three subfamilies. Tandem duplication was the main reason for the expansion of the PtrXTHs. Grafting treatment and Quantitative real time PCR (RT-qPCR) analysis revealed that five XTH genes differentially expressed between self-grafted and reciprocal grafted combinations. Specifically, the high expression levels of SrXTH16, SrXTH17, SrXTH25, PcXTH22 and PcXTH17 may contribute to the high survival rate of the grafted combination with willow scion and poplar rootstock. Subcellular localization identified that the SrXTH16, SrXTH17, SrXTH25, PcXTH17 and PcXTH22 proteins were located on the cell walls. Transcription factors (NAC, MYB and DOF) may regulate the five XTH genes. CONCLUSIONS: This study provides a new understanding of the roles of PcXTH and SrXTH genes and their roles in grafting. Our results will give some hints to explore the molecular mechanisms of PcXTH and SrXTH genes involved in grafting in the future.

Contrasting levels of ß-diversity and underlying phylogenetic trends indicate different paths to chemical diversity in highland and lowland willow species.

Diverse specialised metabolites contributed to the success of vascular plants in colonising most terrestrial habitats. Understanding how distinct aspects of chemical diversity arise through heterogeneous environmental pressures can help us understand the effects of abiotic and biotic stress on plant evolution and community assembly. We examined highland and lowland willow species within a phylogenetic framework to test for trends in their chemical α-diversity (richness) and ß-diversity (variation among species sympatric in elevation). We show that differences in chemistry among willows growing at different elevations occur mainly through shifts in chemical ß-diversity and due to convergence or divergence among species sharing their elevation level. We also detect contrasting phylogenetic trends in concentration and α-diversity of metabolites in highland and lowland willow species. The resulting elevational patterns contribute to the chemical diversity of willows and suggest that variable selective pressure across ecological gradients may, more generally, underpin complex changes in plant chemistry.

Structural variation of a sex-linked region confers monoecy and implicates GATA15 as a master regulator of sex in Salix purpurea.

The Salicaceae, including Populus and Salix, are dioecious perennials that utilize different sex determination systems. This family provides a useful system to better understand the evolution of dioecy and sex chromosomes. Here, a rare monoecious genotype of Salix purpurea, 94003, was self- and cross-pollinated and progeny sex ratios were used to test hypotheses on possible mechanisms of sex determination. To delimit genomic regions associated with monoecious expression, the 94003 genome sequence was assembled and DNA- and RNA-Seq of progeny inflorescences was performed. Based on alignments of progeny shotgun DNA sequences to the haplotype-resolved monoecious 94003 genome assembly and reference male and female genomes, a 1.15 Mb sex-linked region on Chr15W was confirmed to be absent in monecious plants. Inheritance of this structural variation is responsible for the loss of a male-suppressing function in what would otherwise be genetic females (ZW), resulting in monoecy (ZWH or WWH ), or lethality, if homozygous (WH WH ). We present a refined, two-gene sex determination model for Salix purpurea, mediated by ARR17 and GATA15 that is different from the single-gene ARR17-mediated system in the related genus Populus.

The male-heterogametic sex determination system on chromosome 15 of Salix triandra and Salix arbutifolia reveals ancestral male heterogamety and subsequent turnover events in the genus Salix.

Dioecious Salix evolved more than 45 million years ago, but have homomorphic sex chromosomes, suggesting that turnover event(s) prevented major differentiation. Sex chromosome turnover events have been inferred in the sister genus Populus. The genus Salix includes two main clades, Salix and Vetrix, with several previously studied Vetrix clade species having female-heterogametic (ZW) or male-heterogametic (XY) sex-determining systems (SDSs) on chromosome 15, while three Salix clade species have XY SDSs on chromosome 7. We here studied two basal taxa of the Vetrix clade, S. arbutifolia and S. triandra using S. purpurea as the reference genome. Analyses of whole genome resequencing data for genome-wide associations (GWAS) with the sexes and genetic differentiation between the sexes (FST values) showed that both species have male heterogamety with a sex-determining locus on chromosome 15, suggesting an early turnover event within the Vetrix clade, perhaps promoted by sexually antagonistic or (and) sex-ratio selection. Changepoint analysis based on FST values identified small sex-linked regions of ~3.33 Mb and ~2.80 Mb in S. arbutifolia and S. triandra, respectively. The SDS of S. arbutifolia was consistent with recent results that used its own genome as reference. Ancestral state reconstruction of SDS suggests that at least two turnover events occurred in Salix.

De Novo Assembly and Annotation of 11 Diverse Shrub Willow (Salix) Genomes Reveals Novel Gene Organization in Sex-Linked Regions.

Poplar and willow species in the Salicaceae are dioecious, yet have been shown to use different sex determination systems located on different chromosomes. Willows in the subgenus Vetrix are interesting for comparative studies of sex determination systems, yet genomic resources for these species are still quite limited. Only a few annotated reference genome assemblies are available, despite many species in use in breeding programs. Here we present de novo assemblies and annotations of 11 shrub willow genomes from six species. Copy number variation of candidate sex determination genes within each genome was characterized and revealed remarkable differences in putative master regulator gene duplication and deletion. We also analyzed copy number and expression of candidate genes involved in floral secondary metabolism, and identified substantial variation across genotypes, which can be used for parental selection in breeding programs. Lastly, we report on a genotype that produces only female descendants and identified gene presence/absence variation in the mitochondrial genome that may be responsible for this unusual inheritance.

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis.

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

Transcriptome-Wide Identification and Response Pattern Analysis of the Salix integra NAC Transcription Factor in Response to Pb Stress.

The NAC (NAM-ATAF1/2-CUC) transcription factor family is one of the largest plant-specific transcription factor families, playing an important role in plant growth and development and abiotic stress response. As a short-rotation woody plant, Salix integra (S. integra) has high lead (Pb) phytoremediation potential. To understand the role of NAC in S. integra Pb tolerance, 53 SiNAC transcripts were identified using third-generation and next-generation transcriptomic data from S. integra exposed to Pb stress, and a phylogenetic analysis revealed 11 subfamilies. A sequence alignment showed that multiple subfamilies represented by TIP and ATAF had a gene that produced more than one transcript under Pb stress, and different transcripts had different responses to Pb. By analyzing the expression profiles of SiNACs at 9 Pb stress time points, 41 of 53 SiNACs were found to be significantly responsive to Pb. Short time-series expression miner (STEM) analysis revealed that 41 SiNACs had two significant Pb positive response patterns (early and late), both containing 10 SiNACs. The SiNACs with the most significant Pb response were mainly from the ATAF and NAP subfamilies. Therefore, 4 and 3 SiNACs from the ATAF and NAP subfamilies, respectively, were selected as candidate Pb-responsive SiNACs for further structural and functional analysis. The RT-qPCR results of 7 transcripts also confirmed the different Pb response patterns of the ATAF and NAP subfamilies. SiNAC004 and SiNAC120, which were randomly selected from two subfamilies, were confirmed to be nuclear localization proteins by subcellular localization experiments. Functional prediction analysis of the associated transcripts of seven candidate SiNACs showed that the target pathways of ATAF subfamily SiNACs were "sulfur metabolism" and "glutathione metabolism", and the target pathways of NAP subfamily SiNACs were "ribosome" and "phenylpropanoid biosynthesis". This study not only identified two NAC subfamilies with different Pb response patterns but also identified Pb-responsive SiNACs that could provide a basis for subsequent gene function verification.

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.

Somatic Mutation Analysis in Salix suchowensis Reveals Early-Segregated Cell Lineages.

Long-lived plants face the challenge of ever-increasing mutational burden across their long lifespan. Early sequestration of meristematic stem cells is supposed to efficiently slow down this process, but direct measurement of somatic mutations that accompanies segregated cell lineages in plants is still rare. Here, we tracked somatic mutations in 33 leaves and 22 adventitious roots from 22 stem-cuttings across eight major branches of a shrub willow (Salix suchowensis). We found that most mutations propagated separately in leaves and roots, providing clear evidence for early segregation of underlying cell lineages. By combining lineage tracking with allele frequency analysis, our results revealed a set of mutations shared by distinct branches, but were exclusively present in leaves and not in roots. These mutations were likely propagated by rapidly dividing somatic cell lineages which survive several iterations of branching, distinct from the slowly dividing axillary stem cell lineages. Leaf is thus contributed by both slowly and rapidly dividing cell lineages, leading to varied fixation chances of propagated mutations. By contrast, each root likely arises from a single founder cell within the adventitious stem cell lineages. Our findings give straightforward evidence that early segregation of meristems slows down mutation accumulation in axillary meristems, implying a plant "germline" paralog to the germline of animals through convergent evolution.