# Genomic Insights into the Origins, Population Structure, and Local Adaptation of Philippine Visayan Native Cattle

**Authors:** Jorge Michael D. Dominguez, Medino Gedeun N. Yebron, Joy B. Banayo, Ningbo Chen, Agapita J. Salces, Kwan Suk Kim

PMC · DOI: 10.3390/ani16040539 · Animals : an Open Access Journal from MDPI · 2026-02-09

## TL;DR

This study explores the genetic origins and adaptations of native cattle in the Visayas region of the Philippines, revealing their indicine ancestry and traits suited to tropical island conditions.

## Contribution

The study identifies genomic regions linked to local adaptation in Visayan native cattle, offering insights into their ancestry and traits for conservation and breeding.

## Key findings

- Visayan native cattle show indicine ancestry bridging mainland Southeast Asia and China with additional taurine and South Asian indicine components.
- Candidate genes suggest adaptations for small stature, heat tolerance, and reproductive performance in tropical island environments.
- Low genetic diversity and long ROH segments in Siquijor cattle indicate recent inbreeding.

## Abstract

Philippine farmers utilize native cattle for meat, draft power, and income security, but until now their origins and perceived traits have been poorly understood. This study utilized genome-wide SNP to analyze native cattle from the Panay and Siquijor islands, Visayas, comparing them to global breeds to uncover their origin and the genomic regions that potentially allow them to thrive in a tropical island environment. Visayan native cattle mostly belonged to indicine cattle and occupied a middle position between mainland Southeast Asian and Southern Chinese cattle, reflecting past mixing among these regions. Their genetic history is consistent with a complex sequence of ancestry signals that may reflect historical livestock movements associated with Southeast Asian maritime connectivity and later colonial-era introductions. Finally, several candidate regions were identified, and these could be related to the traits perceived by local farmers.

The introduction of domestic cattle to the Philippines is often attributed to Spanish and Chinese sources, yet the origins and adaptive history of Philippine Visayan native cattle remain unclear. This study examined the ancestry, structure, and putative selection signals of the Visayan native cattle from Panay and Siquijor islands (VNC) in a global context. Using genome-wide SNP data, population structure was assessed by PCA, IBS/Nei/FST trees, and ADMIXTURE; historical relationships were explored with migration, f-statistics, and an admixture graph; and positive selection was scanned using commonly used methods such as ROH, Tajima’s D, iHS/XP-EHH, and FST with cross-validation across methods and functional enrichment of the overlapping regions. VNC exhibited low-to-moderate genetic diversity (Ho and He ≈ 0.21; and FIS = 0.01 to 0.02) with Siquijor enriched for long ROH segments indicating recent inbreeding. Across multiple complementary analyses, VNC showed predominantly indicine ancestry and occupied an intermediate bridge-like position between indicine from mainland Southeast Asia and from Southeastern China, with additional components that were most similar to Iberian taurine cattle and South Asian indicine. Moreover, the current study identified putative selection signatures that would possibly provide insights to better understand the local adaptation of VNC under insular tropical conditions of the Philippines: (1) small stature (HOXC cluster, STAC3, NXPH4, STARD13, RTN1), (2) heat tolerance and immune robustness (NDUFA4L2, SHMT2, ATP5MC2, ATF7, R3HDM2, CALCOCO1); (3) early reproductive and maturity reproductive performance (IGF2BP2, KL, LRP1, PDS5B). Overall, the VNC in Panay and Siquijor showed a predominantly indicine ancestry with putatively island-adapted physiology, emphasizing the need for conservation and island-specific breeding that preserves local adaptation while managing inbreeding.

## Linked entities

- **Genes:** HOXC@ (homeobox C cluster) [NCBI Gene 3220], STAC3 (SH3 and cysteine rich domain 3) [NCBI Gene 246329], NXPH4 (neurexophilin 4) [NCBI Gene 11247], STARD13 (StAR related lipid transfer domain containing 13) [NCBI Gene 90627], RTN1 (reticulon 1) [NCBI Gene 6252], COXFA4L2 (cytochrome c oxidase hypoxia associated subunit FA4L2) [NCBI Gene 56901], SHMT2 (serine hydroxymethyltransferase 2) [NCBI Gene 6472], ATP5MC2 (ATP synthase membrane subunit c locus 2) [NCBI Gene 517], ATF7 (activating transcription factor 7) [NCBI Gene 11016], R3HDM2 (R3H domain containing 2) [NCBI Gene 22864], CALCOCO1 (calcium binding and coiled-coil domain 1) [NCBI Gene 57658], IGF2BP2 (insulin like growth factor 2 mRNA binding protein 2) [NCBI Gene 10644], KL (klotho) [NCBI Gene 9365], LRP1 (LDL receptor related protein 1) [NCBI Gene 4035], PDS5B (PDS5 cohesin associated factor B) [NCBI Gene 23047]

## Full-text entities

- **Genes:** SHMT2 (serine hydroxymethyltransferase 2) [NCBI Gene 507197], FGF23 (fibroblast growth factor 23) [NCBI Gene 530239], STARD13 (StAR related lipid transfer domain containing 13) [NCBI Gene 538697], NXPH4 (neurexophilin 4) [NCBI Gene 539535], PDS5B (PDS5 cohesin associated factor B) [NCBI Gene 533890], LRP1 (LDL receptor related protein 1) [NCBI Gene 533894], HOXC12 (homeobox C12) [NCBI Gene 613573], MC1R (melanocortin 1 receptor) [NCBI Gene 281298] {aka MC1-R, MSHR}, ATF7 (activating transcription factor 7) [NCBI Gene 541204], BOLA (MHC class I antigen clone 2) [NCBI Gene 751813] {aka BoLA-N}, KL (klotho) [NCBI Gene 784635], STAC3 (SH3 and cysteine rich domain 3) [NCBI Gene 518234], KIT (KIT proto-oncogene, receptor tyrosine kinase) [NCBI Gene 280832] {aka c-kit}, IGF2BP2 (insulin like growth factor 2 mRNA binding protein 2) [NCBI Gene 519028], HOXC@ (homeobox region 3) [NCBI Gene 281822], HOXC6 (homeobox C6) [NCBI Gene 317660], CALCOCO1 (calcium binding and coiled-coil domain 1) [NCBI Gene 538675], HOXC13 (homeobox C13) [NCBI Gene 538716], R3HDM2 (R3H domain containing 2) [NCBI Gene 613499], INS (insulin) [NCBI Gene 280829], ASIP (agouti signaling protein) [NCBI Gene 404192], NDUFA4L2 (NDUFA4 mitochondrial complex associated like 2) [NCBI Gene 613541] {aka COXFA4L2}, RTN1 (reticulon 1) [NCBI Gene 532565]
- **Diseases:** hypoxia (MESH:D000860), gonadal atrophy (MESH:D001284), HD (MESH:D006816), metabolic disease (MESH:D008659), dwarfism (MESH:D004392), small (MESH:D018288), injury to (MESH:D014947)
- **Chemicals:** PAN (MESH:C041728), indicine (MESH:C018336), oxygen (MESH:D010100), lipid (MESH:D008055), silica (MESH:D012822), vitamin D (MESH:D014807), m6A (MESH:C005955), FROH (-), taurine (MESH:D013654), THS (MESH:D013910)
- **Species:** Homo sapiens (human, species) [taxon 9606], Bos gaurus (gaur, species) [taxon 9904], Capra hircus (domestic goat, species) [taxon 9925], Mus musculus (house mouse, species) [taxon 10090], Bos sauveli (kouprey, species) [taxon 135829], Bos indicus (Indicine cattle, species) [taxon 9915], Bos taurus (bovine, species) [taxon 9913], Brachymeles cebuensis (species) [taxon 979815], Bos frontalis (gayal, species) [taxon 30520], Bos javanicus (banteng, species) [taxon 9906], Bubalus mindorensis (tamaraw, species) [taxon 56639]
- **Mutations:** Rs4402960

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12937209/full.md

## Figures

5 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12937209/full.md

## References

173 references — full list in the complete paper: https://tomesphere.com/paper/PMC12937209/full.md

---
Source: https://tomesphere.com/paper/PMC12937209