To enhance the treatment capacity and operational stability of the anammox process for nutrient-rich wastewater disposal, a hydroxyapatite (HAP)-anammox coupled process was developed in a granular sludge-based upflow bioreactor. The impact of nitrogen loading rate (NLR) on process performance was investigated by gradually increasing the NLR from 1 to 1.5, 3, and 5 kgN/m3/d for 264 days continuous operation. TN and TP removal rates of 82.51 ± 3.58 % and 60.49 ± 2.82 % were achieved at the NLR of 5 kgN/m3/d. The precipitation of Ca and phosphorus resulted in HAP crystals formation evidenced by sludge morphology and component analysis, which enhanced sludge settling properties and facilitated biomass enrichment. The HAP-anammox granular sludge was characterized by the outer layer HAP crystals formed on inner anammox granules, while highly mineralized granular sludge with higher density distributed at the bottom of the reactor is conducive to HAP-based phosphorus recovery. After a long-term accumulation, the abundance of Planctomycetota phylum gradually increased (from 7.7 % to 14.1 %), consisting of three species of anammox bacteria, namely Candidatus_Brocadia_sapporoensis, Candidatus_Brocadia_sinica and Candidatus_Kuenenia_stuttgartiensis. Nitrogen metabolic pathway analysis illustrated that denitrification was also present except the anammox reaction, synergistically contributing to superior nitrogen removal in the HAP-anammox process.