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.

Genome-Wide Comparative Analysis of the Fasciclin-like Arabinogalactan Proteins (FLAs) in Salicacea and Identification of Secondary Tissue Development-Related Genes.

Fasciclin-like arabinogalactan proteins (FLAs) are a subclass of arabinogalactan proteins (AGPs) containing both AGP-like glycated domains and fasciclin (FAS) domains, which are involved in plant growth and development and synthesis of the cell wall. However, these proteins have not been identified or analyzed in willow, Salix, the sister genus of Populus. In this study, we performed a whole genome study of the FLA gene family of Salix suchowensis and compared it with the FLA gene family of Populus deltoides. The results showed the presence of 40 and 46 FLA genes in P. deltoides and S. suchowensis, distributed on 17 and 16 chromosomes, respectively. Four pairs of tandem repeat genes were found in willow, while poplar had no tandem repeat genes. Twelve and thirteen pairs of duplicated gene fragments were identified in poplar and willow, respectively. The multispecies phylogenetic tree showed that the FLA gene family could be divided into four groups (I-IV), with Group 1 showing significant expansion in woody plants. A gene expression analysis showed that PdeFLA19/27 in Group I of poplar was highly expressed, specifically during the secondary growth period of the stem and the rapid elongation of seed hairs. In the Group I genes of S. suchowensis, SsuFLA25/26/28 was also highly expressed during the secondary growth period, whereas increased expression of SsuFLA35 was associated with seed hair tissue. These results provide important clues about the differences in the FLA gene family during the evolution of herbs and woody plants, and suggest that the FLA gene family may play an essential role in regulating the secondary growth of woody plants. It also provides a reference for further studies on the regulation of secondary growth and seed hair development by FLA genes in poplar and willow.

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.

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 F1 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 F1 hybrid common parent population (Salix F1 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.

Complete chloroplast genome sequencing of five Salix species and its application in the phylogeny and taxonomy of the genus.

In this study, whole chloroplast genomes of five Salix species (S. argyracea, S. dasyclados, S. eriocephala, S. integra 'Hakuro Nishiki', and S. suchowensis) were sequenced. These chloroplast genomes were 155 ,605, 155, 763, 155, 552, 155, 538, and 155 ,550 bp in length, harboring 131 genes (77 unigenes), 37 tRNA genes, 8 rRNA genes, and 86 mRNA genes, respectively. The genes ycf1, psaI, ycf2-2, rpoC2, rpl22, atpF, and ndhF were under positive selection among the 21 Salix species. psaI, ycf2-2, atpF, and ycf1-2 were under positive selection between the tree willow and shrub willow, and rpoC2, rpl22, and ycf1-2 were positively selected among the shrub genomes. The gene rps7 was most variable among the genomes. Phylogenetic analysis of 21 Salix species and Chosenia arbutifolia provide evidence that the cp genome data partially support the relationship with traditional taxonomic concepts in the Flora of China. This chloroplast genome elucidates Salix taxonomy and provides evidence for evolutionary research.

The valid publication of Salix suchowensis (Salicaceae).

The nomenclatural problems of Salix suchowensis have been addressed by different authors with varying opinions. However, these efforts were flawed by a lack of observation of relevant specimens. Accordingly, we carefully checked relevant publications and specimens both through internet databases and herbarium visits. Here, we thoroughly review the nomenclatural history of Salix suchowensis in light of the new definition of a gathering in the Shenzhen Code. We conclude that this name was validly published in the original publication in 1963. Furthermore, a lectotype is designated for the precise application of the name. We hope this article will offer guidance for interpreting similar cases.

Genome-Wide Characterization of the sHsp Gene Family in Salix suchowensis Reveals Its Functions under Different Abiotic Stresses.

Small heat shock proteins (sHsps) function mainly as molecular chaperones that play vital roles in response to diverse stresses, especially high temperature. However, little is known about the molecular characteristics and evolutionary history of the sHsp family in Salix suchowensis, an important bioenergy woody plant. In this study, 35 non-redundant sHsp genes were identified in S. suchowensis, and they were divided into four subfamilies (C, CP, PX, and MT) based on their phylogenetic relationships and predicted subcellular localization. Though the gene structure and conserved motif were relatively conserved, the sequences of the Hsp20 domain were diversified. Eight paralogous pairs were identified in the Ssu-sHsp family, in which five pairs were generated by tandem duplication events. Ka/Ks analysis indicated that Ssu-sHsps had undergone purifying selection. The expression profiles analysis showed Ssu-Hsps tissue-specific expression patterns, and they were induced by at least one abiotic stress. The expression correlation between two paralogous pairs (Ssu-sHsp22.2-CV/23.0-CV and 23.8-MT/25.6-MT) were less than 0.6, indicating that they were divergent during the evolution. Various cis-acting elements related to stress responses, hormone or development, were detected in the promoter of Ssu-sHsps. Furthermore, the co-expression network revealed the potential mechanism of Ssu-sHsps under stress tolerance and development. These results provide a foundation for further functional research on the Ssu-sHsp gene family in S. suchowensis.

Assembly and comparative analysis of complete mitochondrial genome sequence of an economic plant Salix suchowensis.

Willow is a widely used dioecious woody plant of Salicaceae family in China. Due to their high biomass yields, willows are promising sources for bioenergy crops. In this study, we assembled the complete mitochondrial (mt) genome sequence of S. suchowensis with the length of 644,437 bp using Roche-454 GS FLX Titanium sequencing technologies. Base composition of the S. suchowensis mt genome is A (27.43%), T (27.59%), C (22.34%), and G (22.64%), which shows a prevalent GC content with that of other angiosperms. This long circular mt genome encodes 58 unique genes (32 protein-coding genes, 23 tRNA genes and 3 rRNA genes), and 9 of the 32 protein-coding genes contain 17 introns. Through the phylogenetic analysis of 35 species based on 23 protein-coding genes, it is supported that Salix as a sister to Populus. With the detailed phylogenetic information and the identification of phylogenetic position, some ribosomal protein genes and succinate dehydrogenase genes are found usually lost during evolution. As a native shrub willow species, this worthwhile research of S. suchowensis mt genome will provide more desirable information for better understanding the genomic breeding and missing pieces of sex determination evolution in the future.

The whole chloroplast genome of shrub willows (Salix suchowensis).

The whole chloroplast genome of Salix suchowensis was determined in this study using next generation sequencing data. The total genome size was 155,214 bp in length, containing a pair of inverted repeats (IRs) of 27,459 bp, which were separated by large single copy (LSC) and small single copy (SSC) of 84,077 bp and 16,219 bp, respectively. The overall GC contents of the chloroplast genome were 36.73%. One hundred and ten unique genes were annotated, including 76 protein coding genes, 30 tRNA genes and 4 rRNA genes. Among these, 18 are duplicated in the inverted repeat regions, 14 genes contained 1 intron, and 3 genes (rps12, clpP and ycf3) comprised of 2 introns.

The heat shock factor gene family in Salix suchowensis: a genome-wide survey and expression profiling during development and abiotic stresses.

Heat shock transcription factors (Hsfs), which act as important transcriptional regulatory proteins, play crucial roles in plant developmental processes, and stress responses. Recently, the genome of the shrub willow Salix suchowensis was fully sequenced. In this study, a total of 27 non-redundant Hsf genes were identified from the S. suchowensis genome. Phylogenetic analysis revealed that the members of the SsuHsf family can be divided into three groups (class A, B, and C) based on their structural characteristics. Promoter analysis indicated that the SsuHsfs promoters included various cis-acting elements related to hormone and/or stress responses. Furthermore, the expression profiles of 27 SsuHsfs were analyzed in different tissues and under various stresses (heat, drought, salt, and ABA treatment) using RT-PCR. The results demonstrated that the SsuHsfs were involved in abiotic stress responses. Our results contribute to a better understanding of the complexity of the SsuHsf gene family, and will facilitate functional characterization in future studies.