Comparing genomic inbreeding of an isolated rhesus macaque study population to wild populations

Long-term studies of isolated animal populations have greatly improved the understanding of various evolutionary processes. However, potentially elevated inbreeding in those compared to wild populations is a common concern. Conventionally, inbreeding has been investigated using reconstructed pedigrees, but nowadays it can be done directly at the genomic level. Here, we utilize genomic data from an intensively studied isolated rhesus macaque (Macaca mulatta) population on the small island Cayo Santiago (Puerto Rico), which was founded in 1938 with wild animals from India. We quantified inbreeding levels by inferring runs of homozygosity (ROH), i.e., long identical haplotypes inherited from both parents. We identified ROH in 97 ~5x-coverage genomes from Cayo Santiago and, for comparison, in 79 rhesus macaque genomes from five wild populations from China. Notably, this conventionally considered low-coverage data proved sufficient to infer ROHs >4 centimorgans long after imputing the genomes using a reference panel. Our results revealed that the ROH-derived effective population size on Cayo Santiago, 420 individuals, falls within the ranges we inferred in wild populations. Moreover, a general scarcity of individuals with long ROH in both the Cayo and wild populations indicates very few cases of close-kin breeding, suggesting that mechanisms to avoid close-kin breeding operate in rhesus macaques, both in wild and isolated populations. Taken together, our results suggest that Cayo Santiago remains a representative study population.