Abstract The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis , P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis , which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica . The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica . Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri , P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.

Forest management changes the physical environments and nutrient dynamics and then regulates the forest productivity. Soil phosphorus (P) availability is critical for productivity in tropical and subtropical forests. However, it was still poorly understood how soil P content and fraction respond to various forest management practices in these regions. Here, we measured the soil total P, available P, and Hedley’s P fractions, including inorganic and organic P (Pi and Po), in subtropical pine plantations treated with understory removal (UR), non-dominant species thinning (NDST) and dominant species thinning (DST) after nine years. Compared to plantations without management (CK), treatments such as UR, NDST, and DST decreased soil total P at 0–10 cm and soil available P at 0–10 cm and 10–20 cm. Increases in resin-Pi, NaOH-Pi, and C.HCl-Pi resulted in a higher total Pi in 0–10 cm ( p < 0.05) in treated plots (UR, NDST, and DST) than in CK plots. UR, NDST, and DST treatments increased NaHCO 3 -Po and NaOH-Po ( p < 0.05) but decreased C.HCl-Po at a depth of 0–10 cm. Regardless of management treatments, soil total P, available P, and P fractions in 0–10 cm showed higher contents than those in 10–20 cm. There were positive relationships between total P and total Po ( p < 0.01) and between available P and total Pi. There were also positive relationships between total P, available P, NaHCO 3 -Pi, and NaOH-Pi ( p < 0.05). In conclusion, forest management such as UR, NDST, and DST decreased soil total P and available P, and transforming soil P fractions to available P will meet the P demand following management in the pine plantations of subtropical China.

To study the interspecific differentiation characteristics of species originating from recent radiation, the genotyping-by-sequencing (GBS) technique was used to explore the kinship, population structure, gene flow, genetic variability, genotype-environment association and selective sweeps of