Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype–environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

Leaf morphology plays a crucial role in predicting the productivity and environmental adaptability of forest trees, making it essential to understand the genetic mechanisms behind leaf variation. In natural populations of Populus cathayana, leaf morphology exhibits rich intraspecific variation due to long-term selection. However, there have been no studies that systematically reveal the genetic mechanisms of leaf variation in P. cathayana. To fill this gap and enhance our understanding of leaf variation in P. cathayana, we collected nine leaf traits from the P. cathayana natural population, consisting of 416 accessions, and conducted the preliminary classification of leaf types with four categories. Subsequently, we conducted an analysis of selective sweep and genome-wide association studies (GWAS) to uncover the genetic basis of leaf traits variation. Most of the leaf traits displayed significant correlations, with broad-sense trait heritability ranging from 0.38 to 0.74. In total, three selective sweep methods ultimately identified 278 positively selected candidate regions and 493 genes associated with leaf size. Single-trait and multi-trait GWAS methods detected 13 and 59 genes, respectively. By integrating the results of selective sweep and GWAS, we further identified a total of nine overlapping genes. These genes may play a role in the leaf development process and are closely associated with leaf size. In particular, the gene CBSCBSPB3 (Pca07G009100) located on chromosome 7, was associated with the response to light stimulation. This study will deepen our understanding of the genetic mechanism of leaf adaptive variation in P. cathayana and provide valuable gene resources.

Light is an important environmental factor for plant growth, and in higher plants, phytochrome A (phyA) is the predominant far-red photoreceptor, involved in various photoresponses. The FAR1/FHY3 transcription factor family, derived from transposases, is able to regulate plant development in response to multiple photosensitizers phytochrome. In total, 51 PtrFRSs were identified in the poplar genome, and were divided into 4 subfamilies. Among them, 47 PtrFRSs are located on 17 chromosomes. Upstream cis-acting elements of the PtrFRS genes were classified into three categories: growth and metabolism, stress and hormone, and the hormone and stress categories contained most of the cis-acting elements. Analysis of the regulatory networks and expression patterns showed that most PtrFRSs responded to changes in light intensity and were involved in the regulation of phytochromes. In this study, 51 PtrFRSs were identified and comprehensively bioinformatically analyzed, and preliminary functional analysis and prediction of PtrFRSs was carried out.

Poplars have a strong response to light, and shade is one of the main environmental factors that limits the growth and development of poplars. Exploring the response mechanism of the developing xylem of poplar under shady conditions is of great reference significance for improving wood yields. In this study, three excellent hybrids of poplar (Populus euramericana ‘Zhonglin46’ (Pe), Populus deltoides ‘27-17’ (Pd), Populus × ‘Wq156’ (Pw) were studied under shady conditions. Based on the phenotypic data and developing a xylem transcriptome analysis, the molecular mechanism of poplars’ response to shade was preliminarily revealed, and the core regulatory genes responding to shade were identified by a weighted co-expression network analysis (WGCNA). The results showed that Pw growth was significantly affected by shade, while Pe growth was slightly affected by shade. An enrichment analysis of 13,675 differentially expressed genes (DEGs) found that shade affected the expression of genes related to the glutathione metabolic pathway. The WGCNA analysis identified two modules (“Brown” and “Purple”) related to the shade response and discovered seven hub genes. These hub genes were related to xylem development, vascular cambium division, stomatal development, and phytochrome A signal transduction. These results provide important basic information for gaining insight into the molecular response to shade in different poplar hybrids.