Evolutionary stasis of the pseudoautosomal boundary in strepsirrhine primates

  1. Rylan Shearn
  2. Alison E Wright
  3. Sylvain Mousset
  4. Corinne Régis
  5. Simon Penel
  6. Jean-François Lemaitre
  7. Guillaume Douay
  8. Brigitte Crouau-Roy
  9. Emilie Lecompte
  10. Gabriel AB Marais  Is a corresponding author
  1. Laboratoire Biométrie et Biologie Evolutive, CNRS / Univ. Lyon 1, France
  2. Department of Animal and Plant Sciences, University of Sheffield, United Kingdom
  3. Faculty of Mathematics, University of Vienna, Austria
  4. Zoo de Lyon, France
  5. Laboratoire Evolution et Diversité Biologique, CNRS / Univ. Toulouse, France
  6. LEAF-Linking Landscape, Environment, Agriculture and Food Dept, Instituto Superior de Agronomia, Universidade de Lisboa, Portugal


Sex chromosomes are typically comprised of a non-recombining region and a recombining pseudoautosomal region. Accurately quantifying the relative size of these regions is critical for sex-chromosome biology both from a functional and evolutionary perspective. The evolution of the pseudoautosomal boundary (PAB) is well documented in haplorrhines (apes and monkeys) but not in strepsirrhines (lemurs and lorises). Here, we studied the PAB of seven species representing the main strepsirrhine lineages by sequencing a male and a female genome in each species and using sex differences in coverage to identify the PAB. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events moved the PAB and shortened the pseudoautosomal region in haplorrhines. Strepsirrhines are well known to have much lower sexual dimorphism than haplorrhines. We suggest that mutations with antagonistic effects between males and females have driven recombination suppression and PAB evolution in haplorrhines


The human sex chromosomes are strongly heteromorphic as they exhibit extensive differences in size, gene number, DNA repeat abundance and heterochromatin composition (Skaletsky et al., 2003; Ross et al., 2005). The X chromosome comprises a large X-specific region recombining only in females whereas the Y comprises a male-specific region that does not recombine at all. Both sex chromosomes share two pseudoautosomal regions (PAR1 and 2) that recombine in both males and females. These sex chromosomes originated from a pair of identical autosomes approximately 150 million years ago, prior to the divergence of placentals and marsupials, with the evolution of Sry – the master male-determining gene in therian mammals – from Sox3 (Lahn and Page, 1999; Skaletsky et al., 2003; Hughes and Rozen, 2012). Since then, at several moments throughout evolutionary history, vast regions of the Y chromosome have stopped recombining with the X, likely through inversions on the Y (Lahn and Page, 1999; Van Laere et al., 2008; Lemaitre et al., 2009; Pandey et al., 2013). These regions show different levels of X-Y divergence and are called evolutionary strata (Lahn and Page, 1999). Strata 1 and 2 are shared among all therians, and stratum three is shared among all placentals (Lahn and Page, 1999; Cortez et al., 2014). The most recent strata (4 and 5) have originated in the history of Catarrhini (Old World monkeys and apes) respectively,~40 and~25 Mya, and now only a very small PAR continues to recombine between X and Y in those primates (Hughes et al., 2012). In humans, PAR1 is the consequence of that process, while PAR2 is a recent addition (Skaletsky et al., 2003).

The process of recombination suppression between sex chromosomes, leading to a reduction in the size of the PAR and formation of evolutionary strata, has been documented in several animal and plant lineages (e.g Nicolas et al., 2005; Zhou et al., 2014; White et al., 2015). Why such a process occurred, however, is unclear. It has been proposed that sexually antagonistic mutations may have favoured the suppression of recombination (Bull, 1983; Rice, 1987; Charlesworth et al., 2005). Theoretical models suggest that if there are male-beneficial/female-detrimental mutations in the PAR, there will be selection to halt recombination, through for example an inversion, to genetically link those mutations to the Y chromosome. Some evidence supporting this hypothesis has recently been found in guppies (Wright et al., 2017), but evidence from a wide range of groups, including primates, is lacking. Furthermore, there are alternative theories for why recombination is halted (reviewed in Charlesworth, 2017; Ponnikas et al., 2018) and so the relative importance of sexual antagonism in sex-chromosome evolution remains unclear.

While previous work on primate sex chromosomes has focused on Haplorrhini (apes, Old and New World monkeys), we studied representatives of the other main primate lineage, the Strepsirrhini (lemurs and lorises). In strepsirrhines, female social dominance (FSD), in which females dominate males, is widespread and likely ancestral (Kappeler and Fichtel, 2015; Petty and Drea, 2015). FSD is associated with increased testosterone production in females, resulting in the masculinization of females, including aspects of their social behaviour and genitalia (Kappeler and Fichtel, 2015; Petty and Drea, 2015). Some species also have rather egalitarian social systems (Pereira and Kappeler, 1997). In addition, sexual size dimorphism is virtually absent among strepsirrhines (Kappeler and Fichtel, 2015; Petty and Drea, 2015). This is in sharp contrast with haplorrhines, where sexual dimorphism is much more pronounced and male-biased; a phenotype that is probably ancient in this group (e.g. Lindenfors, 2002; Kappeler and van Schaik, 2004; Plavcan, 2004). We therefore hypothesized that if male–female differentiation and sexually antagonistic mutations were associated with the degree of X-Y recombination suppression, strepsirrhines should show evidence of less recombination suppression compared to haplorrhines. However, to date, very little is known about the sex chromosomes of strepsirrhines, except that strata 4 and 5 are missing in grey mouse lemurs (Microcebus murinus, see Glaser et al., 1999) preventing previous tests of this hypothesis.

To identify the PAB of strepsirrhines, we used an approach relying on sequencing a male and a female at low/moderate depth, mapping the reads to a reference genome and computing the male:female depth ratio (Vicoso and Bachtrog, 2011; Vicoso et al., 2013a; Vicoso et al., 2013b; Zhou et al., 2014). For autosomes, a M:F depth ratio of 1 is expected as males and females have the same copy number of autosomes. On the X chromosome, a ratio of 1 should indicate the PAR that is shared among sexes, a ratio of 0.5 should indicate the X-specific region as males have only one such region and females two, and the boundary between both would indicate the PAB. Using Illumina short-read sequencing technology, we sequenced a male and a female genome in seven species covering the main strepsirrhine lineages representing 65 My of evolution (Pozzi et al., 2014): four Lemuriformes (Daubentonia madagascariensis - aye-ayes, M. murinus – grey mouse lemur, Eulemur rubriventer – red-bellied lemur, Prolemur simus – greater bamboo lemur) and three Lorisiformes (Otolemur garnettii – northern greater galago, Galago senegalensis – senegal bushbaby, Nyctibebus coucang – slow loris). The sequencing depth of each sample was between 11.8X and 39.1X (assuming a genome size identical to the human genome) with 78% of the samples being between 20X and 40X, that is moderate sequencing depth (Supplementary file 1A). We then mapped the reads onto publicly available reference genomes of two strepsirrhines (using the human X to scaffold the strepsirrhine X chromosomes) and computed a normalised M:F depth ratio to identify the X-specific region and the PAR on the X chromosome (see Materials and methods).


Figure 1A–B shows the results for the grey mouse lemur. Using the human X chromosome to order the grey mouse lemur X scaffolds, we found that the scaffolds corresponding to human PAR1 and strata 4 and 5 have a M:F depth ratio around 1 (Figure 1B), indicating that these regions have remained pseudoautosomal in grey mouse lemurs in agreement with older cytogenetic data (Glaser et al., 1999). The rest of the grey mouse lemur X is X-specific with a M:F ratio close to 0.5. However, five regions in the X-specific region show an elevated ratio. Detailed analysis of these five regions showed that they are fragments of autosomes (see Materials and methods and Supplementary file 1B). It is not clear, however, whether this comes from contamination of the assembly of the X chromosome by autosomal scaffolds or if this has resulted from fusion of autosomal DNA fragments to the PAR during evolution, which are misplaced in the current assembly of the X chromosome. With the fragmented assembly that is available our approach can only reliably identify the PAB, not the size of the PAR. If some autosomal material were translocated to the PAR, and thus enlarging it, it would not be possible to detect it with our approach. Only an improved assembly of the X chromosome in the grey mouse lemur could confirm one of these alternatives. Despite these limitations, it is nonetheless clear that the regions homologous to human PAR1 and strata 4 and 5 are still recombining in grey mouse lemur.

Figure 1 with 2 supplements see all
Identification the PAB in seven strepsirrhine species.

(A) Synteny plot of the human and grey mouse lemur X chromosomes. The human X was used to order the grey mouse lemur scaffolds (see Materials and methods). Black dots represent orthologous genes between the human and grey mouse lemur X chromosomes. Human strata number and boundaries follow Skaletsky et al., 2003 and Hughes and Rozen, 2012. Note that old strata have been split into smaller strata in Pandey et al., 2013. Human strata are indicated by different colours. S4 and S5 are in yellow. PARs are in red.(B) M:F read depth ratio along the grey mouse lemur X chromosome. Inferred PAR is shown in red. Regions of elevated M:F coverage ratio (inferred PAR plus other regions in grey) are indicated on panel A plot. (C) M:F read depth ratio for all seven strepsirrhine species. Inferred PARs for both the lemurs and the lorises are shown in red. Details on the PAR and the grey regions of the lorises can be found in Figure 1—figure supplement 1. In all panels, red lines indicate scaffold boundaries. See Materials and methods and Supplementary file 1B for the detailed analysis of the regions with elevated M:F coverage ratio shown in grey.

We repeated the same analysis for the other six species (Figure 1C). For the lemurs, we used the grey mouse lemur reference genome for the mapping because it is the only one available, and for the lorises, we used the northern greater galago reference genome for the same reason (see Maerials and methods and Figure 1—figure supplement 1 for the dot plot with the human X). Some species are quite distantly related to focal species with the reference genome and so mapping was consequently more difficult. This explains why in some cases the M:F depth ratio is more variable. The results of the aye-ayes analyses are especially noisy because of the large phylogenetic distance to the grey mouse lemur (Figure 1C). However, in all seven species studied here, the pattern is very similar (Figure 1C and a zoom on the PABs in Figure 1—figure supplement 2). All studied strepsirrhines harbour a large pseudoautosomal region including the genes that are in PAR1 and strata 4 and 5 in humans (compare Figure 1A and C for lemurs and Figure 1—figure supplement 1 and C for lorises; both lemur and loris PABs correspond to the boundary between human strata 4 and 3). We can therefore conclude that no suppression of recombination between the X and the Y has occurred in strepsirrhines since the origin of the group >65 millions years ago (Figure 2).

Strata formation in primates.

Data on strata in haplorrhines are from Lahn and Page, 1999, Skaletsky et al., 2003, Ross et al., 2005, Hughes and Rozen, 2012, Hughes et al., 2012, Cortez et al., 2014. Data on strepsirrhines are from this study. The phylogenetic tree and divergence times are from Horvath et al., 2008, Pozzi et al., 2014. Drawings of primates were prepared by Philippe Faivre.

It is possible that the M:F read depth approach missed recently evolved strata in strepsirrhines. Recent strata are indeed more difficult to detect with the M:F read depth approach as sex-chromosome divergence can be so low that both X and Y reads map onto the X chromosome and the ratio is close to 1 (Wright et al., 2017). To identify recent strata, we computed the male:female SNP density ratio, which is expected to more effectively detect the PAB when recent strata are present (Vicoso et al., 2013a; Wright et al., 2017). The M:F SNP density ratio is predicted to be one for the PAR, <1 for old strata due to haploidy in males and >1 for recent strata due to accumulation of fixed X-Y differences (Wright et al., 2017). However, our analyses revealed no recent strata in the seven strepsirrhine species studied here (Figure 3).

SNP density analysis.

M:F SNP density ratio (ln scale) for all seven strepsirrhine species (see Materials and methods for details). Dashed lines indicate the mean M:F SNP density across sliding windows of the same size on chromosome 4, the 97.5 and 2.5% quantiles, to show the variation across the autosomes. See legend of Figure 1 for more details.

Our findings are consistent with the hypothesis that recombination suppression between X and Y chromosomes was driven by sexually antagonistic mutations. However, the rate of strata formation is generally low: in primates two strata (4 and 5) were formed in apes and Old World monkeys, one was formed independently in New World monkeys (4’) based on the species studied so far (Hughes et al., 2012; Cortez et al., 2014, and see Figure 2) and our additional data found no new strata formation in strepsirrhines. This observation is consistent with our hypothesis, but could have happened by chance because of a low common rate of strata formation in both suborders. We designed a statistical test to compare the rates of strata formation (expressed in event per My) taking into account the respective divergence times in the haplorrhine and strepsirrhine parts of the phylogenetic tree of the studied species, but this test was only marginally significant (binomial test, p=0.051 see Materials and methods). Because haplorrhines and strepsirrhines have different generation times, comparing rates on a generation-based timescale might however be more relevant. Rescaling time in generations to compare rates of strata formation (expressed in event per million generations) lead to a significantly higher rate in haplorrhines (binomial test, p=0.01 see Materials and methods), consistent with our hypothesis.

We collected phenotypic data from the literature for our set of 13 primate species and confirmed that our sets of strepsirrhine and haplorrhine species differ significantly in sexual dimorphism (teeth and body size, assuming that they reflect the global level of sexual dimorphism in an organism; see Materials and methods and Table 1) but not in sperm competition (testes size, see Materials and methods and Table 1).

Table 1
Measures of sexual dimorphism and other features in the set of studied haplorrhine and strepsirrhine species.
SpeciesMale canine height in mm (sexual selection)Female canine height in mm (sexual selection)Refs*Combined testes mass in g (sperm competition)Male body mass in g (sperm competition)Refs*Male body mass in g (sexual selection)Female body mass in g (sexual selection)Refs*Social and mating systemRefs*
Callithrix jacchus5.084.81[1]1.3320[1]317324[1]Multimale[1]
D. madagascariensisNANANANANANA26212446[8]Multimale[9]
E. rubriventer10.499.98[2]1.762512[1]19801940[1]Monogamous[1]
G. senegalensis4.013.61[2]1.66210[3]227199[2]Unimale/Polygynous[2]
Gorilla gorilla30.2617.4[2]29.6169000[1]17040071500[1]Polygynous[1]
Homo sapiens10.859.97[2]40.566825[1]7210062100[1]Monogamous/Unimale/
Macaca mulatta16.978.13[2]469200[4]110008800[1]Multimale[1]
M. murinus2.072.08[2]2.4960[5]5963[2]Unimale/Polygynous***[2]
N. coucang7.056.8[2]1.21058[6]679626[2]Unimale/Polygynous[2]
O. garnettii6.536.04[2]8.93320[7]794734[2]Unimale/Polygynous[2]
Pan troglodytes21.7215.26[2]128.944670[1]5970045800[1]Multimale[1 , 10]
Pongo pygmaeus2715.95[2]35.374640[1]7850035800[1]Unimale/Multimale[2]
P. simus**5.945.91[2]NANANA25322248[8]Monogamous[2]
  1. *References are [1] Lüpold et al., 2019, [2] Thorén et al., 2006, [3] Gomendio et al., 2011, [4] Harcourt et al., 1995, [5] Lüpold, 2013, [6] Anderson et al., 1999, [7] Dixson and Anderson, 2004, [8] Taylor and Schwitzer, 2011, [9] Mittermeier et al., 2013, [10] Soulsbury, 2010.**or Hapalemur griseus or H. alaotrensis ***or multimale, see Soulsbury, 2010. Male body mass values for sperm competition and sexual selection analyses of the same species may differ because they come from different sources (even when a single reference is mentioned).


Our work shows that, during primate evolution, the PAB has remained unchanged in strepsirrhines while several X-Y recombination suppression events have shortened the PAR in haplorrhines. We interpreted this as a consequence of differences in sexual dimorphism, and therefore sexual conflict, in both groups. However, strepsirrhines and haplorrhines differ in many ways and it is of course possible that other aspect(s) of their biology drove the pattern that we found. Strata formation may be influenced for example by gene flow (Matsumoto et al., 2017) and meiotic drive (Scott and Otto, 2017) as suggested recently. Previous work has shown that the genetic diversity of strepsirrhines is highly variable (e.g. Perry et al., 2012). It is however unknown whether strepsirrhines and haplorrhines exhibit systematic differences in gene flow rates and meiotic drive dynamics. A limit of this work is the use of a qualitative description of sexual dimorphism and not a quantitative one, which we could have compared to the number of strata. Future work could explore strata formation in more species and gain sufficient statistical power to compare the number of strata to phenotypic data on sexual dimorphism in primates using trait-evolution phylogenetic methods, which requires large datasets.

Evidence for the sexually antagonistic mutations hypothesis has been found in other organisms. In guppies, while the Y chromosome exhibits low levels of divergence from the X (Wright et al., 2017; Bergero et al., 2019; Darolti et al., 2019), populations exhibiting stronger sexual dimorphism seem to have a larger non-recombining region (Wright et al., 2017; Wright et al., 2019; Almeida et al., 2020). In the brown alga Ectocarpus, sexual dimorphism is extremely low and as expected sex chromosomes are homomorphic, with a small non-recombining region, despite being very old (Ahmed et al., 2014). It should be noted, however, that other forces might be driving the process of strata formation in some lineages. In ruminants, the PAR seems to have undergone a process of attrition due to accumulation of DNA repeats (Van Laere et al., 2008; Raudsepp and Chowdhary, 2016). In Microbotryum violaceum, strata are found on the mating-type chromosomes despite the fact that this species only has mating types and not sexes, such that sexual antagonism is absent (Branco et al., 2017). Thus, sexually antagonistic mutation may not be a ubiquitous explanation of strata formation in all organisms.

Although sexual dimorphism is generally low in strepsirrhines, there are some differences among species in this lineage, with the genus Eulemur exhibiting the most pronounced sexual dimorphism (Petty and Drea, 2015). In these species, including the red-bellied lemur (E. rubriventer), which was analysed here, males and females exhibit striking sexual dichromatism, that is they differ in pelage colouration (Rakotonirina et al., 2017). The red-bellied lemur did not show more evidence for recombination suppression than the other species studied here. Sexual dichromatism may rely on sexually antagonistic mutations. The antagonism might have been solved not through Y-linkage but instead through sex-biased expression for example (Ellegren and Parsch, 2007; Gazda et al., 2020). Future research could focus on sex-biased expression in strepsirrhines to test these ideas.

Materials and methods

Research plan

Request a detailed protocol

To test whether recombination suppression is less frequent on strepsirrhine sex chromosomes compared to haplorrhines, we selected strepsirrhine species that would maximise the representation of this group’s diversity, and that were also readily accessible. We then sequenced a male and female of each species and mapped the obtained male and female reads to a reference X chromosome. The male to female depth ratio was then computed along the length of the X chromosome and the PAB was identified as the boundary between zones with a ratio of one (indicative of the PAR) and zones with a ratio of 0.5 (indicative of the non-recombining region).


Request a detailed protocol

We selected seven species covering as much phylogenetic diversity of Strepsirrhini as possible (see Supplementary file 1A). Both infra-orders (Lemuriformes and Lorisiformes) are equally represented. A male and a female individual were sampled for all species (except O. garnettii, the northern greater galago, for which sequence data from a female individual were retrieved from NCBI, see Supplementary file 1A). Blood samples of E. rubriventer (red-bellied lemur) and P. simus (greater bamboo lemur) were collected from living animals at Zoo de Lyon in EDTA blood collection tubes to avoid coagulation. Hair samples (with follicles and roots) of the female Daubentonia madagascarensis (aye-aye) were collected from a living animal at Zoo Frankfurt. Samples of M. murinus belong to the Brunoy laboratory (UMR7179, France; agreement E91-114-1 from the Direction Départementale de la Protection des Populations de l’Essonne): the biopsies were obtained from muscle tissues after the animals’ natural death. Tissues samples of a male D. madagascariensis, and samples of G. senegalensis (Senegal bushbaby), Nycticebus coucang (slow loris) and of a male O. garnettii were obtained from the tissues and cryopreserved cell collection of the National Museum of Natural History (MNHN, Paris, see Supplementary file 1A).

DNA extraction and sequencing

Request a detailed protocol

DNA from E. rubriventer, P. simus and female D. madagascariensis were extracted using two different Macherey Nagel kits. Blood samples were treated with NucleoSpin Blood Quickpure kit. Hair samples were treated with NucleoSpin DNA trace kit after a mechanical crushing of hair bulbs. DNA from the tissues and cells samples (for other species) was extracted using the DNeasy Blood and Tissue kit (Qiagen) following the manufacturer’s instructions. DNA was stored at −20° C and sent on dry ice to the sequencing platform.

A genomic DNA library was constructed for each sample using Illumina kits (TruSeq nano LT for Hiseq 2500 and 3000 sequencing). Paired-end sequencing was conducted using an Illumina Hiseq 2500 (2 × 125 bp) or 3000 (2 × 150 bp) with one or two individuals per lane at Genotoul, the INRA sequencing platform in Toulouse. Sequences were all found to be of high quality (using FastQC, https://www.bioinformatics.babraham.ac.uk/projects/fastqc) and without contamination. Consequently, no trimming was done. Sequence data and coverage are shown in Supplementary file 1A.

Chromosome assembly

Request a detailed protocol

Reference X chromosomes were not available for any species and genome assemblies were only available for two species that were (1) closely related to, or the same as the species being studied, and (2) assembled to an extent that it would be possible to construct a de novo X chromosome. These were M. murinus (grey mouse lemur, Mmur_2.0 version from NCBI) and O. garnettii (northern greater galago, OtoGar4 version from NCBI).

De novo X chromosomes were constructed for these species using scaffolds from whole genome assemblies on NCBI, which were selected, ordered and oriented against the human X chromosome. This was achieved using SynMap, an online software pipeline within the CoGe toolkit (Lyons and Freeling, 2008; Lyons et al., 2008) that identified putative homologous genes between potential X scaffolds and the human X chromosome with a blast comparison (Altschul et al., 1990) using the Last algorithm (a variant of Blastz, see Schwartz et al., 2003). An algorithm within the SynMap pipeline then identified a colinear series of homologous genes between potential X scaffolds and the human X chromosome as regions of synteny, and these were arranged in order accordingly. The relative gene order DAGChainer option was used, with a maximum distance of 20 genes between two matches and a minimum of five aligned pairs of genes. The human X chromosome reference was sourced from the GRCh37.p13 Primary Assembly on NCBI (Reference Sequence: NC_000023.10).

As the results of some of the analyses in this study required normalisation using an autosome from the corresponding species, a reference autosome was constructed using the same process. In this case, the human chromosome four was used to construct a de novo chromosome four for M. murinus and O. garnettii, which was selected for its similar size to the X chromosome.

Read mapping

Request a detailed protocol

Male and female reads for each species were aligned separately to their most closely related de novo X chromosome using Bowtie version 2–2.2.7 (Langmead et al., 2009). The reads were then sorted according to their position on the de novo X chromosome using Samtools version 1.3.1 (Li et al., 2009; Li, 2011).

Coverage analysis

Request a detailed protocol

Read depth was calculated for each sex at each position from the mapped reads on the de novo X using Samtools. The coverage for each sex was then normalised by dividing the depth at each position by the mean coverage depth for that species and sex on an autosome (chromosome four). The ratio of normalised male to female coverage was then calculated at each position and the data was summarised as a sliding window average using a window size of 150 kb sliding at increments of 10 kb or larger windows and increments depending on the species. This data manipulation was performed using AWK version 4.1.3.

Analysis of the regions of the strepsirrhine X chromosomes with unusual male:female coverage ratio

Request a detailed protocol

In Figure 1, both lemur and loris X chromosomes exhibit regions with male:female coverage ratio close to 1 (shown in grey) in their X-specific parts, where a ratio of 0.5 is expected. The grey mouse lemur has five such regions, the northern greater galago three. The dot plots of the strepsirrhine and the human X chromosomes (see Figure 1 and Figure 1—figure supplement 1) clearly show that little or no homologous genes are found in those regions, which suggest that they may be homologous to other human chromosomes. This would be consistent with the male:female coverage ratio of 1, typical of autosomal regions, that we found for these regions. To explore this possibility, we extracted the sequences of those regions and performed a tblastn against all the human proteins (human genome version GRCh38). In case of isoforms, the longest protein was kept so that a human gene was present only once. We then filtered the tblastn results by keeping only hits with >80% similarity (based on average nucleotide divergence between lemurs and humans) and e-value <10–9. From those, we kept human proteins covered by hits to >80% using SiLix (Miele et al., 2011). Only proteins matching to no more than one region were kept. The results of the tblastn are shown in the supplementary file 1B.

For all regions except one, most homologs that we identified are from the human autosomes, which confirms our hypothesis. These homologs are mainly a few sources: chromosomes 1, 8 and 12 for regions 46.8–48, 61.5–63.7, 92.7–93.7 and 41.6–44.1 of the grey mouse lemur X chromosome, and chromosomes 12, 13 and 20 for regions 80–84.5, 116–133 and 49.5–68.5 of the northern greater galago X chromosome. These results can be interpreted two ways. One possibility is that the assemblies of the lemur and loris X chromosomes wrongly include autosomal scaffolds. Another possibility is that during the evolution of strepirrhines, some autosomal fragments have been translocated to the PAR, and the assembly failed to order these fragments correctly. Our approach cannot tell apart these possibilities but in all cases, our results suggest that these regions are probably assembly errors.

Changing tblastn outputs filtering did not change qualitatively the results. With lower %identity thresholds, we detected autosomal homologs for region 30.3–33.2 (for example, with %identity >65, we found two proteins from chrom. 1, one from chrom. 2 and 1 from chrom 19).

SNP density analysis

Request a detailed protocol

To detect potential regions that may have stopped recombining between strepsirrhine X and Y chromosomes relatively recently, the difference in male to female SNP density was examined for all species. For each sex of each species, SNPs were called from the mapped reads using Samtools mpileup and then converted to profiles using sam2pro version 0.8 from the mlRho package (Haubold et al., 2010). Specifically, sites with coverage <5 were excluded from the analysis and SNPs were called when a site had a minor allele frequency of 0.3 times the site coverage. The ratio of male to female SNP density was calculated for 600 kb sliding windows at increments of 10 kb. 0.001 was added to allow for a Log transformation and male to female SNP density was calculated at each window as Log(sum male SNPs) – Log(sum female SNPs). This calculation was performed using R version 3.3.2. We also calculated SNP density across an autosome (chromosome four) using the same approach and computed mean male to female SNP density and 97.5% and 2.5% quantiles across all windows.

Statistical test on strata formation

Request a detailed protocol

We partitioned the phylogenetic tree with total branch length Δt into two subtrees with branch lengths Δt1 and Δt2, Δt = Δt1+Δt2. Assuming a constant rate λ for the formation of new evolutionary stratum, the number S of new strata formed during a time interval Δt is Poisson-distributed with parameter λΔt

(1) P(S=k)=(λΔt)keλΔtk!.

On the subtree i, during the time interval Δti we observe the formation of Si new strata. We wanted to contrast the following two hypotheses:

  • H0: Strata accumulated at a common rate λ0 on both parts of the tree.

  • H1: Strata accumulated at different rates λi during the time intervals Δti.

We used the number S1 of strata formed in the time interval Δt1 as the test statistics and compute the conditional probability to observe a larger value given the total number S1 +S2 of strata formed in the time interval Δt1+Δt2 under the null hypothesis H0:

(2) P(S1k1|S1+S2=k1+k2,Δt1,Δt2)=(k1+k2)!(λ0(Δt1+Δt2))k1+k2eλ0(Δt1+Δt2)×j=k1k1+k2λ0K1+k2Δt1jΔt2k1+k2jeλ0(Δt1+Δt2)j!(k1+k2j)!=j=k1k1+k2(k1+k2j)(Δt1Δt1+Δt2)j(Δt2Δt1+Δt2)k1+k2j,

where we recognized the binomial distribution. Note that this probability is independent of the common rate λ0 of the Poisson process. Applying this test is conceptually equivalent to tossing an unbalanced coin k1+k2 times with a probability p = Δt1 / (Δt1+Δt2) to get a head and computing the probability to obtain at least k1 times a head.

The phylogenetic relationships and mean divergence times for the included primate species were recovered from a previously published primate phylogeny and divergence dates (Pozzi et al., 2014). Detailed phylogenetic relationships among strepsirrhine lineages (Horvath et al., 2008) were used to infer phylogenetic relationships in the cases when species in our analysis were not included in this reference study. The divergence times are shown in a see Supplementary file 1C. Generation times in the studied primate species are highly variable, and we are mostly interested in comparing per generation rather than per year rates of strata formation. For this purpose, the branch lengths in the phylogenetic tree needed to be rescaled by the generation times. We used the age at first reproduction as a proxy for generation time following the example of Gaillard et al., 2005. Ages at first reproduction for the extant species in the phylogenetic trees were obtained from Ernest, 2003 and maximum-likelihood estimates of this trait were obtained for internal nodes of the phylogenetic tree with the fastAnc method implemented in phytools (Revell, 2012). This method assumes that the age at first reproduction evolves neutrally according to a Brownian motion model (Felsenstein, 1973; Schluter et al., 1997). The branch lengths of the phylogenetic tree were rescaled by the generation times. In order to take into account variable generation time along a branch, we used the following method: We denoted g the time counted in generations and t the time counted in years along a phylogenetic branch. The instantaneous generation time (expressed in years per generation) along a given branch at any time t is γ(t) = dt / dg. We assume a linear trend for γ(t) between an ancestral node (at t = ta, for which γ(ta) = γa) and a descendant node (at t = td, for which γ(td) = γd). This assumption raised the following ordinary differential equation:


where λ = d - γa) / (td - ta). The general form for the solutions of this equation is


where K is an integration constant. The number of generations elapsed on the branch between times ta and td is thus

(3) g(td)g(ta)=tdtaγdγalnγdγa.

The ages at first reproduction for the extant species and their maximum-likelihood estimates as well as the rescaled branch lengths in the primate phylogeny are shown in Supplementary file 1D.

The haplorrhine lineages in our sample have evolved for Δt1 = 188.52 My (44.23 million generations) during which S1 = 3 new strata were formed. The strepsirrhine lineages evolved for Δt2 = 321.32 My (158.52 million generations) and no new strata was formed (S2 = 0). Comparing the rates of strata formation expressed in number of events per million year lead to a marginally significant p-value (one-tailed binomial test, p=0.051), this trend became significant when considering the rates expressed in number of events per million generations (one-tailed binomial test, p=0.010).

Statistical analysis of phenotypic differences among primates

Request a detailed protocol

All statistical analyses were conducted with the R statistical software (R Development Core Team, 2019). Sexual dimorphism based on body mass (SSD, size-based sexual dimorphism) or on canine length (CSD, canine height based sexual dimorphism) was quantified as the logarithm of the ratio of the male to the female values (for instance, SSD = ln(male body mass/female body mass), Plavcan, 2004). The relative testes mass (RTM) was computed as the residual of the linear regression ln(combinedtestesmass)∼ln(malebodymass).

In a first approach, the phylogenetic architecture underlying the data was ignored and we simply compared the average dimorphism value between the two groups (haplorrhines vs strepsirrhines). In a second stage, we accounted for the underlying phylogenetic architecture using phylogenetic contrasts in a classical phylogenetic generalised least square analysis (see Symonds and Blomberg, 2014). Two evolutionary models were investigated: a simple Brownian motion (BM) and the Ornstein-Uhlenbeck model (OU) that includes stabilizing selection. The results based on the latter (OU) model should however be considered cautiously as this analysis is certainly over-parameterized considering the very small sample size (between n = 11 and n = 13 species). Analyses accounting for phylogenetic architecture in the data used the following specialized R packages:adephylo (Jombart and Dray, 2010), ape (Paradis and Schliep, 2019), geiger (Harmon et al., 2008) and phytools (Revell, 2012).

Sexual dimorphism based on body mass (SSD, mean ± standard error) was 0.378 ± 0.097 in haplorrhines and 0.062 ± 0.017 in strepsirrhines. This difference based on n = 13 observations was statistically significant only when ignoring phylogenetic inertia (p=0.043) but no longer significant when considering phylogenetic inertia with a Brownian motion model (p=0.66). Analysis involving an OU model would lead to a significant difference between the two groups (p=0.043) but this analysis may either be over-parameterized or suffer from the lack of phylogenetic signal in our data as revealed by the low Pagel’s λ <0.001 (not significantly different from 0) estimated in the Brownian motion model. In such a case, non phylogenetically-corrected analyses should be reported (Freckleton, 2009).

Sexual dimorphism based on canine height (CSD) showed the same kind of pattern: the mean is 0.385 ± 0.076 in haplorrhines and 0.045 ± 0.016 in strepsirrhines. This difference based on n = 12 observations is only significant when ignoring the underlying phylogeny (p=0.013) but no longer significant (p=0.39) when phylogeny is accounted for with a Brownian motion model (leading to a non different from 0 estimate of Pagel’s λ). The OU model leads to a significant difference between groups (p=0.013).

Based on our n = 11 observations, the average relative testes mass did not significantly differ between haplorrhines (0.18 ± 0.24) and strepsirrhines (−0.21 ± 0.32). In order to avoid using residuals of a generalised least square model, we also compared testes mass in an analysis of covariance model (see Lemaître et al., 2009, for an example) including the male body mass as a covariate using the following statistical model in R: ln(combinedtestesmass)∼ln(malebodymass)+group. The results were however qualitatively unchanged (the p-value associated with the ‘group’ factor was p=0.4).

Data and code accessibility

Request a detailed protocol

All the data generated in this study is available at NCBI (project # PRJNA482296). Scripts for the entire coverage analysis pipeline (suitable for compute clusters using Torque job scheduling) are available on GitHub (https://github.com/rylanshearn/sex-read-depthShearn, 2018; copy archived at swh:1:rev:0e33f6b8158f4e1385af58117afeb762576cc0fb). 

Data availability

All the data generated in this study is available at NCBI (project # PRJNA482296).

The following data sets were generated
    1. Shearn R
    (2018) NCBI BioProject
    ID PRJNA482296. Strepsirrhini Raw sequence reads.
The following previously published data sets were used
    1. Broad Institute
    (2006) ENA Project
    ID PRJNA16955. Genome sequencing of Otolemur garnettii.
    1. Broad Institute
    (2007) NCBI BioProject
    ID PRJNA285159. Microcebus murinus RefSeq Genome.
    1. Broad Institute
    (2012) NCBI BioProject
    ID PRJNA169348. Otolemur garnettii RefSeq Genome.


  1. Book
    1. Bull JJ
    Evolution of Sex Determining Mechanisms
    The Benjamin/Cummings Publishing Company, Inc.
    1. Charlesworth D
    (2017) Evolution of recombination rates between sex chromosomes
    Philosophical Transactions of the Royal Society B: Biological Sciences 372:20160456.
    1. Felsenstein J
    Maximum-likelihood estimation of evolutionary trees from continuous characters
    American Journal of Human Genetics 25:471–492.
  2. Book
    1. Kappeler PM
    2. van Schaik CP
    (2004) Sexual selection in primates: review and selective preview
    In: Kappeler P, van Schaik C, editors. Sexual Selection in Primates: New and Comparative Perspectives. Cambridge: University Press. pp. 3–23.
  3. Book
    1. Mittermeier RA
    2. Wilson DE
    3. Rylands AB
    Handbook of the Mammals of the World: Primates
  4. Book
    1. Plavcan JM
    (2004) 13• sexual selection, measures of sexual selection, and sexual dimorphism in primates
    In: Kappeler P, Van Schaik C, editors. Sexual Selection in Primates: New and Comparative Perspectives. Cambridge University Press. pp. 230–252.
  5. Software
    1. R Development Core Team
    (2019) R: A Language and Environment for Statistical Computing
    R Foundation for Statistical Computing, Vienna, Austria.
    1. Ross MT
    2. Grafham DV
    3. Coffey AJ
    4. Scherer S
    5. McLay K
    6. Muzny D
    7. Platzer M
    8. Howell GR
    9. Burrows C
    10. Bird CP
    11. Frankish A
    12. Lovell FL
    13. Howe KL
    14. Ashurst JL
    15. Fulton RS
    16. Sudbrak R
    17. Wen G
    18. Jones MC
    19. Hurles ME
    20. Andrews TD
    21. Scott CE
    22. Searle S
    23. Ramser J
    24. Whittaker A
    25. Deadman R
    26. Carter NP
    27. Hunt SE
    28. Chen R
    29. Cree A
    30. Gunaratne P
    31. Havlak P
    32. Hodgson A
    33. Metzker ML
    34. Richards S
    35. Scott G
    36. Steffen D
    37. Sodergren E
    38. Wheeler DA
    39. Worley KC
    40. Ainscough R
    41. Ambrose KD
    42. Ansari-Lari MA
    43. Aradhya S
    44. Ashwell RIS
    45. Babbage AK
    46. Bagguley CL
    47. Ballabio A
    48. Banerjee R
    49. Barker GE
    50. Barlow KF
    51. Barrett IP
    52. Bates KN
    53. Beare DM
    54. Beasley H
    55. Beasley O
    56. Beck A
    57. Bethel G
    58. Blechschmidt K
    59. Brady N
    60. Bray-Allen S
    61. Bridgeman AM
    62. Brown AJ
    63. Brown MJ
    64. Bonnin D
    65. Bruford EA
    66. Buhay C
    67. Burch P
    68. Burford D
    69. Burgess J
    70. Burrill W
    71. Burton J
    72. Bye JM
    73. Carder C
    74. Carrel L
    75. Chako J
    76. Chapman JC
    77. Chavez D
    78. Chen E
    79. Chen G
    80. Chen Y
    81. Chen Z
    82. Chinault C
    83. Ciccodicola A
    84. Clark SY
    85. Clarke G
    86. Clee CM
    87. Clegg S
    88. Clerc-Blankenburg K
    89. Clifford K
    90. Cobley V
    91. Cole CG
    92. Conquer JS
    93. Corby N
    94. Connor RE
    95. David R
    96. Davies J
    97. Davis C
    98. Davis J
    99. Delgado O
    100. DeShazo D
    101. Dhami P
    102. Ding Y
    103. Dinh H
    104. Dodsworth S
    105. Draper H
    106. Dugan-Rocha S
    107. Dunham A
    108. Dunn M
    109. Durbin KJ
    110. Dutta I
    111. Eades T
    112. Ellwood M
    113. Emery-Cohen A
    114. Errington H
    115. Evans KL
    116. Faulkner L
    117. Francis F
    118. Frankland J
    119. Fraser AE
    120. Galgoczy P
    121. Gilbert J
    122. Gill R
    123. Glöckner G
    124. Gregory SG
    125. Gribble S
    126. Griffiths C
    127. Grocock R
    128. Gu Y
    129. Gwilliam R
    130. Hamilton C
    131. Hart EA
    132. Hawes A
    133. Heath PD
    134. Heitmann K
    135. Hennig S
    136. Hernandez J
    137. Hinzmann B
    138. Ho S
    139. Hoffs M
    140. Howden PJ
    141. Huckle EJ
    142. Hume J
    143. Hunt PJ
    144. Hunt AR
    145. Isherwood J
    146. Jacob L
    147. Johnson D
    148. Jones S
    149. de Jong PJ
    150. Joseph SS
    151. Keenan S
    152. Kelly S
    153. Kershaw JK
    154. Khan Z
    155. Kioschis P
    156. Klages S
    157. Knights AJ
    158. Kosiura A
    159. Kovar-Smith C
    160. Laird GK
    161. Langford C
    162. Lawlor S
    163. Leversha M
    164. Lewis L
    165. Liu W
    166. Lloyd C
    167. Lloyd DM
    168. Loulseged H
    169. Loveland JE
    170. Lovell JD
    171. Lozado R
    172. Lu J
    173. Lyne R
    174. Ma J
    175. Maheshwari M
    176. Matthews LH
    177. McDowall J
    178. McLaren S
    179. McMurray A
    180. Meidl P
    181. Meitinger T
    182. Milne S
    183. Miner G
    184. Mistry SL
    185. Morgan M
    186. Morris S
    187. Müller I
    188. Mullikin JC
    189. Nguyen N
    190. Nordsiek G
    191. Nyakatura G
    192. O'Dell CN
    193. Okwuonu G
    194. Palmer S
    195. Pandian R
    196. Parker D
    197. Parrish J
    198. Pasternak S
    199. Patel D
    200. Pearce AV
    201. Pearson DM
    202. Pelan SE
    203. Perez L
    204. Porter KM
    205. Ramsey Y
    206. Reichwald K
    207. Rhodes S
    208. Ridler KA
    209. Schlessinger D
    210. Schueler MG
    211. Sehra HK
    212. Shaw-Smith C
    213. Shen H
    214. Sheridan EM
    215. Shownkeen R
    216. Skuce CD
    217. Smith ML
    218. Sotheran EC
    219. Steingruber HE
    220. Steward CA
    221. Storey R
    222. Swann RM
    223. Swarbreck D
    224. Tabor PE
    225. Taudien S
    226. Taylor T
    227. Teague B
    228. Thomas K
    229. Thorpe A
    230. Timms K
    231. Tracey A
    232. Trevanion S
    233. Tromans AC
    234. d'Urso M
    235. Verduzco D
    236. Villasana D
    237. Waldron L
    238. Wall M
    239. Wang Q
    240. Warren J
    241. Warry GL
    242. Wei X
    243. West A
    244. Whitehead SL
    245. Whiteley MN
    246. Wilkinson JE
    247. Willey DL
    248. Williams G
    249. Williams L
    250. Williamson A
    251. Williamson H
    252. Wilming L
    253. Woodmansey RL
    254. Wray PW
    255. Yen J
    256. Zhang J
    257. Zhou J
    258. Zoghbi H
    259. Zorilla S
    260. Buck D
    261. Reinhardt R
    262. Poustka A
    263. Rosenthal A
    264. Lehrach H
    265. Meindl A
    266. Minx PJ
    267. Hillier LW
    268. Willard HF
    269. Wilson RK
    270. Waterston RH
    271. Rice CM
    272. Vaudin M
    273. Coulson A
    274. Nelson DL
    275. Weinstock G
    276. Sulston JE
    277. Durbin R
    278. Hubbard T
    279. Gibbs RA
    280. Beck S
    281. Rogers J
    282. Bentley DR
    (2005) The DNA sequence of the human X chromosome
    Nature 434:325–337.
  6. Book
    1. Symonds MR
    2. Blomberg SP
    (2014) A primer on phylogenetic generalised least squares
    In: Garamszegi L, editors. Modern Phylogenetic Comparative Methods and Their Application in Evolutionary Biology. Springer. pp. 105–130.

Article and author information

Author details

  1. Rylan Shearn

    Laboratoire Biométrie et Biologie Evolutive, CNRS / Univ. Lyon 1, Villeurbanne, France
    Data curation, Software, Formal analysis, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing
    Competing interests
    No competing interests declared
  2. Alison E Wright

    Department of Animal and Plant Sciences, University of Sheffield, Sheffield, United Kingdom
    Software, Formal analysis, Validation, Investigation, Visualization, Writing - review and editing
    Competing interests
    No competing interests declared
  3. Sylvain Mousset

    1. Laboratoire Biométrie et Biologie Evolutive, CNRS / Univ. Lyon 1, Villeurbanne, France
    2. Faculty of Mathematics, University of Vienna, Vienna, Austria
    Software, Formal analysis, Investigation, Visualization, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
  4. Corinne Régis

    Laboratoire Biométrie et Biologie Evolutive, CNRS / Univ. Lyon 1, Villeurbanne, France
    Competing interests
    No competing interests declared
  5. Simon Penel

    Laboratoire Biométrie et Biologie Evolutive, CNRS / Univ. Lyon 1, Villeurbanne, France
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
  6. Jean-François Lemaitre

    Laboratoire Biométrie et Biologie Evolutive, CNRS / Univ. Lyon 1, Villeurbanne, France
    Software, Formal analysis, Investigation, Methodology
    Competing interests
    No competing interests declared
  7. Guillaume Douay

    Zoo de Lyon, Lyon, France
    Competing interests
    No competing interests declared
  8. Brigitte Crouau-Roy

    Laboratoire Evolution et Diversité Biologique, CNRS / Univ. Toulouse, Toulouse, France
    Conceptualization, Resources, Funding acquisition
    Competing interests
    No competing interests declared
  9. Emilie Lecompte

    Laboratoire Evolution et Diversité Biologique, CNRS / Univ. Toulouse, Toulouse, France
    Conceptualization, Resources, Data curation, Funding acquisition, Investigation, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5711-7395
  10. Gabriel AB Marais

    1. Laboratoire Biométrie et Biologie Evolutive, CNRS / Univ. Lyon 1, Villeurbanne, France
    2. LEAF-Linking Landscape, Environment, Agriculture and Food Dept, Instituto Superior de Agronomia, Universidade de Lisboa, Lisbon, Portugal
    Conceptualization, Supervision, Funding acquisition, Investigation, Visualization, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2134-5967


Agence Nationale de la Recherche (ANR-‐12-‐ BSV7-‐0002-‐04)

  • Gabriel AB Marais

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.


We thank the Muséum National d'Histoire Naturelle de Paris for access to their collections of tissues and cryopreserved cells, Fabienne Aujard (UMR7179 MECADEV – CNRS, MNHN, Brunoy, France) for providing the grey mouse lemur samples and Christina Geiger from Zoo Frankfurt for providing the female aye-ayes samples. This work was performed using the computing facilities of the CC LBBE/PRABI. We thank Bruno Spataro and Stéphane Delmotte for the CC LBBE/PRABI management. We thank Peter Kappeler, Qi Zhou, Mathieu Joron and three anonymous referees for suggestions to improve this manuscript and Alice Baniel, Lounès Chikhi, Laurent Duret and Hervé Philippe for discussions. RS, EL, BCR, and GABM acknowledge the financial support by ANR (grant number ANR-12- BSV7-0002-04).

Version history

  1. Received: October 1, 2020
  2. Accepted: November 16, 2020
  3. Accepted Manuscript published: November 18, 2020 (version 1)
  4. Accepted Manuscript updated: November 19, 2020 (version 2)
  5. Version of Record published: December 4, 2020 (version 3)


© 2020, Shearn et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.


  • 1,434
  • 132
  • 8

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Rylan Shearn
  2. Alison E Wright
  3. Sylvain Mousset
  4. Corinne Régis
  5. Simon Penel
  6. Jean-François Lemaitre
  7. Guillaume Douay
  8. Brigitte Crouau-Roy
  9. Emilie Lecompte
  10. Gabriel AB Marais
Evolutionary stasis of the pseudoautosomal boundary in strepsirrhine primates
eLife 9:e63650.

Share this article


Further reading

    1. Evolutionary Biology
    Maryline Blin, Louis Valay ... Sylvie Rétaux
    Research Article

    Animals are adapted to their natural habitats and lifestyles. Their brains perceive the external world via their sensory systems, compute information together with that of internal states and autonomous activity, and generate appropriate behavioral outputs. However, how do these processes evolve across evolution? Here, focusing on the sense of olfaction, we have studied the evolution in olfactory sensitivity, preferences, and behavioral responses to six different food-related amino acid odors in the two eco-morphs of the fish Astyanax mexicanus. To this end, we have developed a high-throughput behavioral setup and pipeline of quantitative and qualitative behavior analysis, and we have tested 489 six-week-old Astyanax larvae. The blind, dark-adapted morphs of the species showed markedly distinct basal swimming patterns and behavioral responses to odors, higher olfactory sensitivity, and a strong preference for alanine, as compared to their river-dwelling eyed conspecifics. In addition, we discovered that fish have an individual ‘swimming personality’, and that this personality influences their capability to respond efficiently to odors and find the source. Importantly, the personality traits that favored significant responses to odors were different in surface fish and cavefish. Moreover, the responses displayed by second-generation cave × surface F2 hybrids suggested that olfactory-driven behavior and olfactory sensitivity is a quantitative genetic trait. Our findings show that olfactory processing has rapidly evolved in cavefish at several levels: detection threshold, odor preference, and foraging behavior strategy. Cavefish is therefore an outstanding model to understand the genetic, molecular, and neurophysiological basis of sensory specialization in response to environmental change.

    1. Evolutionary Biology
    2. Genetics and Genomics
    Yannick Schäfer, Katja Palitzsch ... Jaanus Suurväli
    Research Article

    Copy number variation in large gene families is well characterized for plant resistance genes, but similar studies are rare in animals. The zebrafish (Danio rerio) has hundreds of NLR immune genes, making this species ideal for studying this phenomenon. By sequencing 93 zebrafish from multiple wild and laboratory populations we identified a total of 1,513 NLRs, many more than the previously known 400. Approximately half of those are present in all wild populations, but only 4% were found in 80% or more of the individual fish. Wild fish have up to two times as many NLRs per individual and up to four times as many NLRs per population than laboratory strains. In contrast to the massive variability of gene copies, nucleotide diversity in zebrafish NLR genes is very low: around half of the copies are monomorphic and the remaining ones have very few polymorphisms, likely a signature of purifying selection.