Contrasted sex chromosome evolution in primates with and without sexual dimorphism

: The human sex chromosomes provide a striking example of how sex chromosomes can become modified in gene content and structure (i.e. become heteromorphic). Why sex chromosomes in humans and other lineages can reach such levels of heteromorphy is not clear. One hypothesis 5 posits that mutations with antagonistic effects between males and females can drive recombination suppression and subsequent divergence between X and Y chromosomes. Here, we tested this hypothesis by focusing on strepsirrhine primates (lemurs and lorises) that have much lower sexual dimorphism and opportunities for sexually antagonistic mutations than haplorrhines (apes and monkeys). We sequenced seven strepsirrhine species and identified the 10 pseudoautosomal boundary (PAB), the border between recombining and non-recombining regions on the sex chromosomes. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events have occurred in haplorrhines, suggesting that sexually antagonistic mutations might have driven sex chromosomes heteromorphy in primates.


Abstract:
The human sex chromosomes provide a striking example of how sex chromosomes can become modified in gene content and structure (i.e. become heteromorphic). Why sex chromosomes in humans and other lineages can reach such levels of heteromorphy is not clear. One hypothesis 5 posits that mutations with antagonistic effects between males and females can drive recombination suppression and subsequent divergence between X and Y chromosomes. Here, we tested this hypothesis by focusing on strepsirrhine primates (lemurs and lorises) that have much lower sexual dimorphism and opportunities for sexually antagonistic mutations than haplorrhines (apes and monkeys). We sequenced seven strepsirrhine species and identified the 10 pseudoautosomal boundary (PAB), the border between recombining and non-recombining regions on the sex chromosomes. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events have occurred in haplorrhines, suggesting that sexually antagonistic mutations might have driven sex chromosomes heteromorphy in primates. 15 The human sex chromosomes are strongly heteromorphic as they exhibit extensive differences in size, gene number, DNA repeat abundance and heterochromatin composition 1, 2 . 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 5 two regions called pseudoautosomal regions (PAR1 and 2) that recombine in both males and females. Strikingly, our strongly heteromorphic sex chromosomes originated from a pair of nearly identical autosomes ~150 millions years ago 1, 3, 4 . It all started with the emergence of Sry -the master male-determining gene in mammals-from Sox3 just prior to the divergence of placentals and marsupials 1, 3, 4 . Since then, at several moments throughout evolutionary history, 10 vast regions of the Y chromosome have at once stopped recombining with the X, likely through inversions on the Y, forming the so-called evolutionary strata 3, 5, 6, 7 . 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 8 . In humans, PAR1 is the consequence of that process, while PAR2 is a recent 15 addition 1 .
Such a process of gradual suppression of recombination between sex chromosomes, of PAR reduction, and formation of evolutionary strata has been documented in several animal and plant lineages (e.g 9, 10, 11 ). Why such a process occurred, however, is unclear. It has been proposed that sexually antagonistic mutations may have favoured the suppression of recombination 12, 13 . 20 Theoretical models suggest that if there are male-beneficial/female-detrimental mutations in the PAR, genetically linking those mutations to the Y chromosome, through for example an inversion, will be selected. Some evidence supporting this hypothesis has recently been found in certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which was not this version posted October 16, 2018. . https://doi.org/10.1101/445072 doi: bioRxiv preprint guppies 14 , but in many groups including primates where it has been proposed to explain sex chromosome evolution, evidence is lacking. 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 5 females dominate males, is widespread and likely ancestral 15,16 . FSD is associated with increased testosterone production in females, resulting in the masculinization of females, including aspects of their social behaviour and genitalia 15,16 . Some species also have rather egalitarian social systems 17 . In addition, sexual size dimorphism is virtually absent among strepsirrhines 15,16 . This is in sharp contrast with haplorrhines, where sexual dimorphism is much 10 more pronounced and male-biased 18 . 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 gray mouse lemurs 15 (Microcebus murinus, 19 ) 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 10, 20,21,22 . For autosomes, a M:F depth ratio of 1 is expected as males and females have the same number of autosomes. On the X chromosome, a ratio of 1 should 20 indicate the PAR that is shared among sexes, a ratio of 0.5 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 certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which was not this version posted October 16, 2018. . https://doi.org/10.1101/445072 doi: bioRxiv preprint seven species covering the main strepsirrhine lineages representing 65 My of evolution 23 : four Lemuriformes (aye-ayes, gray mouse lemur, red-bellied lemur, greater bamboo lemur) and three Lorisiformes (northern greater galago, senegal bushbaby, slow loris). The sequencing depth of each sample was between 11.8X and 39.1X (of human genome equivalent) with 78% of the samples being between 20X and 40X, i.e. moderate sequencing depth (Table S1). We then 5 mapped the reads onto publicly available reference genomes of two strepsirrhines (using the human X to scaffold the strepsirrhine X chromosome) and computed a normalized M:F depth ratio to identify the X-specific region and the PAR on the X chromosome (see Methods).

Results
10 Figure 1A-B shows the results for the gray mouse lemur. Using the human X chromosome to order the gray mouse lemur X scaffolds, we found that the scaffolds corresponding to human strata 4, 5 and PAR1 have a M:F depth ratio around 1 (Fig. 1B), indicating that these regions have remained pseudoautosomal in gray mouse lemurs and agreeing with older cytogenetic data 19 . The rest of the gray mouse lemur X is X-specific with a ratio close to 0.5. However, five 15 regions in the X-specific region show an elevated ratio. Detailed analysis of these five regions showed that they are fragments of autosomes (see Text S1). 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. Only an improved assembly of the 20 X chromosome in the gray mouse lemur could confirm one of these alternatives. With the fragmented assembly that is available our approach only reliably identifies the PAB, not the size of the PAR. If some autosomal material were translocated to the PAR, and thus enlarging it, it certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which was not this version posted October 16, 2018. . https://doi.org/10.1101/445072 doi: bioRxiv preprint would not be possible to detect it with our approach. Despite these limitations, it is nonetheless clear that the regions homologous to strata 4, 5 and PAR1 of humans are still recombining in gray mouse lemur and that the PAB is still the same PAB as in the ancestor of all primates. We repeated the same analysis for the other six species (Fig. 1C). For the lemurs, we used the gray mouse lemur reference genome for the mapping because it is the only available one, and for 5 the lorises, we used the northern greater galago reference genome for the same reason (see Methods and Fig. S1 for the dot plot with the human X). Some species are quite distant from the species with the reference genome and 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 species with the reference 10 genome (Fig. 1C). However, in all seven species studied here, the pattern is very similar. All studied strepsirrhines harbor a large pseudoautosomal region including the genes that are in PAR1, strata 4 and 5 in humans. We can therefore conclude that no recombination suppressing event has occurred in strepsirrhines since the origin of the group >65 millions years ago (Fig. 2).
Our findings are consistent with the hypothesis that recombination suppression between X and Y 15 chromosomes is driven by sexually antagonistic mutations. However, the rate of strata formation is generally low and additional strata in strepsirrhines might be absent just by chance. In haplorrhines, at least two strata (4 and 5) were formed in apes and Old World monkeys, and one was formed independently in New World monkeys (4') based on the species studied so far ( 8, 24 and see Fig. 2). We designed a statistical test to compare the rate of strata formation in both 20 suborders using the phylogeny of the species in which the presence of strata was examined. Our test points to a significant difference between both lineages (see Text S2). It is possible that the M:F read depth approach missed recently evolved strata in strepsirrhines. Recent strata are certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which was not this version posted October 16, 2018. . https://doi.org/10.1101/445072 doi: bioRxiv preprint 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 14 . 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 14,21 . The M:F SNP density ratio is predicted to be 1 for the PAR, 0.5 for the old strata (using the same rationale as for the M:F 5 read depth ratio) and to be increased (>1) for the recent strata due to accumulation of fixed X-Y differences 14 . However, our analyses revealed no recent strata in the seven strepsirrhines species studied here ( Figure S2).

10
Evidence for the sexually antagonistic mutations hypothesis has been found in other organisms.
In guppies, populations exhibiting stronger sexual dimorphism seem to have more strata 14 . 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 25 . It should be noted, however, that other forces might be driving the process of strata formation in some 15 lineages. In ruminants, the PAB seems to have moved through a process of erosion due to DNA repeats accumulation 5 . Moreover, in fungal species having mating types instead of sexes, strata have formed when genes important for a given mating type were linked to the mating locus 26 .
Thus, sexually antagonistic mutation may not be a ubiquitous explanation of strata formation in all organisms. 20 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 16 . In these species, including the red-bellied lemur (Eulemur rubriventer), which certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which was not this version posted October 16, 2018. . https://doi.org/10.1101/445072 doi: bioRxiv preprint was analysed here, males and females exhibit striking sexual dichromatism, i.e. they differ in pelage colouration 27 . The red-bellied lemur showed no more evidence for recombination suppression than the other species studied here. Sexual dichromatism may imply sexually antagonistic mutations. The antagonism might have been solved not through Y-linkage but instead through sex-biased expression for example 28 . Similarly, in Emu (Dromaius 5 novaehollandiae), sex chromosomes are homomorphic and sexual dimorphism is thought to be have evolved mainly through sex-biased expression 22 . Future research could focus on sex-biased expression in strepsirrhines to test these ideas.

10
Research plan. To examine whether strepsirrhine sex chromosomes have undergone less recombination suppression 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 15 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).

20
Sampling. We selected seven species covering as much phylogenetic diversity of Strepsirrhini as possible (see Table S1). Both infra-orders (Lemuriformes and Lorisiformes) are equally represented. A male and a female individual were sampled for all species (except Otolemur certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which was not this version posted October 16, 2018. . https://doi.org/10.1101/445072 doi: bioRxiv preprint garnetti, the northern greater galago, for which sequence data from a female individual were retrieved from NCBI, see Table S1). Blood samples of Eulemur rubriventer (red-bellied lemur) and Prolemur simus (greater bamboo lemur) were collected from living animals at Zoo de Lyon on EDTA blood collection tubes to avoid coagulation. Hair samples (with follicles and roots) of  Table S1).  Table S1.
Chromosome assembly. 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 5 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 Microcebus murinus (gray mouse lemur) and Otolemur garnettii (northern greater galago).
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 10 chromosome. This was achieved using SynMap, an online software pipeline within the CoGe toolkit 29, 30 that identified putative homologous genes between potential X scaffolds and the human X chromosome with a blast comparison 31 using the Last algorithm (a variant of Blastz; 32 ). 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 15 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 20 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 certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which was not this version posted October 16, 2018. . https://doi.org/10.1101/445072 doi: bioRxiv preprint four for Microcebus murinus and Otolemur garnettii, which was selected for its similar size to the X chromosome).

Read mapping.
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 33 . The reads were then 5 sorted according to their position on the de novo X chromosome using Samtools version 1.3.1 34,

.
Coverage analysis. 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 10 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. This data manipulation was performed using AWK version 4.1.3. Scripts for the entire coverage analysis pipeline (suitable for compute 15 clusters using Torque job scheduling) are available on GitHub (https://github.com/rylanshearn/sex-read-depth).

SNP density analysis.
To detect potential areas that may have stopped recombining between strepsirrhine X and Y chromosomes relatively recently, the difference in male to female SNP 20 diversity was also 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 36 . Average SNP density was then calculated at each position certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which was not this version posted October 16, 2018. . https://doi.org/10.1101/445072 doi: bioRxiv preprint following 14 as the sum of SNPs divided by four, and one was added to allow for a Log2 transformation. The male to female SNP density was then calculated at each site as Log2(average male SNP density) -Log2(average female SNP density). This calculation was performed using AWK version 4.1.3 and R version 3.4.0. Scripts for the SNP density analysis are also available on the GitHub repository (https://github.com/rylanshearn/sex-read-depth). 5 Data availability: all the data generated in this study is available at NCBI (project # PRJNA482296).

5.
Van      certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.