•Phylogenomic analysis uncovers widespread discordance in the extended Picea likiangensis complex.•Introgression (54.99%) and incomplete lineage sorting (ILS; 33.12%) are key drivers of this incongruity.•Recombination rates shape ILS and introgression, with high rates correlating with elevated levels.•Genes linked to abiotic stress responses exhibit significant introgression and ILS, suggesting adaptive evolution.•Lower recombination rates improve accuracy in species relationships.
© 2024 Kunming Institute of Botany, Chinese Academy of Sciences. Publishing services by Elsevier B.V. on behalf of KeAi Communications Co., Ltd.