Lysine crotonylation, as a novel evolutionarily conserved type of post-translational modifications, is ubiquitous and essential in cell biology. However, its functions in tea plants are largely unknown, and the full functions of lysine crotonylated proteins of tea plants in nitrogen absorption and assimilation remains unclear. Our study attempts to describe the global profiling of nonhistone lysine crotonylation in tea leaves and to explore how ammonium (NH4+) triggers the response mechanism of lysine crotonylome in tea plants. Here, we performed the global analysis of crotonylome in tea leaves under NH4+ deficiency/resupply using high-resolution LC-MS/MS coupled with highly sensitive immune-antibody. A total of 2288 lysine crotonylation sites on 971 proteins were identified, of which contained in 15 types of crotonylated motifs. Most of crotonylated proteins were located in chloroplast (37%) and cytoplasm (33%). Compared with NH4+ deficiency, 120 and 151 crotonylated proteins were significantly changed at 3 h and 3 days of NH4+ resupply, respectively. Bioinformatics analysis showed that differentially expressed crotonylated proteins participated in diverse biological processes such as photosynthesis (PsbO, PsbP, PsbQ, Pbs27, PsaN, PsaF, FNR and ATPase), carbon fixation (rbcs, rbcl, TK, ALDO, PGK and PRK) and amino acid metabolism (SGAT, GGAT2, SHMT4 and GDC), suggesting that lysine crotonylation played important roles in these processes. Moreover, the protein-protein interaction analysis revealed that the interactions of identified crotonylated proteins diversely involved in photosynthesis, carbon fixation and amino acid metabolism. Interestingly, a large number of enzymes were crotonylated, such as Rubisco, TK, SGAT and GGAT, and their activities and crotonylation levels changed significantly by sensing ammonium, indicating a potential function of crotonylation in the regulation of enzyme activities. The results indicated that the crotonylated proteins had a profound influence on metabolic process of tea leaves in response to NH4+ deficiency/resupply, which mainly involved in diverse aspects of primary metabolic processes by sensing NH4+, especially in photosynthesis, carbon fixation and amino acid metabolism. The data might serve as important resources for exploring the roles of lysine crotonylation in N metabolism of tea plants. Data were available via ProteomeXchange with identifier PXD011610.

Despite the high economic and ecological importance of forests, our knowledge of the adaptive evolution of leaf traits remains very limited. Euphrates poplar (Populus euphratica), which has high tolerance to arid environment, has evolved four heteromorphic leaf forms, including narrow (linear and lanceolate) and broad (ovate and broad-ovate) leaves on different crowns. Here, we revealed the significant functional divergence of four P. euphratica heteromorphic leaves at physiological and cytological levels. Through global analysis of transcriptome and DNA methylation across tree and leaf developmental stages, we revealed that gene expression and DNA epigenetics differentially regulated key processes involving development and functional adaptation of heteromorphic leaves, such as hormone signaling pathways, cell division, and photosynthesis. Combined analysis of gene expression, methylation, ATAC-seq, and Hi-C-seq revealed longer interaction of 3D genome, hypomethylation, and open chromatin state upregulates IAA-related genes (such as PIN-FORMED1 and ANGUSTIFOLIA3) and promotes the occurrence of broad leaves while narrow leaves were associated with highly concentrated heterochromatin, hypermethylation, and upregulated abscisic acid pathway genes (such as Pyrabactin Resistance1-like10). Therefore, development of P. euphratica heteromorphic leaves along with functional divergence was regulated by differentially expressed genes, DNA methylation, chromatin accessibility, and 3D genome remodeling to adapt to the arid desert. This study advances our understanding of differential regulation on development and functional divergence of heteromorphic leaves in P. euphratica at the multi-omics level and provides a valuable resource for investigating the adaptive evolution of heteromorphic leaves in Populus.

The Populus pruinosa is a relic plant that has managed to survive in extremely harsh desert environments. Owing to intensifying global warming and desertification, research into ecological adaptation and speciation of P. pruinosa has attracted considerable interest, but the lack of a chromosome-scale genome has limited adaptive evolution research. Here, a 521.09 Mb chromosome-level reference genome of P. pruinosa was reported. Genome evolution and comparative genomic analysis revealed that tandemly duplicated genes and expanded gene families in P. pruinosa contributed to adaptability to extreme desert environments (especially high salinity and drought). The long terminal repeat retrotransposons (LTR-RTs) inserted genes in the gene body region might drive the adaptive evolution of P. pruinosa and species differentiation in saline-alkali desert environments. We recovered genetic differentiation in the populations of the northern Tianshan Mountain and southern Tianshan Mountain through whole-genome resequencing of 156 P. pruinosa individuals from 25 populations in China. Further analyses revealed that precipitation drove the local adaptation of P. pruinosa populations via some genetic sites, such as MAG2-interacting protein 2 (MIP2) and SET domain protein 25 (SDG25). This study will provide broad implications for adaptative evolution and population studies by integrating internal genetic and external environmental factors in P. pruinosa.

NIN-like protein (NLP) transcription factors are key regulators of plant nitrate signaling and stress responses. Although extensively studied in Arabidopsis thaliana and various crops, it has rarely been reported in woody plants, particularly in drought-tolerant tree species. In this study, 10 PeNLP genes were identified in the drought-tolerant tree Populus euphratica Oliv. through comparative genomics. These genes were unevenly distributed across seven chromosomes, and the gene-family expansion was mainly driven by whole-genome duplication (WGD). Analysis of conserved domains showed that PeNLPs contained 4–10 characteristic motifs, and most members possessed the typical RWP-RK and PB1-related domains. Collinearity analysis identified 18 NLP orthologous gene pairs between P. euphratica and its relatives (Populus pruinosa and Salix sinopurpurea), which exceeded the 15 pairs detected between P. euphratica and A. thaliana, indicating that the NLP family is more conserved within the Salicaceae. Phylogenetic analysis divided PeNLPs into three subfamilies, and their promoter regions harbored diverse cis-acting elements associated with hormone signaling, environmental stress, growth, and light response. Transcriptome and qRT-PCR analyses further demonstrated that PeNLPs were generally downregulated under drought stress. Overall, this study systematically characterized the evolution, structure, and drought responsiveness of the PeNLPs, providing a theoretical basis and genetic resources for improving nitrogen use efficiency and drought resistance in trees.

Populus euphratica Oliv. typically has four kinds of heteromorphic leaves: linear (Li), lanceolate (La), ovate (Ov) and broad ovate (Bo). Heteromorphic leaves help P. euphratica adapt to extreme desert environments and further contribute to protection against land desertification in Northwest China. In the authors’ previous research, growth-regulating factors (GRFs) were speculated to be related to the development of P. euphratica heteromorphic leaves via multi-omics integrated analysis. However, the genomic features and biological role of the P. euphratica GRF gene family in heteromorphic leaves are still unclear. In this study, 19 PeGRF genes were genome-widely identified and characterized in P. euphratica, and their physicochemical properties, gene structure and phylogenetic evolution were analyzed. An analysis of the research showed that PeGRFs were unevenly distributed on 11 chromosomes and that PeGRF proteins contained conserved motif 1 (WRC) and motif 2 (QLQ). Moreover, 19, 15, 19 and 22 GRFs were identified in Populus deltoides Marshall, Populus pruinosa Schrenk, Salix sinopurpurea C. Wang et C. Y. Yang and Salix suchowensis W. C. Cheng, respectively. A collinearity analysis showed that the PeGRF family evolved slowly within Populus species. A phylogenetic tree of the GRF family was also constructed, and GRFs were divided into four subfamilies. A large number of cis-acting elements were related to plant growth and development, plant hormone response and stress response on the promoter of PeGRFs. The expression pattern of PeGRFs showed significant up-regulation in broad leaves (Ov and Bo) compared with narrow leaves (Li and La). In combination with the predicted gene regulatory network, PeGRF9 (PeuTF06G01147.1) may have an important contribution to the leaf shape development of P. euphratica. The heterologous expression of PeGRF9 in wild-type plants (Col-0) of Arabidopsis thaliana (L.) Heynh was also studied, showing a significant increase in the leaf area of overexpressing plants compared with the wild type. Nineteen PeGRF gene members were identified and characterized in P. euphratica, and a comparison of the genomic analysis of Populus GRF members revealed their evolutionary features. The further overexpression of PeGRF9 in A. thaliana revealed its biological role in the heteromorphic leaves of P. euphratica. This study not only provides new insights into the evolution and function of PeGRFs in P. euphratica heteromorphic leaves but also helps in an understanding of the adaptive evolution of P. euphratica in drought desert environments.

Abscisic acid (ABA) is a key phytohormone involved in regulating plant growth and responses to environmental stress. As receptors of ABA, pyrabactin resistance 1 (PYR)/PYR1-like (PYL) proteins play a central role in initiating ABA signal transduction. In this study, a total of 30 PopPYL genes were identified and classified into three sub-families (PYL I–III) in the pan-genome of 17 Populus species, through phylogenetic analysis. Among these subfamilies, the PYL I subfamily was the largest, comprising 21 members, whereas PYL III was the smallest, with only four members. To elucidate the evolutionary dynamics of these genes, we conducted synteny and Ka/Ks analyses. Results indicated that most PopPYL genes had undergone purifying selection (Ka/Ks < 1), while a few were subject to positive selection (Ka/Ks > 1). Promoter analysis revealed 258 cis-regulatory elements in the PYL genes of Populus euphratica (EUP) and Populus pruinosa (PRU), including 127 elements responsive to abiotic stress and 33 ABA-related elements. Furthermore, six structural variations (SVs) were detected in PYL_EUP genes and significantly influenced gene expression levels (p < 0.05). To further explore the functional roles of PYL genes, we analyzed tissue-specific expression profiles of 17 PYL_EUP genes under drought stress conditions. PYL6_EUP was predominantly expressed in roots, PYL17_EUP exhibited leaf-specific expression, and PYL1_EUP showed elevated expression in stems. These findings suggest that the drought response of PYL_EUP genes is tissue-specific. Overall, this study highlights the utility of pan-genomics in elucidating gene family evolution and suggests that PYL_EUP genes contribute to the regulation of drought stress responses in EUP, offering valuable genetic resources for functional characterization of PYL genes.