Abstract
The world’s largest butterfly is the microendemic Papua New Guinean Ornithoptera alexandrae. Despite years of conservation efforts to protect its habitat and breed this up-to-28-cm butterfly, this species still figures as endangered in the IUCN Red List and is only known from two allopatric populations occupying a total of only ∼140 km². Here we aim at assembling reference genomes for this species to investigate its genomic diversity, historical demography and determining whether the population is structured, which could provide guidance for conservation programs attempting to (inter)breed the two populations. Using a combination of long and short DNA reads and RNA sequencing, we assembled six reference genomes of the tribe Troidini, with four annotated genomes of O. alexandrae and two genomes of related species O. priamus and Troides oblongomaculatus. We estimated the genomic diversity of the three species, and we proposed scenarios for the historical population demography using two polymorphism-based methods taking into account the characteristics of low-polymorphic invertebrates. Indeed, chromosome-scale assemblies reveal very low levels of nuclear heterozygosity across Troidini, which appears to be exceptionally low for O. alexandrae (lower than 0.01%). Demographic analyses demonstrate low and steadily declining Ne throughout O. alexandrae history, with a divergence into two distinct populations about 10,000 years ago. These results suggest that O. alexandrae distribution has been microendemic for a long time. It should also make local conservation programs aware of the genomic divergence of the two populations, which should not be ignored if any attempt is made to cross the two populations.