Leucine-rich repeat receptor-like protein kinases (LRR-RLKs) are the largest group of receptor-like kinases in plants and play crucial roles in development and stress responses. The evolutionary relationships among LRR-RLK genes have been investigated in flowering plants; however, no comprehensive studies have been performed for these genes in more ancestral groups. The subfamily classification of LRR-RLK genes in plants, the evolutionary history and driving force for the evolution of each LRR-RLK subfamily remain to be understood. We identified 119 LRR-RLK genes in the Physcomitrella patens moss genome, 67 LRR-RLK genes in the Selaginella moellendorffii lycophyte genome, and no LRR-RLK genes in five green algae genomes. Furthermore, these LRR-RLK sequences, along with previously reported LRR-RLK sequences from Arabidopsis thaliana and Oryza sativa, were subjected to evolutionary analyses. Phylogenetic analyses revealed that plant LRR-RLKs belong to 19 subfamilies, eighteen of which were established in early land plants, and one of which evolved in flowering plants. More importantly, we found that the basic structures of LRR-RLK genes for most subfamilies are established in early land plants and conserved within subfamilies and across different plant lineages, but divergent among subfamilies. In addition, most members of the same subfamily had common protein motif compositions, whereas members of different subfamilies showed variations in protein motif compositions. The unique gene structure and protein motif compositions of each subfamily differentiate the subfamily classifications and, more importantly, provide evidence for functional divergence among LRR-RLK subfamilies. Maximum likelihood analyses showed that some sites within four subfamilies were under positive selection. Much of the diversity of plant LRR-RLK genes was established in early land plants. Positive selection contributed to the evolution of a few LRR-RLK subfamilies.

Leucine-rich repeat receptor-like protein kinases (LRR-RLKs) are the largest group of receptor-like kinases, which are one of the largest protein superfamilies in plants, and play crucial roles in development and stress responses. Although the evolution of LRR-RLK families has been investigated in some eudicot and monocot plants, no comprehensive evolutionary studies have been performed for these genes in basal angiosperms like Amborella trichopoda. In this study, we identified 94 LRR-RLK genes in the genome of A. trichopoda. The number of LRR-RLK genes in the genome of A. trichopoda is only 17–50% of that of several eudicot and monocot species. Tandem duplication and whole-genome duplication have made limited contributions to the expansion of LRR-RLK genes in A. trichopoda. According to the phylogenetic analysis, all A. trichopoda LRR-RLK genes can be organized into 18 subfamilies, which roughly correspond to the LRR-RLK subfamilies defined in Arabidopsis thaliana. Most LRR-RLK subfamilies are characterized by highly conserved protein structures, motif compositions, and gene structures. The unique gene structure, protein structures, and protein motif compositions of each subfamily provide evidence for functional divergence among LRR-RLK subfamilies. Moreover, the expression data of LRR-RLK genes provided further evidence for the functional diversification of them. In addition, selection analyses showed that most LRR-RLK protein sites are subject to purifying selection. Our results contribute to a better understanding of the evolution of LRR-RLK gene family in angiosperm and provide a framework for further functional investigation on A. trichopoda LRR-RLKs.

Amur honeysuckle (Lonicera maackii) is a widely used medicinal plant of the Caprifoliaceae family that produces chlorogenic acid. Research on this plant mainly focuses on its ornamental value and medicinal compounds, but a reference genome sequence and molecular resources for accelerated breeding are currently lacking. Herein, Nanopore sequencing and high-throughput chromosome conformation capture (Hi-C) allowed a chromosome-level genome assembly of Lonicera maackii (2n = 18). A global view of the gene regulatory network involved in the biosynthesis of chlorogenic acid and the dynamics of fruit coloration in L. maackii was established through metabolite profiling and transcriptome analyses. Moreover, we identified the genes encoding hydroxycinnamoyl-CoA quinate transferase (LmHQT) and hydroxycinnamoyl-CoA shikimic/quinate transferase (LmHCT), which localized to the cytosol and nucleus. Heterologous overexpression of these genes in Nicotiana benthamiana leaves resulted in elevated chlorogenic acid contents. Importantly, HPLC analyses revealed that LmHCT and LmHQTs recombinant proteins modulate the accumulation of chlorogenic acid (CGA) using quinic acid and caffeoyl CoA as substrates, highlighting the importance of LmHQT and LmHCT in chlorogenic acid biosynthesis. These results confirmed that LmHQTs and LmHCT catalyze the biosynthesis of CGA in vitro. The genomic data presented in this study will offer a valuable resource for the elucidation of CGA biosynthesis and facilitating selective molecular breeding.

Chinese sweetgum (Liquidambar formosana) is valued as a source of resin and timber and is an important ornamental tree due to its showy fall foliage. Here, we report the chromosome-level assembly of the Chinese sweetgum genome. Phylogenomic analyses showed the basal phylogenetic position of Chinese sweetgum in core eudicots. Comparative genomic analyses revealed that the well-known gamma event in the common ancestors of core eudicots is evident in the Chinese sweetgum genome, and ancestral triplicated blocks resulting from that event are more intact in Chinese sweetgum than in grapevine (Vitis vinifera). Because of its conserved genome structure, very slow rate of evolution, and basal phylogenetic position, the Chinese sweetgum genome is a good reference for comparative genome studies. Further, we reconstructed the entire metabolic pathway for anthocyanins and potential regulatory networks of autumn leaf coloration of this species via metabolomics and transcriptomics. The transcription factors LfMYB69, basic helix–loop–helix (LfbHLH4), and WD40-repeat (LfWDR1) may collectively regulate the transcription of anthocyanin biosynthetic genes. The regulation of chalcone synthase genes (LfCHS1-3) and dihydroflavonol 4-reductase genes (LfDFR1-2) by the LfMYB69–LfbHLH4–LfWDR1 complex was confirmed by luciferase assays. Epigenomic analyses revealed that 5 structural genes, including LfCHS1, and 2 regulatory LfMYBs are epigenetically regulated. This study expands our understanding of autumn leaf coloration and provides valuable genomic resources for comparative biology, breeding, and biotechnology.