Elm (Ulmus) has a long history of use as a high-quality heavy hardwood famous for its resistance to drought, cold, and salt. It grows in temperate, warm temperate, and subtropical regions. This is the first report of Ulmaceae chloroplast genomes by de novo sequencing. The Ulmus chloroplast genomes exhibited a typical quadripartite structure with two single-copy regions (long single copy [LSC] and short single copy [SSC] sections) separated by a pair of inverted repeats (IRs). The lengths of the chloroplast genomes from five Ulmus ranged from 158,953 to 159,453 bp, with the largest observed in Ulmus davidiana and the smallest in Ulmus laciniata. The genomes contained 137–145 protein-coding genes, of which Ulmus davidiana var. japonica and U. davidiana had the most and U. pumila had the fewest. The five Ulmus species exhibited different evolutionary routes, as some genes had been lost. In total, 18 genes contained introns, 13 of which (trnL-TAA⁺, trnL-TAA⁻, rpoC1^-, rpl2^-, ndhA^-, ycf1, rps12^-, rps12⁺, trnA-TGC⁺, trnA-TGC^-, trnV-TAC^-, trnI-GAT⁺, and trnI-GAT) were shared among all five species. The intron of ycf1 was the longest (5,675bp) while that of trnF-AAA was the smallest (53bp). All Ulmus species except U. davidiana exhibited the same degree of amplification in the IR region. To determine the phylogenetic positions of the Ulmus species, we performed phylogenetic analyses using common protein-coding genes in chloroplast sequences of 42 other species published in NCBI. The cluster results showed the closest plants to Ulmaceae were Moraceae and Cannabaceae, followed by Rosaceae. Ulmaceae and Moraceae both belonged to Urticales, and the chloroplast genome clustering results were consistent with their traditional taxonomy. The results strongly supported the position of Ulmaceae as a member of the order Urticales. In addition, we found a potential error in the traditional taxonomies of U. davidiana and U. davidiana var. japonica, which should be confirmed with a further analysis of their nuclear genomes. This study is the first report on Ulmus chloroplast genomes, which has significance for understanding photosynthesis, evolution, and chloroplast transgenic engineering.

The WRKY transcription factor (TF) family, named for its iconic WRKY domain, is among the largest and most functionally diverse TF families in higher plants. WRKY TFs typically interact with the W-box of the target gene promoter to activate or inhibit the expression of downstream genes; these TFs are involved in the regulation of various physiological responses. Analyses of WRKY TFs in numerous woody plant species have revealed that WRKY family members are broadly involved in plant growth and development, as well as responses to biotic and abiotic stresses. Here, we review the origin, distribution, structure, and classification of WRKY TFs, along with their mechanisms of action, the regulatory networks in which they are involved, and their biological functions in woody plants. We consider methods currently used to investigate WRKY TFs in woody plants, discuss outstanding problems, and propose several new research directions. Our objective is to understand the current progress in this field and provide new perspectives to accelerate the pace of research that enable greater exploration of the biological functions of WRKY TFs.

In this study, we assembled and annotated the chloroplast (cp) genomes of four Ligustrum species, L. sinense, L. obtusifolium, L. vicaryi, and L. ovalifolium ‘Aureum’. Including six other published Ligustrum species, we compared various characteristics such as gene structure, sequence alignment, codon preference, and nucleic acid diversity, and performed positive-selection genes screening and phylogenetic analysis. The results showed that the cp genome of Ligustrum was 162,185–166,800 bp in length, with a circular tetrad structure, including a large single-copy region (86,885–90,106 bp), a small single-copy region (11,446–11,499 bp), and a pair of IRa and IRb sequences with the same coding but in opposite directions (31,608–32,624 bp). This structure is similar to the cp genomes of most angiosperms. We found 132–137 genes in the cp genome of Ligustrum, including 89–90 protein-coding genes, 35–39 tRNAs, and 8 rRNAs. The GC content was 37.93–38.06% and varied among regions, with the IR region having the highest content. The single-nucleotide (A/T)n was dominant in simple-sequence repeats of the Ligustrum cp genome, with an obvious A/T preference. Six hotspot regions were identified from multiple sequence alignment of Ligustrum; the ycf1 gene region and the clpP1 exon region can be used as potential DNA barcodes for the identification and phylogeny of the genus Ligustrum. Branch-site model and Bayes empirical Bayes (BEB) analysis showed that four protein-coding genes (accD, clpP, ycf1, and ycf2) were positively selected, and BEB analysis showed that accD and rpl20 had positively selected sites. A phylogenetic tree of Oleaceae species was constructed based on the whole cp genomes, and the results were consistent with the traditional taxonomic results. The phylogenetic results showed that genus Ligustrum is most closely related to genus Syringa. Our study provides important genetic information to support further investigations of the phylogenetic development and adaptive evolution of Ligustrum species.