Evaluation of Four Methods to Identify the Homozygotic Sex Chromosome in Small Populations

BackgroundWhole genomes are commonly assembled into a collection of scaffolds and often lack annotations of autosomes, sex chromosomes, and organelle genomes (i.e., mitochondrial and chloroplast). As these chromosome types differ in effective population size and can have highly disparate evolutionary histories, it is imperative to take this information into account when analysing genomic variation. Here we assessed the accuracy of four methods for identifying the homogametic sex chromosome in a small population using two whole genome sequences (WGS) and 133 RAD sequences of white-tailed eagles (Haliaeetus albicilla): i) difference in read depth per scaffold in a male and a female, ii) heterozygosity per scaffold in a male and a female, iii) mapping to a reference genome of a related species (chicken) with identified sex chromosomes, and iv) analysis of SNP-loadings from a principal components analysis (PCA), based on the low-depth RADseq data. ResultsThe best performing approach was the reference mapping (method iii), which identified 98.12% of the expected homogametic sex chromosome (Z). The read depth per scaffold (method i) identified 86.41% of the homogametic sex chromosome with few false positives. The SNP-loading scores (method iv) found 78.6% of the Z-chromosome and had a false positive discovery rate of more than 10%. The heterozygosity per scaffold (method ii) did not provide clear results due to a lack of diversity in both the Z and autosomal chromosomes, and potential interference from the heterogametic sex chromosome (W). The evaluation of these methods also revealed 10 Mb of likely PAR and gametologous regions.ConclusionIdentification of the homogametic sex chromosome in a small population is best accomplished by reference mapping or examining read depth differences between sexes.