Comparative analysis shows high level of lineage sorting in genomic regions with low recombination in the extended Picea likiangensis species complex

•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.