Abstract
The Asian Houbara Bustard (Chlamydotis macqueenii), a partially migratory bird from the western and Central Asian steppes, is listed as vulnerable on the IUCN Red List. This study reassesses the species’ genetic structure using modern genomics to identify evolutionary significant units (ESUs). Following the generation of a de novo reference assembly and resequencing data (114 birds, 10 locations), we integrated genetic results, migratory behaviour, and geography to identify eight hierarchically structured ESUs: four near range edges (Yemen, Mongolia, Eastern Kazakhstan, Israel) and four within the central range (Central-Eastern, Central-Western, North Iran, South Iran). Low genetic diversity and recent inbreeding make ESUs on the range periphery (Israel, Mongolia, Yemen) the most genetically threatened, consistent with the central-marginal hypothesis. ESUs do not cluster according to their migrant/non-migrant status. Geographic distance significantly shaped genetic structure, with longitudinal separation (isolation-by-distance along an east–west axis) emerging as the strongest predictor of differentiation, particularly among high-latitude migrant populations. Our findings underscore the importance of integrating genomic, geographic and behavioural criteria to define intraspecific units that effectively address the conservation needs of widespread species with complex evolutionary dynamics.
Data availability
The sequencing data are accessible in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under BioProject accession number PRJNA1338100 (sample accessions SRR35788338–SRR35788451).
References
Allendorf, F. W., Luikart, G. & Aitken, S. N. Conservation and the Genetics of Populations. 2nd ed. i-xviii, 1–602 (2013).
Ryder, O. A. Species conservation and systematics—The dilemma of subspecies. Trends Ecol. Evol. 1(1), 9–10 (1986).
Waples, R. S. Pacific salmon, Oncorhynchus spp., and the definition of “species” under the Endangered Species Act. Mar. Fish. Rev. 53, 11–22 (1991).
Moritz, C. Defining “Evolutionary Significant Units” for conservation. Trends Ecol. Evol. 9(10), 373–375 (1994).
Crandall, K. A., Bininda-Emonds, O. R. P., Mace, G. M. & Wayne, R. K. Considering evolutionary processes in conservation biology. Trends Ecol. Evol. 15(7), 290–295 (2000).
Fraser, D. J. & Bernatchez, L. Adaptive evolutionary conservation: towards a unified concept for defining conservation units. Mol. Ecol. 10(12), 2741–2752 (2001).
Geue, J. et al. Practical genetic diversity protection: an accessible framework for IUCN subpopulation and Evolutionarily Significant Unit identification (2025).
Hanski, I. Metapopulation dynamics. Nature 396(6706), 41–49 (1998).
Waples, R. S. & Gaggiotti, O. What is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity. Mol. Ecol. 15(6), 1419–1439 (2006).
Eckert, C. G., Samis, K. E. & Lougheed, S. C. Genetic variation across species’ geographical ranges: the central-marginal hypothesis and beyond. Mol. Ecol. 17(5), 1170–1188 (2008).
Alda, F. et al. Genetic diversity, structure and conservation of the endangered Cantabrian Capercaillie in a unique peripheral habitat. Eur. J. Wildl. Res. 59(5), 719–728 (2013).
Mathur, S., Tomecek, J. M., Heniff, A., Luna, R. & DeWoody, J. A. Evidence of genetic erosion in a peripheral population of a North American game bird: the Montezuma quail (Cyrtonyx montezumae). Conserv. Genet. 20(6), 1369–1381 (2019).
Kirchman, J. J., Ross, A. M. & Johnson, G. Historical decline of genetic diversity in a range-periphery population of Spruce Grouse (Falcipennis canadensis) inhabiting the Adirondack Mountains. Conserv. Genet. 21(2), 373–380 (2020).
Davis, D. M., Reese, K. P., Gardner, S. C. & Bird, K. L. Genetic structure of Greater Sage-Grouse (Centrocercus urophasianus) in a declining, peripheral population. Condor 117(4), 530–544 (2015).
Luna, L. W. et al. Historical climate change as driver of populational range expansion and differentiation in a rare and partially migratory Neotropical bird. J. Ornithol. 163(2), 495–507 (2022).
Sexton, J. P., McIntyre, P. J., Angert, A. L. & Rice, K. J. Evolution and ecology of species range limits. Annu. Rev. Ecol. Evol. Syst. 40, 415–436 (2009).
Rehm, E. M., Olivas, P., Stroud, J. & Feeley, K. J. Losing your edge: climate change and the conservation value of range-edge populations. Ecol. Evol. 5(19), 4315–4326 (2015).
IUCN. Chlamydotis macqueenii (2024).
Combreau, O., Riou, S., Judas, J., Lawrence, M. & Launay, F. Migratory pathways and connectivity in Asian Houbara Bustards: Evidence from 15 years of satellite tracking. PLoS ONE 6(6), e20570 (2011).
Combreau, O., Launay, F., Al Bowardi, M. & Gubin, B. Outward migration of Houbara Bustards from two breeding areas in Kazakhstan. Condor 101(1), 159–164 (1999).
Madon, B., Le Nuz, E., Ferlat, C. & Hingrat, Y. Insights into the phenology of migration and survival of a long migrant land bird. bioRxiv (2015).
Azar, J. F., Ferlat, C., Le Nuz, E., Madon, B. & Hingrat, Y. PY -. Variations in the annual migration cycle of bustards: the interplay of geography, age, and sex. [Available from the corresponding author upon request]; in prep.
Madon, B. & Hingrat, Y. Deciphering behavioral changes in animal movement with a “multiple change point algorithm-classification tree” framework. Front. Ecol. Evol. 2, 30 (2014).
Granjon, L., Gaucher, P., Greth, A., Paillat, P. & Vassart, M. Allozyme study of two subspecies of houbara bustard (Chlamydotis undulata undulata and C. u. macqueenii). Biochem. Syst. Ecol. 22(8), 775–779 (1994).
Gaucher,. Taxonomy of the Houbara Bustard Chlamydotis undulata subspecies considered on the basis of sexual display and genetic divergence (vol 138, pg 273, 1996). Ibis 138(3), 570 (1996).
Idaghdour, Y., Broderick, D., Korrida, A. & Chbel, F. Mitochondrial control region diversity of the houbara bustard Chlamydotis undulata complex and genetic structure along the Atlantic seaboard of North Africa. Mol. Ecol. 13(1), 43–54 (2004).
Pitra, C., D’Aloia, M. A., Lieckfeldt, D. & Combreau, O. Genetic variation across the current range of the Asian houbara bustard (Chlamydotis undulata macqueenii). Conserv. Genet. 5(2), 205–215 (2004).
Riou, S. et al. Genetic differentiation among migrant and resident populations of the threatened Asian Houbara Bustard. J. Hered. 103(1), 64–70 (2012).
Gu, Z. R., Dixon, A. & Zhan, X. J. Genetics and evolution of bird migration. Annu. Rev. Anim. Biosci. 12, 21–43 (2024).
Saino, N. et al. Polymorphism at the Clock gene predicts phenology of long-distance migration in birds. Mol. Ecol. 24(8), 1758–1773 (2015).
Berthold, P., Helbig, A. J., Mohr, G. & Querner, U. Rapid microevolution of migratory behavior in a wild bird species. Nature 360(6405), 668–670 (1992).
Zink, R. M. The evolution of avian migration. Biol. J. Lin. Soc. 104(2), 237–250 (2011).
Delmore, K. et al. The evolutionary history and genomics of European blackcap migration. Elife 9, e54462 (2020).
Jombart, T., Devillard, S. & Balloux, F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11 (2010).
Frankham, R. Genetics and extinction. Biol. Conserv. 126(2), 131–140 (2005).
Combreau, O. et al. Breeding success in a Houbara Bustard Chlamydotis [undulata] macqueenii population on the eastern fringe of the Jungar Basin, People’s Republic of China. Ibis 144(2), E45–E56 (2002).
Seddon, P. J. & vanHeezik, Y. Seasonal changes in houbara bustard Chlamydotis undulata macqueenii numbers in Harrat al Harrah, Saudi Arabia: Implications for managing a remnant population. Biol. Conserv. 75(2), 139–146 (1996).
Goriup, P. D. The world status of the houbara bustard Chlamydotis undulata. Bird Conserv. Int. 7, 373–397 (1997).
Tourenq, C. et al. Alarming houbara bustard population trends in Asia. Biol. Conserv. 121(1), 1–8 (2005).
Porter, R. & Eriksen, J. Field Guide to the Birds of Oman (Bloomsbury Publishing, 2017).
Andrews, I. J. The Birds of the Hashemite Kingdom of Jordan. Privately published (1995).
Fitzpatrick, S. W. et al. Genomic and fitness consequences of genetic rescue in wild populations. Curr. Biol. 30(3), 517 (2020).
Gu, Z. R. et al. Climate-driven flyway changes and memory-based long-distance migration. Nature 591(7849), 259 (2021).
Burnside, R. J., Collar, N. J., Scotland, K. M. & Dolman, P. M. Survival rates of captive-bred Asian Houbara Chlamydotis macqueenii in a hunted migratory population. Ibis 158(2), 353–361 (2016).
Hailer, F. et al. Bottlenecked but long-lived high genetic diversity retained in white-tailed eagles upon recovery from population decline. Biol. Lett. 2(2), 316–319 (2006).
Lesobre, L., Azar, J. F., Hoareau, T. B., van Heerden, A. & Hingrat, Y. Genetic determinants of migration in the Asian Houbara Bustard: Insights from radiotracking and genomic techniques. Genet. Migr. (2024).
BirdLife International. Review of the Global Conservation Status of the Asian Houbara Bustard Chlamydotis macqueenii. (BirdLife International, 2014).
de Coster, W., D’Hert, S., Schultz, D. T., Cruts, M. & van Broeckhoven, C. NanoPack: visualizing and processing long-read sequencing data. Bioinformatics 34(15), 2666–2669 (2018).
Kolmogorov, M., Yuan, J., Lin, Y. & Pevzner, P. A. Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 37(5), 540 (2019).
Vaser, R., Sovic, I., Nagarajan, N. & Sikic, M. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 27(5), 737–746 (2017).
Guan, D. F. et al. Identifying and removing haplotypic duplication in primary genome assemblies. Bioinformatics 36(9), 2896–2898 (2020).
Simao, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31(19), 3210–3212 (2015).
gfastats. https://github.com/vgl-hub/gfastats (2022).
Marcais, G. & Kingsford, C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics 27(6), 764–770 (2011).
Vurture, G. W. et al. GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics 33(14), 2202–2204 (2017).
Seddon, P. J., Launay, F., van Heezik, Y. & Al, B. M. Methods for live trapping Houbara Bustards. J. Field Ornithol. 70(2), 169–181 (1999).
FastQC. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (2018).
BBTools bioinformatics tools, including BBMap v.5.85. https://github.com/BioInfoTools/BBMap/releases (2016).
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25(14), 1754–1760 (2009).
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25(16), 2078–2079 (2009).
DePristo, M. A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43(5), 491 (2011).
Poplin, R. et al. Scaling accurate genetic variant discovery to tens of thousands of samples. bioRxiv (2017).
Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27(15), 2156–2158 (2011).
Danecek, P. et al. Twelve years of SAMtools and BCFtools. Gigascience 10(2), giab008 (2021).
Li, H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27(21), 2987–2993 (2011).
Kawakami, T. et al. A high-density linkage map enables a second-generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution. Mol. Ecol. 23(16), 4035–4058 (2014).
Goudet, J. HIERFSTAT, a package for R to compute and test hierarchical F-statistics. Mol. Ecol. Notes 5(1), 184–186 (2005).
Jombart, T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24(11), 1403–1405 (2008).
Thia, J. A. Guidelines for standardizing the application of discriminant analysis of principal components to genotype data. Mol. Ecol. Resour. 23(3), 523–538 (2023).
Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19(9), 1655–1664 (2009).
Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14(8), 2611–2620 (2005).
Korneliussen, T. S., Albrechtsen, A. & Nielsen, R. ANGSD: Analysis of next generation sequencing data. BMC Bioinform. 15, 356 (2014).
Ceballos, F. C., Joshi, P. K., Clark, D. W., Ramsay, M. & Wilson, J. F. Runs of homozygosity: windows into population history and trait architecture. Nat. Rev. Genet. 19(4), 220 (2018).
Purcell, S. et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81(3), 559–575 (2007).
Brüniche-Olsen, A., Kellner, K. F., Anderson, C. J. & DeWoody, J. A. Runs of homozygosity have utility in mammalian conservation and evolutionary studies. Conserv. Genet. 19(6), 1295–1307 (2018).
Meyermans, R., Gorssen, W., Buys, N. & Janssens, S. How to study runs of homozygosity using PLINK? A guide for analyzing medium density SNP data in livestock and pet species. BMC Genomics 21(1), 94 (2020).
Acknowledgements
We are grateful to HH Sheikh Mohamed bin Zayed Al Nahyan, Crown Prince of Abu Dhabi and Founder of the IFHC, HH Sheikh Theyab Bin Mohamed Al Nahyan, Chairman of the IFHC, and HE Mohammed Ahmed Al Bowardi, Deputy Chairman, for their support. This study was conducted under the guidance of Reneco International Wildlife Consultants LLC, a company that manages the IFHC’s conservation programs, such as the National Avian Research Centre. The authors wish to thank the numerous field ecologists from the following organisations for their invaluable efforts in collecting samples over the years: the Emirates Centre for the Conservation of Houbara in Uzbekistan, the Wildlife Science and Conservation Center of Mongolia, the Israel Nature and Heritage Foundation (INHF), and the Israel Nature and Parks Authority (INPA). Special thanks to Joseph Azar for providing ecological insights on the species’ migration, and to Sandra Berthou and Manal Alnaqbi for creating the distribution map.
Funding
Funding was provided by International Funds for Houbara Conservation, Abu Dhabi, UAE.
Author information
Authors and Affiliations
Contributions
L.L. identified the need for this research and secured funding. T.B.H., LL and YH conceptualised and designed the study with inputs from M.M.J. and K.C. Samples were facilitated or provided by E.L.N., N.B., H.O., A.B. and A.M. T.B.H. and K.C. conducted population genetic analyses. T.B.H. interpreted the results and wrote the manuscript, with inputs from L.L., Y.H., M.M.J. and K.C. All authors contributed to the final revisions of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Supplementary Information.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and permissions
About this article
Cite this article
Hoareau, T.B., Collier, K.A., Miller, M.J. et al. Genome wide data recover hierarchical genetic structure and help define conservation units for the threatened Asian Houbara.
Sci Rep (2026). https://doi.org/10.1038/s41598-025-33691-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-025-33691-3
Keywords
- Adaptive evolutionary conservation
- Conservation strategy
- DAPC
- Partially migratory species
- Runs of homozygosity
- Whole genome sequencing
Source: Ecology - nature.com
