† We deeply regret the passing of our esteemed colleague Paul Botriell during the course of the project.

words)Cheetahs (Acinonyx jubatus) are majestic carnivores and the fastest land animals; yet, they are quickly heading towards an uncertain future. Threatened by habitat loss, human-interactions and illegal trafficking, there are only approximately 7,100 individuals remaining in the wild. Cheetahs used to roam large parts of Africa, and Western and Southern Asia. Today they are confined to about 9% of their original distribution. To investigate their genetic diversity and conservation status, we generated genome-wide data from historical and modern samples of all four currently recognized subspecies, along with mitochondrial DNA (mtDNA) and major histocompatibility complex (MHC) data. We found clear genetic differentiation between the subspecies, thus refuting earlier assumptions that cheetahs show only little population differentiation. Our analyses further showed that cheetahs from East Africa are more closely related to A. j. soemmeringii than they are to A . j. juabtus. This raises the need for a reclassification of the cheetah subspecies, as East African cheetahs are currently included into A . j. juabtus. We detected stronger inbreeding in the Critically Endangered A. j. venaticus (Iran) and A. j. hecki (Northwest Africa), and show that overall genome-wide heterozygosity in cheetah is lower than that reported for other threatened and endangered felids, such as tigers and lions. Furthermore, we show that MHC class II diversity in cheetahs is generally higher than previously reported, but still lower than in other felids. Our results provide new and important information for efficient genetic monitoring, subspecies assignments and evidence-based conservation policy decisions.

The cheetah (Acinonyx jubatus, SCHREBER 1775) is a large felid and is considered the fastest land animal. Historically, it inhabited open grassland across Africa, the Arabian Peninsula, and southwestern Asia; however, only small and fragmented populations remain today. Here, we present a de novo genome assembly of the cheetah based on PacBio continuous long reads and Hi-C proximity ligation data. The final assembly (VMU_Ajub_asm_v1.0) has a total length of 2.38 Gb, of which 99.7% are anchored into the expected 19 chromosome-scale scaffolds. The contig and scaffold N50 values of 96.8 Mb and 144.4 Mb, respectively, a BUSCO completeness of 95.4% and a k-mer completeness of 98.4%, emphasize the high quality of the assembly. Furthermore, annotation of the assembly identified 23,622 genes and a repeat content of 40.4%. This new highly contiguous and chromosome-scale assembly will greatly benefit conservation and evolutionary genomic analyses and will be a valuable resource, e.g., to gain a detailed understanding of the function and diversity of immune response genes in felids.

There are only about 7,100 adolescent and adult cheetahs remaining in the wild. With the majority occurring outside protected areas, their numbers are rapidly declining. Evidence-based conservation measures are essential for the survival of this species. Genetic data is routinely used to inform conservation strategies, e.g., by establishing conservation units (CU). A commonly used marker in conservation genetics is mitochondrial DNA (mtDNA). Here, we investigated the cheetah's phylogeography using a large-scale mtDNA data set to re ne subspecies distributions and better assign individuals to CUs. Our dataset mostly consisted of historic samples to cover the cheetah's whole range as the species has been extinct in most of its former distribution. While our genetic data largely agree with geography-based subspecies assignments, several geographic regions show con icting mtDNA signals. Our analyses support previous ndings that evolutionary forces such as incomplete lineage sorting or mitochondrial capture likely confound the mitochondrial phylogeography of this species, especially in East and, to some extent, in Northeast Africa. We caution that subspecies assignments solely based on mtDNA should be treated carefully and argue for an additional standardized nuclear single nucleotide polymorphism (SNP) marker set for subspecies identi cation and monitoring. However, the detection of the A. j. soemmeringii speci c haplogroup by a newly designed Ampli cation-Refractory Mutation System (ARMS) can already provide support for conservation measures.