A Reference Genome for the Critically Endangered Philippine Eagle (Pithecophaga jefferyi), the National Bird of the Philippines

Abstract

The worlds largest and rarest eagle, the Philippine Eagle (Pithecophaga jefferyi), also known as the monkey-eating eagle, is the national bird of the Philippines. This raptor species is endemic to the Philippine archipelago, with populations on the islands of Luzon, Leyte, Samar, and Mindanao. It is critically endangered, with an average estimated population of 392 potentially breeding pairs or 784 mature individuals. In this paper, we describe a reference genome of the Philippine Eagle (Pithecophaga jefferyi) from a female juvenile from the province of Nueva Ecija on the island of Luzon. We generated a de novo genome assembly with high contiguity and completeness, comprising 178 contigs totaling 1.345 Gbp. The genome was sequenced at a coverage of 75.2x, and Benchmarking Universal Single-Copy Orthologs (BUSCO)/Compleasm analysis yielded a BUSCO score of 99.92% (aves_odb12), corresponding to 99.7% single-copy, 0.21% duplicated, and 0.08% fragmented genes. A consensus mitogenome sequence of 19,377 bp was also generated. The genome assembly included 23,847 putative genes, and our annotation estimated that 15.78% of the genome consisted of repetitive elements. Genome heterozygosity (H) was estimated to be 0.020%, in comparison to other birds with genome heterozygosity values ranging from 0.0103% to 0.923%. Whole-genome comparisons with publicly available genomes suggest that the Philippine eagle belongs to the snake-eagle subfamily (Circaetinae) rather than the harpy-eagle subfamily (Harpiinae). Pairwise sequentially Markovian coalescent (PSMC) analysis suggests that the effective population size was around 4,000 individuals from about 100 KYA to about 1 KYA. Finally, we constructed a minimum spanning network, which revealed that our mitogenome from the northern island of Luzon occupies a peripheral position, separated from the dominant haplotype cluster found in the southern island of Mindanao by multiple mutational steps, indicating substantial mitochondrial divergence.